Hadronic D and Ds Meson Decays 

Anders RycO 

Laboratory of Elementary-Particle Physics, Cornell University, Ithaca, NY 14853 USA 
Alexey A Petro\ll] 

Department of Physics and Astronomy, Wayne State University, Detroit, Ml 48201 and 
Michigan Center for Theoretical Physics, University of Michigan, Ann Arbor, Ml 48109 

We provide a comprehensive review of hadronic decays of D and Ds mesons. We discuss current 
theoretical and experimental challenges and successes in understanding of hadronic transitions 
of those mesons. A brief overview of the theoretical and experimental tools are given before 
discussing the absolute branching fractions for D and Ds mesons. Cabibbo suppressed and rare 
hadronic decays are discussed and compared with theory before discussing our understanding of 
hadronic multibody decays. 
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I. INTRODUCTION 



The discovery of charmed meson states in 1974 sig- 
naled a new era in particle physics. The arrival of the first 
heavy quark has solidified the evidence that the Standard 
Model (SM) provides a correct low-energy description of 
particle physics. Three decades later, the charm quark 
still plays an important role in studies of strong and weak 
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interactions. It also serves as an important tool for ex- 
ploring physics beyond the Standard Model, indirectly 
probing energy scales well above several TeV, which will 
be directly probed by the Large Hadron Collider (LHC). 
In some cases, charm transitions provide possibilities for 
almost background-free studies of low-energy signals of 
new physics (NP) . For example, signals of CP violation in 
the charm system predicted within the Standard Model 
are very small, so any observation of CP violation in 
the current round of experiments would rather unam- 
biguously signal presence of new physics. Charm is also 
rather unique in that it is the only up-type quark that 
can have flavor oscillations. 

A distinctive feature of all charmed hadrons is that 
their masses, 0(2 GeV), place them in the middle of the 
region where non-perturbative hadronic physics is oper- 
ative. While this fact does not markedly affect theoret- 
ical description of leptonic and semileptonic decays of 
charmed hadrons, it poses significant challenges in the 
analyses of their hadronic transitions. There is a great 
deal of optimism, however, that abundant experimental 
data would provide some hints on the structure of charm 
hadronic decays, so those problems will eventually be 
overcome. 

The data on charm transitions originate from several 
different types of experiments. Experiments at e^e~ ma- 
chines operating at the -0(3770) and 0(4140) resonances, 
such as CLEO-c and BES III, have several important ad- 
vantages. First, the final state is extremely simple, be- 
ing essentially just a DD pair. Second, the cross-section 
for charm production is relatively high, (j(D^D^) = 
3.66 ±0.03 ±0.06 nb anAcj{D+D-) = 2.91 ± 0.03 ± 0.05 
nb. In conjunction with low multiplicity of the final state, 
this allows for measurements of absolute branching frac- 
tions for several reference modes. Finally, in those experi- 
ments, the DD pairs are produced in a quantum-coherent 
state, which allows for unique probes of the structure of 
decay amplitudes and phases, as well as novel measure- 
ments of mixing and CP violation. 

The B factory e+e" experiments BABAR and Belle, 
operating at the T(4S) center-of-mass energy, produce 
significant amount of charm data. In fact, at the res- 
onance center-of-mass energy, a{bb) 1.1 nb, while 
a{cc) ~ 1.3 nb. The very large integrated luminosities 
of these experiments have produced large samples of re- 
constructed charm. The higher operating energy makes 
possible the production of charmed baryons. 

Experiments at hadron machines, such as CDF and 
D0, and fixed targcd facilities are plagued by even 
higher backgrounds. However, much higher production 
cross-section, combined with a relatively long lifetime 
of charmed hadrons, provides a possibility to trigger on 
charm decay events with displaced vertices. This tech- 
nique allowed for hadron machines to be major players 
in charm physics. New results from the Large Hadron 
ColHder (LHC) experiments LHCb, ATLAS, and CMS 
will continue to supply us with new data. 

This paper provides a comprehensive review of 



hadronic decays of D and Ds mesons. In this review 
we ado pt the average s perfo rmed by the Particle Data 
Group (|Amsler et all I2008D . Only if there are newer 
measurements that are not included in the review by the 
Particle Data Group we will do our own averaging. 

This review is organized as follows. Section [TT] con- 
tains a brief discussion of the discovery of open charm 
followed in Section IIIII by a discussion of the experimen- 
tal techniques used for studying charm decays. This in- 
cludes a brief discussion of the main experiments that 
have contributed to our understanding of D decays and 
the production mechanisms employed in these studies. 
Final state radiation is discussed in this section as it is 
an important effect in many of the precision measure- 
ments discussed in this review. In Section IIVI the the- 
oretical description of hadronic D decays is provided. 
This includes discussion of SU(3)f flavor symmetry, the 
flavor-flow-diagram approach, and factorization. These 
are common tools used to analyze and interpret hadronic 
D decay data. Sections |V] and IVII discuss the determi- 
nation of the absolute branching fractions for for D and 
Dg decays. Rare and suppressed modes are discussed in 
Section IVIII Multibody decays and Dalitz plot studies 
are discussed in Section IIXI This review concludes in 
Section |X] with a summary and outlook. 



II. DISCOVERY OF OPEN CHARM 



The arrival of t he quark model in 1964 ( Gell-Mannl . 
Il964 IZweiel Il964l ) greatly simplified the description of 
elementary particles. The idea that all observed parti- 
cles are made of the three quarks, u, d, and s, was gain- 
ing acceptance. By the early 1970's, the proton struc- 
ture was probed and the quarks were found to be real 
particles. Further development of perturbative Quantum 
Chromodynamics and the concept of asymptotic freedom 
allowed consistent explanation of those experiments in 
terms of those three quark fiavors. The possible exis- 
tence of a f ourth quark had been theore tically discussed 
in the 60's ( Biorken and Glashowl . [T"964f) . however it was 
not required. 

Hints of the incompleteness of the current picture came 
after experimental observation of rare, electroweak, de- 
cays of kaons. The observed rate for turned 
out to be smaller than predicted. Similarly, the Kg-K^ 
mass difference did not agree with predictions based on 
only having the u, d, and s quarks. To solve those prob- 
lems, Glashow, Iliopoulos, a nd Maiani (GIM) proposed 
an elegant mechanism (Glashow et all \l970) , which in- 
volved adding the forth quark, c. The resulting mech- 
anism not only established the absence of the tree-level 
flavor-changing neutral currents in the Standard Model, 
but also provided for reduced rates for IJ-'^ de- 

cays by requiring cancelations with additional diagrams 
involving intermediate charm quarks. Using the observed 



rate for K'j^ 



and Kg~Kl mass difference, it was 



estimated that the charm quark would have a mass in the 
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FIG. 1 The invariant mass distributions observed by the 
SLAC-LBL experiment for two and four hadrons in the fi- 
nal state, (a) tt+tt" assigning pion mass to all tracks, (b) 
K^Tv^ assigning kaon and pion masses to all tracks, (c) 
K^K~ assigning kaon mass to all tracks, (d) tt'^tt~ weighted 
by Tvn time of flight probability, (e) K^tt^ weighted by Ktv 
time of flight, (f) K+ K~ weighted by KK time of flight, 
(g) 7r+7r~7r^7r~ weighted by An time of flight probability (h) 
K^'K^n^Tv~ weighted by K3n time of flight probability (i) 
K^K ^n^n' weighted by K K-ktt time of flight probability. 
From lGoldhaber et al\ (| 19761 '). 



range 1 to 3 GeV (ICxaillard and Led . ll974l : laaillard et al\ . 
I1975D . The existence of the new quark impHed that it 
would form bound states with its own anti-quark, as well 
as with the lighter quarks, which could be observable ex- 
perimentally. 

These bound states were experimentally discovered in 
Novem ber 1974 by two ind ependent res earch groups at 
SLAG (jAubert et all Il974f ) and BNL (lAugustin et all 
The mass of the observed J /ip resonance of about 
3.1 GeV was in the range where a cc bound state was 
expected. In addition, the very small width, of about 93 
keV, was very different from other high mass resonances 
observed. The interpretation of the J /ip as a cc bound 
state was confirmed when " open charm" states were dis - 
covered a little late r, first the (Goldh aber et a/.l . ll976l ) 
and then the (jPeruzzi et aLl . il976[ ). The first obser- 
vation of the was made in the final states K~'k^ and 
K^TT^TT^Ti^ . The observed invariant mass distributions 
are shown in Fig. [TJ 

After the observation of the D'^ and £)+ mesons it took 



a little longer to establish the D+. There were several 
candidates observed before the Dt ■, originally called the 
F meson, was observed by CLEG ( Ghen e^oH . Il983l ). 

It is also interesting to note that there were hints 
of the existence of open charm states in photoemul- 
sion experiments even before the J/ i /) had been discov- 
ered (|Hoshino et all Il975l: iNiu et all Il97ll ) . 



III. GENERAL REMARKS ON EXPERIMENTAL 
FACILITIES AND TECHNIQUES 

Gharm has been studied in a large number of differ- 
ent experiments. In e~^e~ collisions charm decays have 
been studied from threshold to the Z pole. There has 
also been a number of fixed target experiment, either 
using hadroproduction or photoproduction. The e^e^ 
and fixed target experiments dominate the literature on 
charm meson decays. In addition, there are also studies 
using proton-anti-proton collisions. 

In this section we review some of the basic properties 
of the different types of production mechanisms and the 
experiments used to collect the data. First, e+e" exper- 
iments are discussed and then fixed target. For e^e~ 
experiments, where typically triggering is very open and 
most of the produced events are recorded, we compare 
the luminosity and the produced number of cc events. A 
summary of e"'"e~ experiments is given in Table IH For 
fixed target experiments a similar comparison is made in 
Table |TT] for the number of exclusively reconstructed D 
mesons. At threshold the final state charm mesons are 
produced without any additional hadrons. The GLEO- 
c experiment is described in some more detail as it is 
the experiment operating near threshold with the largest 
data samples to date. At higher e'*'e~ center-of-mass en- 
ergy the charm hadrons are produced either in fragmen- 
tation or in decays of heavier particles such as hadrons 
containing 6-quarks. Last, fixed target experiments are 
reviewed. Fixed target experiments can be categorized as 
photoproduction or hadroproduction experiments based 
on the particle type incident on the target. 



A. Experiments using e^e annihilation near threshold 

At threshold D meson pairs are produced without any 
additional hadrons. This provides the experiments op- 
erating at threshold with a very clean environment for 
studying charm decays. As will be discussed in Sec- 
tion IIII.A.3I the initial electron or positron may radiate 
low energy photons, initial state radiation (ISR), such 
that the total energy of the produced charm hadrons is 
less than the center-of-mass energy in the e"'"e~ initial 
state. 

Experiments that studied charm decays at 
threshol d include the Mark L II, and III e xperi- 
ments (lAbranis et dl Il979at Uugustin et all Il975t 
iBernstein et all Il984^ at SPEAR; EES I, EES II, 
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TABLE I Summary of charm samples produced in e^e colliding beam experiments. 



Experiment 


Year 




J ^ 


Produced Charm 


Mark III 


1982-1988 


3.77 GeV 


7^ ^ 1 

9 pb ^ 


28,000 D"D" 










20,000 D+D- 






4.14 GeV 


6.3 pb-i 




EES 




4.04 GeV 


9 pb-i 


6,000 D+DJ 


EES II 










CLEO-c 


2003-2008 


3.77 GeV 


818 pb^^ 


3.0 X lO'* D°D° 










2.4 X 10" D+D- 






4.17 GeV 


589 pb"^ 


0.58 X lO'' D^D*^ 


CLEO 


1979-1988 


^ 10.5 GeV 


314 pb"' 


0.41 X 10" cc 


CLEO II 


1989-1994 


« 10.5 GeV 


4.7 fb~^ 


6.1 X 10" cc 


CLEO II.V 


1995-1999 


^ 10.5 GeV 


9.1 fb-i 


12 X 10" cc 


CLEO III 


2000-2003 


« 10.5 GeV 


15 fb-i 


19 X 10" cc 


ARGUS 


1982-1992 


« 10.5 GeV 


514 pb-i 


0.67 X 10" cc 


BABAR 


1999-2008 


^ 10.5 GeV 


531 fb-i 


0.69 X 10^ cc 


Belle" 


1999- 


^ 10.5 GeV 


> 700 fb-i 


0.91 X 10^ cc 


MRS 


1982-1986 


29 GeV 


300 pb~' 


52,000 cc 


LEP 


1989-1996 


^ 91 GeV 


4.2 X 10^ Z's 


220,000 cc 



per experiment per experiment 



"As of Sept. 1, 2009 



TABLE II The number of reconstructed charm mesons for 
different fixed target experiments. 



Experiment 


Year 


Events 


Reconstructed 






Recorded/lO" 


Charm Decays 


Photoproduction: 








E691 


1985 


100 


10,000 


E687 


1992 


500 


100,000 


FOCUS (E831) 


1996 


7,000 


1.2 X 10" 


Hadroproduction: 








WA75 


1984 


2 


350 


NA32 


1986 


17 


1,300 


WA82 


1989 


10 


3,000 


E653 


1988 


10 


1,000 


E769 


1988 


500 


4,000 


E791 


1992 


20,000 


200,000 



BES III l|Bai et all [1994 12001"; 'Collaboration', '2009") at 
BEP C and CLEO-c jArtuso et al . 2003; Kubota et d}, 
I1992I: iPeterson et oZI . |2002| ) at CESR-c. For studies 
of D° and D+ decays experiments have run at the 
i/'(3770). The total hadronic cross-section at the •i/'(3770 
resonance has been measured by CLEO-c ( Besson et 1 
[200g) 



hadrons) = (6.38 ± 0.08; 



nb 

-0.30^ 



been measured by CLEO-c (jPobbs et aLU2007^ 
cr(e 



a(e~^e 



D+D-) 



(3.66 ±0.03 ±0.06) nb, 
(2.91 ±0.03 ±0.05) nb. 



The total cross-section for DD production at the -0(3770) 
is cr(e+e- DD) = (6.57 ± 0.04 ± 0.10) nb. This is 
larger than, but consistent with, the inclusive hadronic 
cross-section discussed above. These results indicates 
th at the majority of the 0(377 0) decays t o DD. CLEO- 
c (|Adam et a/.l . l2005[) and BES (|Bai et aLl . [20051 ) have ob- 
served some non-DD decays of the 0(3770). The largest 
of these decays is the radiative transition 0(3770) 7XcO 
with a branching fraction of (0.73 ± 0.09)%. Summing 
the observed branching fractions for non-DD decays 
we obtain 1.4 ± 0.1%, consiste nt with the cro s s-secti on 
measureme nts above. BES ( Ablikim et all . l2006al lbl. 

12008?) has performed direct measurements of the 
cross-section for 0(3770) — > non—DD" final states as 
well as measuremen ts of the DD cross-sections. The 
PDG ( Amsler et all . l2008f ) average these measurements 
and finds that (14.7±3.2)% of 0(3770) resonances decays 
to non-DD final states. This result is inconsistent with 
the CLEO-c results at the 2a level. 

For studies of Ds mesons different e^e" center-of-mass 
energies have been used. The cross-sections for pro- 
du cing -P(s), or D*^^ mes ons, as measured by CLEO- 



c ( Cronin-Hennessv et all . l2008t ). are shown in Fig. O 
BES collected data at 4.03 GeV. At this energy D+D- 
mesons pairs are produced. CLEO-c on the other hand 
ran at a higher energy, about 4.17 GeV. At this en- 
ergy pairs of DfDf^ mesons are produced. The D^ 
meson decays to either Dsj or D^i^^, with branching 
fractio ns of (94.2 ± 0.7)% and (5.8 ± .7)%, respec- 
tively (jAmsler et all . 120081 : lAubert et all . l2005dl ). The 
advantage of the higher energy is the larger cros s-section. 
CLEO-c reports (|C ronm- Hennessv etM . |2008[ ) a cross- 
section of (0.27 ± 0.03) nb at 4.03 GeV for D+Dj pro- 
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FIG. 2 The measured cross-section s for different DD final 
states. From ICronin-Hennessv et ad (|2008l l. 



auction and (0.92 ± 0.05) nb at 4.17 GeV for D;^DJ 
production. For most analyses the larger cross-section 
outweighs the complication of the additional particles in 
the final state. 



1. Quantum coherence 

Threshold production of DD pairs can be explored to 
understand the phase structure of hadronic decay ampli- 
tudes of D° mesons. Here one can use the fact that neu- 
tral charm mesons D'^ and D^ mix. mixing arises 
from electroweak or New Physics |AC| = 2 interactions 
that generate off- diagonal terni s in the neutral D mass 
matr ix (see, e.g. ( Artuso et all 120081 : iBergmann et all 
I2OOOD for more information) 



M- 



A p2 
g2 A 



(1) 



where A parameterizes masses and lifetimes of D^ and 
D^ states and the complex parameters and param- 
eterize contributions from |AC| = 2 interactions. The 
non-diagonal structure of the mixing matrix of Eq. ([1]) 
leads to the (physical) mass eigenstates of a Hamiltonian 
of Eq. Il]) Di and D2 becoming superpositions of the 
flavor eigenstates and I?", 



\D. 



p\D^) ± q\D^ 



(2) 



where |pp 



1. A simplified assumption can 



below CP violation may be neglected. This could 
be justified in the Standard Model by noting that 
CP- violating contributions are always suppressed by 
small values of the third-generation Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements. In general, smallness 
of CP-violating contributions in charm transitions can 
be deduced from tight experimental constraints on CP- 
violating asymmetries (Artuso et al., 2008). In such case 
p — so mass eigenstates also become eigenstates of 
CP, 



1 



D±) = -j= [\D^) ± 



(3) 



It follows then that these CP eigenstates \D±) do not 
evolve with time. Their mass and lifetime differences 
can be observed. 



r 



Ar 



y 



D 



2V 



(4) 



be made that in the studies of strong phases described 



where F — (F+ -I- F_) /2 is the average lifetime of mass 
and CP eigenstates. 

At threshold e+e" experiments, such as BES and 
CLEO-c, D'^D'^ pairs are produced through resonances 
of specific charge conjugation. The D^D° will there- 
fore be in an entangled state with the same quantum 
numbers as the parent resonance. In particular, since 
both mesons are pseudoscalars, charge conjugation reads 
C = (— 1)'^, if the produced resonance has angular mo- 
mentum L. This implies that the quantum mechanical 
state at the time of D^D^ production is 

= i= {\D\\.^)D\\.^)) + C\D\\^^)D\\.^))] . (5) 

where ki and k2 are the momenta of the mesons. Rewrit- 
ing this in terms of the CP basis we arrive at 

*C=+1 = -i={|D+(ki)D+(k2))-|i?-(ki)D_(k2))} 
*C=-1 = -^{|i?-(ki)^+(k2))-f |7^+(ki)7?_(k2))} 

(6) 

Thus in the L = odd; C = — 1 case, which would ap- 
ply to the experimentally important '0(3770) resonance, 
the CP eigenstates of the D mesons are anti-correlated 
while if P = even; C = -1-1 the eigenstates are corre- 
lated. This can happen when D^D^ pair is produced in 
the decays ^'(4140) DDj of the more massive char- 
monium state -0(4140). In either case the CP conserva- 
tion implies that correlation between the eigenstates is 
independent of when they decay. In this way, if D{ki) 
decays to the final state which is also a CP-eigenstate, 
then the CP eigenvalue of the meson D{k2) is therefore 
determined: it is either the same as D{ki) for C = -1-1 
or opposite, as in the case of C = — 1. The use of this 
eigenstate correlation as a tool to inve stigate CP viola- 
tion has been suggested in P-physics ( Falk and Petrovl 
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TABLE III Correlated branching ratios for various processes. 
Correlated results are presented for C — 1 and normalized to 
the multiple of the uncorrelated branching fractions. CP- 
violation is neglected. 

Decay modes Correlated branching fractions 



CLEO-c 



Sotenoid Coll 



K~-K^ vs. K~n 
K~n^ vs. K^n 
K^n^ vs. S± 
K-n+ vs. i± 
S± vs. S± 
S± vs. 5*^ 
S± vs. L± 



{l + Rvjs,? --irsirs+y) 

1 + -Rws ±2rs±y 
1 — VR{y cos 5 + X sin 5) 



4 

l±y 



Barrel 
Cslorimei«r 
Ring Imaging Cherenkov 
/ Detector 

□ rift 
Chamtwr 
Inner Drttt Cltamber / 
'Beampipe 




Barrel Muon 
Chambers 



l200Cl[ ). In charm physics this method of C P -tagging can 
be used to study relative strong phases of Z?°-meson am- 
plitudes. Such measurements are needed for studies of 
D"U°-mixing. 

To illustrate the method, the amplitude for the CP- 
tagged eigenstate decaying to, say, Kir final state can be 
written as 

V2A{D± K-TT+) = AiD'^ K-n+)±A{W K-n+) 

(7) 

which follows from Eq. ([3]). This relation implies that 



1 ±2cosSVR = 2 



K- 



B(DO K-7r+) ' 



(8) 



where i? is a small ratio of doubly-Cabibbo suppressed 
(DCS) decay rate to Cabibbo favored (CF) one (see Sec- 
tion lIVp , and S is the strong phase difference between 
those amplitudes, AiD° K-7t+)/A{D° K-tt+) = 



-/Re 



-iS 



Eq. 



can be used to extract 6 if the CP- 



tagged branching ratio is m easured (jAtwood and Petrovl 
I2Q05I: iGronau et all [200lh . 

The method of quantum correlations can be used to 
stud y the multitude of parameters of D° decay and mix- 
ing ner and "Sm] . l2006HAtwood and Petrovil2005t ). In 
particular, correlated decays of i?-mesons into CP-mixed 
final states (such as K~tt^), CP-specific final states S± 
(such as 5*+ = K^K~ or 5„ = ^Tstt^), or a flavor spe- 
cific semi-leptonic decay into a state containing 
can probe various combinations of mixing and decay pa- 
rameters (see Table imi- Note that P™ = {x^ + y'^)/2 
and Pws = R + VR{y cos S — xsinS) + Rm, and rs — 
y/R cos 5. The quantum-correlated rates are clearly dif- 
ferent from the singly-tagged (ST) rates, i.e. when only 
one of the D^s is reconstructed. For example, the ST rate 
for the wrong-sign (e.g. — > K^it~) decay is given by 
Pws- 

Besides the discussed studies of the phases of hadronic 
decay amplitudes, the results summarized in Table IIIII 
can be used to extract D'^ — D -mixing parameters. The 
current status of charm mix i ng wi ll be discussed else- 
where (see also ( Artuso et aLl . lioOSh ). 



FIG. 3 The CLEO-c detector. The charged particle tracking 
system consists of an inner drift chamber near the interaction 
point and the main drift chamber for the momentum measure- 
ment. Radially outside the main drift chamber is the CLEO-c 
RICH detector for charged hadron identification followed by 
the Csl electromagnetic calorimeter. The instrumented flux 
return for muon detectors is outside the super conducting 
solenoid coil. 



2. Experiments at threshold 

The CLEO-c experiment plays a unique role here as 
it has a very large data sample collected at thresh- 
old. The CLEO-c detector is an evolution of the 
CLEO III detector where the silicon-strip vertex detec- 
tor has been r eplaced with a low - mass inner six-layer 
drift chamber (I Artuso et all 120031: iKubota et aZI . Il992l: 
iPeterson et al\ . 20021 ). The CLEO-c experiment is shown 
schematically in Fig. [31 The wires in the inner drift 
chamber are at a small stereo angle with respect to the 
drift chamber axis. This allows determining the z posi- 
tion of charged particles. The charged particle tracking 
system in CLEO-c also includes the 47-layer main drift 
chamber, operating in a 1.0 T magnetic field along the 
drift chamber axis. The CLEO-c tracking system pro- 
vides a momentum resolution of about 0.6% for tracks 
with a momentum of 1 GeV that traverses all layers of 
the drift chamber. CLEO-c has excellent electromagnetic 
calorimetry from the approximately 7800 CsI(Tl) crystal 
calorimeter. For energies of 1 GeV the calorimeter has 
an energy resolution of about 2%. For energies of 100 
MeV the resolution is about 5%. The excellent energy 
resolution and coverage allow CLEO-c to efficiently re- 
construct TT^ and rj mesons in the 77 final state. The 



tt" mass resolution obtained is about 6 MeV. Charged 
hadrons are identified by a combination of specific ion- 
ization, dE/dx, in the drift chamber for particles with 
momenta below about 700 MeV. For higher momenta, 
where dE/dx is less powerful, CLEO-c uses the RICH 
detector to separate kaons from pions. 

The BESIII ( Collaborationll2009f) detector constitutes 
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FIG. 4 Event display from CLEO-c showing a candidate 
Dt'^Dj' event with D^"'' — » Df'y and both Df candidates 
decaying to 07r'*'. 



reconstructed in the tt+tt" final state. For the tt+tt" 
pairs used to form Kg candidates the usual usual track 
quality criteria are relaxed and no particle identification 
criteria are applied. 

To extract the signal in fully reconstructed hadronic D 
decays it is typically required that the reconstructed D 
candidate energy is consistent with the beam energy, as 
each D in the final state will carry half of the center-of- 
mass energy. Specifically, 



where 



cand 



E 



is the energy of the D candidate. For correctly recon- 
structed D candidates the AE distribution peak at zero. 
The resolution on AE is mode dependent and the actual 
criteria applied varies between different analyses depend- 
ing on the backgrounds and cleanliness of the signal that 
is desired. 

After applying a mode dependent AE selection criteria 
the beam constrained mass is formed 



a substantial upgrade of the earlier BES II detector. 
Among the new features are a 1 T magnetic field gen- 
erate by a superconducting coil, a new drift chamber, 
and a CsI(Tl) dopped electromagnetic calorimeter. The 
time-of-flight system provides tt-K separation at 0.9 GeV 
with a 2cr separation. The operation of the BES experi- 
ment has just started with a first run at the ■0(25'). 



3. Experimental features at threshold 

At threshold D mesons are produced in pairs. A very 
powerful analysis technique involves reconstructing one 
D meson exclusively. This allows us to infer the existence 
of another D mesons in the event. This 'tagging' tech- 
niq ue, or 'double tag' technique, was first used b y MARK 
III (|Adler et aLl . llQSSl : iBaltrusaitis et a/] . ll986l) . but due 
to their relatively small sample of tags the technique was 
of limited use. With much larger samples, and a more 
modern detector, the CLEO-c experiment has made great 
use of this tagging technique. The event environment at 
threshold is very clean. The DD signal is produced with 
no additional hadrons. An example from CLEO-c of a 
fully reconstructed D*^Df is shown in Fig. 21 

Many analyses make use of fully reconstructed D can- 
didates. The D candidates arc built from charged kaons 
and pions, neutral pious, 77 and Kg mesons. CLEO-c 
typically require that kaon and pion candidates are con- 
sistent with charged hadron particle identification based 
on energy loss in the drift chamber and Cherenkov ra- 
diation in the RICH detector 



BC 



E 



beam 



Here the candidate energy has been replaced by the beam 
energy which typically is much better known. 

A typical plot of the Mbc distribution is shown in 
Fig. m The signal yield is determined by fitting the 
Mbc distribution to a background shape plus a signal 
shape. The background shape is due to combinatorial 
backgrounds either from other D decays or from contin- 
uum. Th e background is typic ally fit using an 'ARGUS' 
function (jAlbrecht et all . llQQOf ) 



a(MBc;mo,^,p) = AMbc 1 



The Kg candidates are 



(9) 

This function describes the phase space distribution ex- 
pected near threshold for p = 1/2 and ^ = 0. By allowing 
p and ^ to take on different values a more general function 
which can describe the data better is obtained. 

For the signal shape CLEO-c has used several different 
parameterizations. The most deta i led de scription is that 
used for example in iDobbs et al\ ( 2QQ3) ■ This form in- 
corporates the effects of detector resolution, beam energy 
distribution, initial state radiation, and the line shape of 
the ■0(3770). The beam energy distribution, initial state 
radiation, and the ■0(3770) lineshape control the energy 
of the produced D-mesons. The effect of ISR is to pro- 
duce the ^(3770) with an energy below the nominal e^e~ 
center-of-mass energy. This produces a tail on the high 
side of the Mbc distribution as seen in Fig. [3 The detec- 
tor resolution effects lead to a smearing of the measured 
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FIG. 5 The Msg distribution. The dotted line shows the 
contribution from the ARGUS function that describes the 
combinatorial background. The solid and the dashed lines 
show the contributions to the signal shape for two different 
detector resolution functions. The tail on the high side for 
the signal shape is due to initial state radiation that lowers 
the energy of the produced D mesons. 



D^D^ or D+D events are being fit. The partial widths 
are 



Q'L, 1 



ql l + (r(?oM)^ 



T+{E) = r^s+ 



%M 1 + No)2 

1 + {rq+M) 
9+M + 



and 



(14) 



respectively. In these expressions, is the measured 
width of the V'(3770), Bo{B+) is the branching fraction 
for the decay of the ^(3770) to DW°{D+D-) pairs, 
go (g+) is the momentum of a D^{D^) of energy E/2, 
and qoM{q+M) is the momentum of a D°{D~^) of energy 
CLEO-c used Bo = 0.57 and B+ = 0.43. The pa- 
rameter r is the Blatt-Weisskopf interaction radius taken 
to be r = 12.3 GeV"^ = 2.4 fm. 

The energy distribution of the V'(3770) mesons that 
are produced is obtained by multiplying the e+e~ energy 
distribution f^+e- (E) with the cross-section for ip{3770) 
production. 



momentum. Following iDobbs et al\ (120071 ) . a brief de- 
scription of the lineshape is given below. 

The distribution^ of the energ y of ISR photons is taken 
to be ( Kuraev and Fadinl . llQSSD 



where 



IT 



2 In 



7 ' 



E, 



(10) 



(11) 



At the -0(3770) resonance, /3 w 0.078. The energy distri- 
bution /e+e- (E) of the e"^ and when they collide is 
obtained from an integration of the beam energy spread 
and the ISR photon energy distribution, 

/>oo 

fe+e-{E)^ / h(E^)gE(E + E^)dE^. (12) 
Jo 

CLEO-C has taken the i/'(3770) natural line shape to be 

TiE) 



fBw{E) = 



{E^-Mlf + {M^TT{E)f 



(13) 



where is the mass of the 0(3770). The total width 
Tt{E) is the sum of the partial widths for neutral and 
charged DD pairs, Tt{E) = T ^{E) + T +{E) . The numer- 
ator T[E) is either ro(£') or r+(i?) depending on whether 



^ For simplicity, t hese distribution functions a re not normalized to 
1. The RooFit llVerkerke and Kirkbvl l2003t) fitting package used 
used by CLEO-c takes care of the overall normalization of the 
distribution functions used in fits. 



U{E) - fBw{E) / h{E^)gE{E + E^) dE^. (15) 
Jo 

The 0(3770) energy E is related to g, the magnitude 
(|q|) of the momentum of the produced D and D, by 

E : 



2V? 



Hence, ftpiE) can be transformed into 



a distribution function U]j{q) for the D momentum. 



UD{q) - U{E) 



dE 
dq 



(16) 



The measured D momentum p differs from q due to 
detector resolution. The resolution distribution is de- 
scribed by the sum of three-dimensional Gaussian reso- 
lution functions. Each term in this sum is given by 



S'p(p;q,crp) 



1 



-(p-q)V(2<T^ 



(17) 



(2^)3/V3 

where q is the momentum of the D meson, p is the re- 
constructed momentum, and ap is the momentum res- 
olution, assumed to be the same for both longitudinal 
and transverse components of p relative to the direction 
of q. The D is reconstructed from multiple final-state 
particles, and the vector sum of their momenta tends to 
average out any directional dependence. In the discus- 
sion below, we consider smearing with a single Gaussian, 
for simplicity. 

Since the line shape distribution uo{q) depends only 
on the magnitude q ~ |q| of the D meson momentum, we 
reduce the three-dimensional momentum resolution func- 
tion gp(p; q, Up) to a one-dimensional resolution function 
r(p; g, Up) for the probability distribution of the measured 
value of p = IpI given the produced value of g. This re- 
quires integrating (7p(p; q, (Jp) dpdfl over angles trans- 
verse to q. In this expression, dp dfl is the usual spher- 
ical coordinate volume element and the polar and az- 
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imuthal angles of d^l are relative to the vector q. There- 
fore, 



q \j2-nGp 



(18) 



The distribution of the reconstructed D momentum, 
voip), is then determined by smearing the distribution of 
the true D momentum, U£i{q) of Eq. (|16p . with r{p; q, dp). 



voip) 



r{p;q,(Tp)uD{q)dq 



2711 £ 



r{p-q{E),ap)U{E)dE. 



(19) 



Since the measured value of A/bc is a function of the 
reconstructed momentum p, the distribution function 
w_d(Mbc) of M^c is related to vd{p) by 



wd(Mbc) 

Mbc 
P 



dp 



vd{p) 



dMsc 
r{p;q{E),ap)UiE)dE. 



(20) 



An example of a fit using this form is shown in Fig. [S] 
The tail on the high side is due to ISR and the two com- 
ponents used for the momentum resolution are shown. 

The distribution for double tags, i.e. for Mbc = 
Mbc(-D) and Mbc = Mbc{D), is similar to the form 
developed above for a single Mbc distribution. Since 
both D mesons are produced with the same momentum 
q, Eq. generalizes to the following probability dis- 
tribution for reconstructing the DD pair with measured 
momenta p and p given resolutions ap and (jp, 

Vjj-^{p,p) = J r{p; q, <Tp)r{p] q, ap)uD{q) dq. (21) 

Written in terms of A/bc and Mbc, we have 



w^;o(Mbc,Mbc) 



— , Mbc Mbc 



DD 

oo 



P P 

r{p; q{E),ap)r{p; q{E), dp)f^{E) dE. 



(22) 



In the single tag fits it is hard to separate the effects 
of beam energy smearing and detector resolution. In the 
double tag fits these contributions can be separated as 
the effects of detector resolution is uncorrelated amongst 
the two D candidates, while the beam energy smearing 
is strongly correlated among the two D candidates. 



4. Systematic uncertainties 

Many of the analyses discussed in this review are lim- 
ited by systematic uncertainties. This applies in par- 
ticular to the determination of the Cabibbo favored 



and D"*" reference branching fractions that are discussed 
in Sect. [V] A substantial effort has been put into un- 
derstanding the systematic uncertainties associated with 
track finding, reconstruction, particle identification, 
and TT*^ reconstruction. At the '0(3770) resonance many 
of these uncertainties can be evaluated using hadronic de- 
cays in an event environment very similar to the channels 
studied. This gives confidence in the sometimes small 
systematic uncertainties obtained in these studies. The 
most detailed system atic studies carried out by CLEO-c 
arc described in Ref. ( Dobbs et aZ.I . [20071 ) . As the results 
of these studies are important for many results discussed 
in this review some of these studies are discussed below. 

Track finding has been studied in CLEO-c using a miss- 
ing mass technique where all particles in an event are 
reconstructed except for one particle which we are inter- 
ested in studying. As an example consider the use of the 
kaon in K~tt^ to measure the kaon tracking effi- 

ciency. In this case the opposite l)° in the event would 
be fully reconstructed in some channel and the 7r+ from 
decay looked for. Given the JD" and 7r+ candidates 
the missing mass in the event can be calculated 



M-^ 



(Ptot -Fd - Pother)^ 



(23) 



where pjy is the four-momentum of the reconstructed D, 
Pother is the four-momentum of the other particles that 
were combined with the tag D, in this example the tt"*", 
andptot is the four-momentum of the initial e'*"e~ pair. In 
the missing mass squared calculation, the D momentum 
is rescaled to the momentum magnitude expected from 
the beam energy, but its direction is left unchanged. This 
constraint improves the M^-^^ resolution. 

Candidates where the missing momentum vector fails 
the polar angle requirement |cos0| < 0.9 are rejected. 
This eliminates candidates in which the missing particle 
is expected to be outside the tracking fiducial volume. 
This requirement is tighter than the angular acceptance 
of the CLEO-c detector, | cos 6*1 < 0.93. A tighter re- 
quirement is used to compensate for the resolution in the 
predicted direction of the missing momentum. Later an 
uncertainty is added to compensate for the extrapolation 
to the full tracking fiducial. 

For each missing mass candidate all remaining tracks 
in the event are considered, applying the standard track 
quality criteria. If a track candidate is found such that 
it forms a good D candidate when combined with the 
other particles in the missing mass candidate the missing 
particle is said to have been found. The requirements for 
a good D candidate are |AfBC — Mo\ < 0.01 GeV/c^ and 
\AE\ < 0.05 GeV/c. If a good D candidate is not found 
the missing particle is said not to be found. 

The missing mass candidates are separated into two 
samples; the sample where the missing particle was found 
and the remaining events where the missing particle was 
not found. An example is shown in Fig. [6l The case 
where the missing particle is found corresponds to a 
fully reconstructed 0(3770) event and is very clean. The 
events in this sample are fit to a signal shape using a sum 
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of two Gaussians. A small background component is also 
included in the fit. For the sample where the missing par- 
ticle is not found a clear peak can be seen corresponding 
to the events where we had an inefficiency. In addition to 
this peak there are also substantial backgrounds. These 
backgrounds include semileptonic decays as well as higher 
multiplicity hadronic D decays. These backgrounds are 
parameterized using Monte Carlo simula t ed ev ents. 

As described in detail in iDobbs et all (|2007l ) CLEO-c 
measures the tracking efficiency for both kaons and pions 
in three momentum ranges (0.2 < p < 0.5 GeV/c, 0.5 < 
p < 0.7 GeV/c, and p > 0.7 GeV/c). CLEO-c evaluates 
the tracking efficiency and find agreement between data 
and the Monte Carlo simulation and assigns a per track 
systematic uncertainty of ±0.3% per track for pions. For 
kaons an additional uncertainty of ±0.6% is added due to 
evidence for a tracking efhciency difference between 
and . 
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FIG. 6 Histograms of and fits to M^^^^ distributions from 
D'^ — > K~Tr'^-K^ decays to determine the cliarged pion ef- 
ficiency for > 0.2 GeV/c. Figures (a) and (c) are from 
events in data, and (b) and (d) are from events in Monte Carlo 
simulation. Figures (a) and (b) are from decays in which the 
pion was found, while (c) and (d) are from decays in which 
the pion was not found. The solid curves are fits to the data 
or Monte Carlo sample; the dash ed curves in (c) an d (d) are 
background contributions. From iDobbs et all (|2007l ). 

The Kg TT'^TT" reconstruction efficiency is studied 
in D'^ or D° decays to i^r^Tr+Tr" decays using a technique 
similar to what was used for the tracking efficiencies. One 
tag D is fully reconstructed and two charged pions are 
required to be found. To factor out the track finding 
efficiency and also to reject KItt^tt^ and Kg — > tt^tt^ 
decays it is required that two additional tracks are found 
in the event. These tracks are required to satisfy loose 
consistency requirements with coming from a Kg decay. 
The invariant mass of the two tracks are required to be 
in the range from 0.2 to 0.7 GeV/c^. In addition the dif- 
ference between the missing momentum vector and the 
momentum vector of the sum of the two charged tracks is 




FIG. 7 Histograms of and fits to M^^^^ distributions to deter- 
mine the Kg efficiency. Figures (a) and (c) are from events in 
data, and (b) and (d) are from events in Monte Carlo simula- 
tion. Figures (a) and (b) are from decays in which the Kg was 
found, while (c) and (d) are from decays in which the Kg was 
not found. The background peak and deficit are determined 
by searching for Kg candidates in high and low sidebands of 
the Kg mass. In Figs, (a) and (b), the dashed curves are 
the contributions from fake Kg candidates. In Figs, (c) and 
(d), the dashed curve is the background — a linear function 
with a deficit due to events in which a fake Kg candidate was 
found — and the solid curve is the total fit function including 
the signal peak. The area between the curves is proportional 
to the number of Kg mesons not found. From iDobbs et al\ 
(|2007l ). 



required to be less than 60 MeV/c. Candidates that sat- 
isfy these requirements are searched for a Kg candidate 
found using the standard Kg vertex finder. Similar to 
the tracking studies the candidates are separated into two 
categories; where the Kg was found and where it was not 
found. Compared to the tracking systematics study de- 
scribed above the Kg study is more complicated because 
there are fake Kg candidates from wrong 7r"''7r~ tracks 
in either Kgir^n^ or 7r+7r^7r+7r^ events. This gives rise 
to a 'hole' in the events where the Kg candidate was not 
found because combinatorial background got promoted 
to signal. This is illustrated in Fig. [T] Using this tech- 
nique CLEO-c assigns a systematic uncertainty of ±1.8% 
for the Kg finding efficiency. 

The efficiency for tt*^ —^ 77 reconstruction has been 
studied using a missing mass technique in tpi'^S) — > 
Tpn'^TT^ events recorded at Eqm — "^i/j(2S)- There are 
not any really useful hadronic D decays for tt^ efficiency 
studies. One candidate is the 1)° K^tt^tt'^ decay. If 
the same missing mass technique as described above for 
the charged particle tracking efficiency is applied where 
all particles except for the tt" are reconstructed there is 
a background from the semileptonic decay K^i^v 



which peaks at M, 



2 

miss 



which can not be sepa- 



rated from the signal peaking at Af, 
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miss 



The only 
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mode of some use is £)" — > K'^tt'^. But the relatively small 
branching fraction and reconstruction efficiency makes it 
marginally useful with the data samples currently avail- 
able. 

To study the tt'^ efficiency using the ip{2S) J / ipir'-' tt'^ 
channel CLEO-c reconstruct the J/tp and one 7r°. The 
J/ijj is reconstructed in the e~^e~ and ^i'^ final states. 
The di-lepton invariant mass is required to be within 
50 MeV of the known J/ip mass. The J/ip candidate 
is combined with a tt" candidate and the M^j^^ is cal- 
culated. For signal events, ip{2S) J / i^ir^ ir^ , M^.^^ 
peaks at rn'^o- To remove backgrounds, primarily from 
-0(25*) tAtt", it is required that < 500 MeV/c and 
Pt^o < 500 MeV/c. In addition, to select the kinematic 
region populated by ip{2S) J/t/^tt^tt" it is required that 

(pIo +pLss) - b'o -pLss)V(0.5 GeV2) > 0.10 GeV^ 



1200 



3970407-023 

' — r 



and 



iplo 



-/^miss ) 



- b^o -pLss)V(2 GcV^) < 0.17 GeV^ 



Next a second is looked for in the event. If a second 
7r° exists that satisfies M{J/%Ptt°tt°) - M{J/ip) within 
50 MeV/c^ of the nominal iJj{2S)-J/ip mass difference, 
the 7r° is considered to be found. The distributions of 
the -M^jgg is shown in Fig. [8] for the case (a) where the 
second tt" was found and for the case (b) where the 
second 7r° was not found. Using this study CLEO-c 
finds a difference between the data and the Monte Carlo 
simulation in this sample of 77 = (edata/cMc) ^ 1 = 
(—4.37 ± 0.72 ± 0.41)%. However the average momen- 
tum of the 7r°s in the iIj{2S) — > V'tt^tt" sample is about 
250 MeV/c whereas the typical momentum in D de- 
cays such as K^TT^TT^, K~tt^tt^tt^, and 
D+ /C^Tr+TT^ is higher, typically around 450 MeV/c. 
CLEO-c study the tt" efficiency as a function of the tt*^ 
momentum. The relative data to Monte Carlo efficiency, 
?7, as a function of momentum was fit to a straight line 
in order to extrapolate the efficiency correction to a tt" 
momentum of 450 MeV/c. CLEO-c finds a 3.9% correc- 
tion and assigns a ±2.0% uncertainty to the corrected 
efficiency. 



B. cc production in e+e above threshold 

At energies above charm threshold, charm hadrons are 
produced in fragmentation and are part of a jet, or are 
produced as secondary particles in decays of 6-hadrons. 
The largest charm samples are those produced at the B 
factories at e+e~ center-of-mass energies near 10.58 GeV 
corresponding to the T(45) resonance. The large cross- 
section, about 1.3 nb, combined with the large integrated 
luminosities recorded by CLEO, BABAR, and Belle have 
produced these very large samples. 

At even higher energy, the LEP operated near the Z 
resonance and produced over 4 million Z bosons per ex- 
periment. The jet nature of the events here is more clear 
than at the T(4S'). 
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FIG. 8 Distributions of tt missing mass squared in candi- 
date iP{2S) J/V'ttV events for data (points) and Monte 
Carlo events (histogram). The predicted background level is 
also shown. The vertical arrows demarcate the signal region. 
Events in which the second tt'' was found are shown in (a) 
whereas the events where the second tt" was not found are 
shown in (b). From lDobbs et al.\ (|2007l '). 



The CLEO, BABAR, and Belle experiments were de- 
signed to study B meson decays but they are also well 
suited for studying charm. These experiments all have 
excellent charged particle tracking capabilities and ver- 
tex detectors capable of detecting the separated ver- 
tices from the relatively long lived charm and beauty 
hadrons. All three experiments have CsI(Tl) electro- 
magnetic calorimeters with excellent photon detection 
capabilities and electron identification using E/p. De- 
tection of muons in all three experiments are done using 
an instrumented flux return. Also key for these exper- 
iments is the identification of charged hadrons, particu- 
larly K-n separation. The three experiments chose dif- 
ferent technologies here. BABAR used a DIRC (Detector 
of Internally Reflected Cherenkov light), CLEO-III used 
a RICH (Ring Imagine Cherenkov Detector), and Belle 
uses aerogel Cherenkov counters. All three different types 
of charged hadron particle identification detectors have 
worked well. 

The BABAR and Belle experiments were built 
around an energy asymmetric collision designed to 
allow resolving the time evolution of the produced 
B mesons , as discussed in the BABAR Physics 
Book (Har rison and Ouinnl . ll998[) . The energy asymmet- 
ric collisions are reflected in the design of the detector; 
the interaction point is offset to optimize the acceptance 
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FIG. 9 The FOCUS (E831) spectrometer. 



due to the boost of the collision center-of-mass. 



C. Fixed target experiments 

Charm mesons are sufficiently light that they can be 
produced efficiently in fixed target experiments. The 
main experimental challenge is to separate charm produc- 
tion from the large non-charm rate. The development of 
silicon based tracking detectors enabled experiments to 
effectively identify the long lived charmed hadrons. The 
pioneering Fermilab photoproduction experiment E691 
was the first experiment to produced large samples of re- 
constructed charm hadrons. In this experiment a beam of 
photons with an average energy around 180 GeV was inci- 
dent on a Beryllium target. The cross-section for charm 
production was measured to be about 0.5 /ib. This is 
about 0.5% of the 100 ^b total hadronic cross-section. 
The most powerful tool for identifying the charm signal 
is to make use of the relatively long charm-hadron life- 
times, from (410.1 ± 1.5) fs for the D° to (1040 ± 7) fs 
for the Using the silicon vertex detectors it is pos- 
sible to separate the long lived charm-hadrons from the 
prompt backgrounds. A series of fixed target experiments 
for charm physics are summarized in Table [ill The latest 
of these experiments at Fermilab, FOCUS or E831, re- 
constructed over 1.2 million exclusive charm decays. The 
FOCUS spectrometer is shown in Fig. [51 



D. Final state radiation 

The treatment of final state radiation (FSR) is com- 
mon to many analyses and will be discussed here. In 
many earlier measurements the effects of final state radi- 
ation was often omitted, but as the measurements have 
become increasingly more precise this has become an im- 
portant effect that can not be ignored. In the latest mea- 
surements of the branching fraction for — ^ K~-k'^ the 



size of the radiative correction is larger than the com- 
bined statistical and systematic uncertainties. 

Any reactio n involving charged particles will also ra- 
diate photons ( Bloch and Nordsieckl Il937( ). In fact, an 
arbitrary number of photons will be produced, though 
most of these are very soft. In general, when we dis- 
cuss a branching fraction for a process, like for exam- 
ple B{D^ if~7r+), this includes final states with ad- 
ditional (soft) photons. Experimentally, if photons are 
emitted with an energy that is smaller than the experi- 
mental resolution these events are automatically included 
in the measurement. However, sometimes the photon 
energies are larger, and the energy carried away by the 
photon will make the event fail the selection criteria. In 
order to account for this, and provide a measurement 
of a physically meaningful quantity, experiments sim- 
ulate the effect of final state radiation in their Monte 
Carlo simulations. This has been a common practice 
for semileptonic decays, in particular with electrons in 
the final state for quite some time. For hadronic fi- 
nal states this is not yet universally done. In D decays 
the first experiment that co nsidered FSR corrections was 
CLEO ( Akerib et aZ.I . Il993l ). Today most measurements 
of hadronic D decays include FSR corrections. 

For simulation of final state radiation in hadronic de- 
cays the most commonly used tool is the PHOTOS pack- 
age ( Barberio and Waa Il994l ). In the measurement of 



the -D° K~TT^ branching fraction CLEO-c uses ver- 
sion 2.15 with interference enabled. The effect of interfer- 
ence, here referring to interference between photons ra- 
diated from different charged particles in the final state, 
is important. For the final state D'~' — > K^tt^ the effect 
of including interference changes the fraction of events 
that radiate more than 30 MeV from 2.0% to 2.8%. Ear- 
lier versions of PHOTOS were only able to simulate the 
interference for decays to final state with a particle — anti- 
particle pair. PHOTOS has been compared with calcu- 
lations to higher order in a and found to produce the 
amount of energy radiated very well in semilepton ic de- 
cays of B mesons and decays of r leptons (Richt cr-Wasl . 
I1993D . However, for hadronic final states there is an ad- 
ditional uncertainty introduced by the fact that the fi- 
nal state particles, kaons and pions, are not point like. 
This uncertainty affects in particular higher energy pho- 
tons that probe the structure of the final state particles. 
Higher energy photons could also be radiated directly 
from the quarks; this effect is not included in the simula- 
tion. CLEO-c includes a 30% systematic uncertainty on 
the correction to the branching fraction due to including 
final state radiation. Given the excellent agreement be- 
tween exact calculations and next order calculations in a 
this systematic uncertainty is probably conservative. 

For many earlier measurements it is not always clear 
what was done to correct for the FSR effects. If the effect 
is not included it is hard to correct for it after-the-fact as 
the effect of FSR depends on the selection criteria used 
and how strongly they would reject events with radiation. 
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IV. THEORETICAL DESCRIPTION OF D DECAYS 

Hadronic decays of D-mesons involve transitions of the 
initial-state D-meson into several final state mesons or 
baryons. Thus, they are described by an effective Hamil- 
tonian containing four-quark operators. The theoretical 
description of hadronic decays of charmed mesons is sig- 
nificantly more complicated than leptonic or semileptonic 
ones, although relevant effective Hamiltonians look sim- 
ilar. 

Charmed hadronic decays are usually classified by the 
degree of Cabibbo-Kobayashi-Maskawa (CKM) matrix 
element suppression. Least suppressed, where the quark 
level transitions are c sud are labeled "Cabibbo fa- 
vored" (CF) decays and governed by 



Ti-CF 
O2 



V2 



h.c, 



(24) 
(25) 



where C„(/i) are the Wilson coefficients obtained by per- 
turbative QCD running from Mw scale to the scale /i 
relevant for hadronic decay, and the Latin indices de- 
note quark color. Gp is a Fermi constant, and F^ — 
It, (1 - 75)- 

The "Cabibbo suppressed" (CS) transitions are driven 
by c ^ dud or c ^ sus quark processes. Due to the 
presence of the quark-antiquark pair of the same fiavor 
in the final state, the effective Hamiltonian takes much 
more elaborate form, 

= ^ E Vu,v:AC^{^i)o\ + C2{^i)ol] 



q—s.d 



G 



Ox 

02 



n— 3 



(26) 



where q = d, s, and O3-G are the so-called "pen- 
guin" operators of the type {uc)v-Aj2ni9l)v±A (see, 



e.g. Ref. (|Buccella et aLl . ll995M It is often easy to de- 
note the degree of suppression by powers of the Wolfen- 
stein parameter A — sin^c — Vus — 0.22, there 6c is a 
Cabibbo angle. 

The "Doubly Cabibbo suppressed" (DCS) decay is the 
one in which c dus quark transition drives the decay. 
The effective Hamiltonian for DCS decay can be obtained 
from Eq. (j24p by interchanging s d. 

Calculations of hadronic decay rates governed by these 
transitions are quite complicated and model-dependent. 
Most often, simp l ified assumptions, such as factoriza- 
tion ( Bauer et all Il987t iBuras et~al. I, UHi) are used to 
estimate the needed branching ratios. Some dynam- 
ical approaches, such as QCD su m rules, have been 
used to justify those assumptions ( Blok and Shifmanl . 
I1993D . The main problem with reliable calculations of 



charmed meson decays is that they populate the energy 
range where non-perturbative quark dynamics is active. 
This leads to resonance effect s that affect the phases 
of hadronic decay amplitudes ( Falk et "01. I. I1999D . which 
makes predictions based on factorization quite unreliable. 

Instead of predicting an absolute decay rate, it is of- 
ten useful to obtain relations among several decay rates. 
These relations are helpful when some decay rates in a 
relation are measured, and some are unknown. This al- 
lows for a relation to be used to predict the unknown 
transition rate(s). The relations can be built b ased on 
some symmetries, such as standard flavor SU{3) (jSavaed . 
Il99ll ). or on overcomple t e set of universal qua rk-level am- 
plitudes (|Gronau et all Il994 iRosneH. Il999t) . We shall 
discuss those methods below. 

The partial width for a specific two-body decay of a 
charmed meson depends on both the invariant amplitude 
A and a phase space factor. For a specific two-body decay 
into a PP final state. 



T{D -> PP) = 



P\\A{D^PP)\\ (27) 



SttM 



D 



where |p| is a center-of-mass 3- momentum of each final 
state particle. For a decay into a PV final state, 



T{D PV) 



\A{D~^PV)\' 



(28) 



Note that in the case of PP final state the final state 
mesons are in the S-wave, while in the case of PV final 
state they are in a P-wave. This is why \A{D PP)\ 
has dimension of energy, while |yl(-D Py)\ is dimen- 
sionless. 



A. 5'(7(3)f flavor symmetries 

One popular approach that was adopted for studies 
of hadronic charm decays involves application of approx- 
imate flavor symmetries, such as flavor S'[/(3)_f. This 
approach is based on the fact that the QCD Lagrangian 
acquires that symmetry in the limit where masses of all 
light quarks are the same. The SU (3) p analysis of decay 
amplitudes cannot predict their absolute values. How- 
ever, at least in the symmetry limit, this approach can 
relate transition amplitudes for different decays, which 
could prove quite useful for an experimental analysis. 
One potential difficulty with this approach is related to 
the fact that available experimental data show that fia- 
vor SU{2))f symmetry is broken in charm transitions, so 
symm etry-breaking corrections shou l d be taken in to ac- 
count (|Hinchliffe and Kaedind . Il996t ISavagd . ll99lD . 

In the fiavor-symmetry approach all particles are de- 
noted by their SU{i)F representations. Charm quark 
transforms as singlet under fiavor S'C/(3). The funda- 
mental representation of SU{i)F is a triplet, 3, so the 
light quarks w, d, and s belong to this representation 
with (1,2,3) = (u,d,s). The operator Di that creates a 
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D-meson is of the form cu, so it also transforms in the 
fundamental representation of SU(3). In the hadronic 
decay of a charm meson the final state mesons are made 
of M, d, and s quarks, so they either form an octet 8 rep- 
resentation of SU{3)f (pseudoscalars tt*, 7r°, K^, K'^, 
K , and vectors p"^ , p° , K*^ , K*'^, K* , ujs), e.g. 



V2' 



K- 



(29) 



'78 



or an SU{2>)f singlet (771 and loi). The physical states 
77, 77', (j), and m are linear combinations of 7718 and wi.s 
states respectively. 

The AC = — 1 part of the weak Hamiltonian has the 
flavor structure {qic){qjqk) (see Eq. ([M]) ). so its matrix 
representation is written with a fundamental index and 
two antifundamentals, H^j^ . This operator is a sum of 
irreducible representations contained in the product 3 x 
3x3 = 15 + 6 + 3-1-3. In the limit in which the third 
generation is neglected, H^^ is traceless, so only the 15 
(symmetric on i and j) and 6 (antisymmetric on i and 
j) representations appear. That is, the AC = —1 part of 
Hw may be decomposed as ^{Oj^ + Oe), where 



{sc){ud) + (wc)(sc?) + si{dc){ud) 
+si{uc){dd) — si{sc){us) — si{uc){ss) 
— sl{dc){us) — Si{uc){ds) , 
Oq — {sc){ud) — {uc){sd) + si{dc){ud) 

—si{uc){dd) — si{sc){us) + si{uc){ss) 
-si (dc) {us) + sj (uc) (ds) , (30) 

and si — sinOc ~ 0.22. The matrix representations 
i?(15)^'' and H{6y^ have nonzero elements 



7?(T5)^ : 




= i/|i = l, 
= Hi' - -si , 


Hi- 
Hi- 


- Hf = 

- Hi' - 


si , 




iJf 
HP 


= -i7f = 1, 
= -Hi' = -s, , 


Hi- 
Hl- 


= -H-' 
= -H-' 


= Sl 



(31) 

In the SU{3)f limit the effective Hamiltonian for the 
hadronic decays to two pseudoscalars D PP can be 
written as 

^effsj/O) = aj^D,Hil5y^P^Pt + bj^D,PiH{15y^P^' 

+ c^DM^riPlPi (32) 

There are a number of amplitude relations that can be 
obtained from Eq. ((32| . In particular, it can be seen that 
it implies that |A£)0^^+^- | — \ A]jo^^+^- \. In practice, 
the corresponding branching fractions differ by a factor of 
three (see Table IXXl below) . Clearly, SU{3)f symmetry 
is broken in Z?-decays. 

A consistent approach should then include SU{3)f- 
breaking corrections, which could consistently be in- 
cluded in the analysis. For example, one could assume 



that SU{3)f breaking is proportional to light quark 
masses. In this case, it can be included in the analy- 
sis as a perturbation that transforms as 8 + 1, as the 
quark mass operator belongs to the matrix representa- 
tion M* = diag(77i„, 77i£i, 771s), which is an 8. Note that 
the SU{3)f breaking term that transforms as a triplet 
3 also breaks isospin, so it is usually neglected in all 
analyses. A complete analysis with broken SUi3) F is 
possible (|Hinchliffe and Kaedin}d . ll996HSavagel . ll99l[) . al- 
though is not quite useful due to a large number of un- 
known amplitudes. 

In some cases one does not need to employ the full 
formalism of SU{3)f, but only rely on its subgroups. 
An example of such subgroup is isospin. Isospin rela- 
tions among decay amplitudes are much more robust, as 
isospin breaking is believed to be quite small in charm 
decays. For example, the di-pion modes, —>■ tt+tt*^, 
TT+TT^ and D'-' tt^tt" are related by two isospin 
amplitudes Aq and A2 corresponding, respectively, to the 
5-wave di-pion isospin 7 = and 1 = 2 states produced 



^00 



A2. 
Ao 



A+- 

'a. 



Ao 



A. 



(33) 



Some conclusions about strong interaction dynamics in 
D-meson decays can be reached by extracting these am- 
plitudes from experimental information. The phases of 
amplitudes in Eq. p3p give an indication of the size of 
strong interactions among decay products in those de- 
cays. Following the p rocedure out l ined i n ( Selen et~ai\ . 
Il993h . CLEO obtains (jRubin et al\ . 120061 ) from theh re- 
sults IA2M0I = 0.420 + 0.014 + 0.01 and arg(^2Mo) = 
(86.4 + 2.8 + 3.3)°. As one can see, the phase is rather 
large. It is thus clear that final state interactions play an 
important role in D-decays. 

Other subgroups of the SU{3)f also offer useful pre- 
dictions. For example, the [/-spin, a symmetry of the 
Lagrangian with respect to s — > d quark interchange, 
can be employed to obtain several useful relations. For 
example, for the decays of meson into final states 
containing Af*^ = tt'^, 77, and 77', one can obtain 



^ = - tan-" Be- 



(34) 



Equation derives from the following argument. The 
initial state, D° contains c and u quarks, and so is a U- 
spin singlet. The CF transition c — > sud and DCS tran- 
sition c — !■ dus produce U = I finals states with opposite 
C/3 = 1 in the decays of meson. The final state me- 
son M'^ form a linear combination of [/-spin singlet and 
triplet states, while neutral kaons are U = 1 , U3 = +1 
states. Thus, [/-spin triplet part of Al'^ cannot be pro- 
duced, as it leads to the U = 2 final state. Thus, only 
the singlet part of can contribute to the transition, 
which leads to Eq. 
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B. Flavor-flow (topological) diagram approach 

Another useful approach to tackle hadronic decays of 
charmed mesons, equivalent to the SU{3)f amplitude 
method described above, is the flavor-flow (or topological 
SU{3) approach ), which involves an overcomplete set of 
quark diagrams ( Gronau et aLl . [l99l : lRosne l ll999ri . The 
application of this method to D-decays can even prove 
advantageous compared to flavor SU{3) approach, as the 
number of unknown amplitudes grows rapidly if SU{3) p- 
breaking is taken into account. 

In the topological flavor-flow approach each decay am- 
plitude is parametrized according to the topology of 
Feynman diagrams (see Fig. fTU]) : a color- favored tree am- 
plitude (usually denoted by T), a color-suppressed tree 
amplitude (C), an exchange amplitude (E), and an an- 
nihilation amplitude (A). This set of amplitudes is suf- 
ficient for description of CF and DCS de cays. For SCS 
decay s other amplitudes must be added (IChiang et all . 

horn . 

In order to describe charm meson decays in terms of 
these amplitudes, it is convenient to decompose initial 
and final states according to their isospin structure. For 
instance, in the notation of ( Rosnei, [199^ , the following 
phase conventions are used: 

1. Charmed mesons: D'^ = —cu, Z)+ = cd, and Dg = 
cs. 

2. Pseudoscalar mesons: 7r+ = ud, 7r° = 
{uu — dd) /v^, TT^ = —dm,, = us, — ds, 

TC^ — sd, = —su, rj = (ss — — dd) /VS, 

and rj' = (uu + dd — 2ss) / \/6. 

3. Vector mesons: p+ = ud, p° — (uu — dd) /\/2, 

-du, cj" = (uu -I- dd) /V^, K*+ = us, 

--. — >^o — 

K = ds, K = sd, K*^ = —su, and (f> = ss. 

As with the SU{3)f approach, this method does not pro- 
vide absolute predictions for the branching fractions in D- 
meson decays. However, it provides relations among sev- 
eral decay amplitudes by matching the quark-level "fla- 
vor topology" graphs with the final states defined above. 
For example, a DCS transition 1)° K'^tt~ can proceed 
via a tree- level amplitude T(c uld) and an exchange 
amplitude E{cu 'sd). Matching those with the initial 
state meson 13° = —cu and final state mesons — us 
and 7r~ = —du, one obtains the following amplitude re- 
lation, 

A(I30 ^ K+n-) =T + E= ^KdK* + S) , (35) 

where we use calligraphic notation for the amplitudes 
with Gf/V^ and CKM-factors removed. Similarly, for 
other transitions one obtains 



A(i?° ^ i^V") = -^(C-E) 



V2V2' 



V^sV:, (C" - £■"), 



-L^K.k;(C-.), 



A{D+ K%+) 



C + A 



Gp 

7f 



A{D+ ^ k\+) = T + G= ^VudV:, (T + C), 



A{D" K%) 



1 



-.C 



1 Gi 



VusV:^ C", (36) 



and so on. Note that in Eq. ((36)) we denoted DCS am- 
plitudes with double primes. Singly-Cabibbo-suppressed 
amplitudes are conventionally denoted by a single prime. 
CF amplitudes can be related to SCS and DCS ampli- 
tudes by proper scaling with ianOc. We shall give par- 
ticular examples below. 

One reason for the employed phase convention is a re- 
quirement that SU{2i)f sum rules are satisfied. For ex- 
ample, for transitions 13+ K^n^ , K^V, and 
Z3+ K~^ri' , a sum rule 

3\/2A(i^+7r") + 4:V3A{K+rj) + V6A{K+r]') = (37) 

can be written. With the flavor-flow parameterization. 



AiD+^K^TT^) = —(T-A) 
' \/3 



A{D'' 



(38) 



^ (T + 3A) 
V6 



the above sum rule gives 3(T - A) - AT + {T + 3A) ^ 0. 

Thus, provided that a sufficient number of decay modes 
is measured, one can predict both branching fractions 
and amplitude phases for a number of transitions. Still, 
no prediction for absolute branching ratios are possible 
in this approach. 



C. Factorization ansatz 

The simplest way to estimate an absolute decay rate 
of a charmed meson is to employ a factorization ansatz. 
This ansatz implies that the amplitude for the hadronic 
transition can be written as a product of known form- 
factors. Schematically, 

AiDg^MiAh) - {Mi,M2\n\Dg) (39) 
~ (Mil (ukT'^qk) |0) X (M2I (g.r^c) |D,) 

This is a clear simplification, as the first non-perturbative 
parameter (Mi| (ukT^qk) |0) can be written in terms of a 
meson decay constant /mi , 



{M,\urj5q\0)^tfMA,, 



(40) 
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which parameterizes the amplitude of probabihty for 
quarks to "find each other" in a fight mesons and can 
be measured in leptonic decays of Mi, 



r(Mi ^ Iv) = ^flj^mlmM, ( 1 



'Ml 



IK, 



(41) 

wfiere mj^i-^ is the Mi mass, is tfie mass of the final 
state lepton, and \Vuq \ is the CKM matrix element associ- 
ated with the q —* u transition. The decay constants can 
also be computed in lattice gauge t heories or usin g other 
non-perturbative approaches (see ( Artuso et all |2008| ) 
for review). 

The second non-perturbative parameter, 
(M2I (q^Tf^Ci) \Dq), is related to form-factors that 
can be extracted from semileptonic Dq decays. 



dT{D A'hev^ 



247r3 



iPM.n/+(9')r (42) 



where pAf^ is the hadron 3-momentum in the D rest 
frame. 

Theoretical parameterizations of semileptonic decays 
involve two non-perturbative quantities parameterizing 
the matrix element of a single hadronic current. Tradi- 
tionally, the hadronic matrix elements for transitions to 
pseudoscalar hadrons are described in terms of two form 
factors, f+{q^) and f-{q^), 



{M2\qT^'c\D) 



(43) 



where P = pjj + a-nd q = pjj — ■ An alternative 
parameterization is often used, 

+ "^"r''^ g^/cf-^'^-(g^), (44) 

q^ 



with/--^^'^^(g^) = / 



+ 



{q^)q^/{m 



D 



m\j^). Form factors have been evaluated at specific q 
points in a variety of phenomenological models, where 
the sh ape is typically assu med from some model argu- 
ments (| Artuso aZ.l . l2008f) . 

Clearly, naive factorization of Eq. ((39|) . while conve- 
nient, cannot be correct, as it assumes that scale and 
scheme dependence of a product of quark bilinears is the 
same as that of a four-fermion operator, which it is not. 
The situation can in principle be corrected, at least in the 
heavy-quark limit. In _B-dccays, a QCD factorization for- 
mula has been writte n that takes into acc ount perturba- 
tive QCD corrections ( Beneke et a^.l . [l999f) . It is however 
not clear that this approach is applicable to charm de- 
cays, as charm quark might be too light for this approach 
to be applicable. Nevertheless, even naive factorization 
provides a convenient way to estimate D-meson decay 
rates. 

Besides decay amplitudes for £)-mesons, which can be 
computed using the factorization arguments above, both 



flavor-flow and SU (3) p amplitudes can also be estimated. 
For example, contrary to the relation Eq. p4|) . the cor- 
responding relation for charged Z)-meson decays, 



^ ,^ C" + A" C + {C2lC^)£ 
tan t^C" — 



c + r 



c + r 



A{D+ K"7r+) 

AiD+ ^lfTT+) 

(45) 

cannot be flxed by symmetry arguments alone. However, 
the factorization approach can be used to estimate this 
ratio. In particular. 



r = Mr 



rrijj — m 



1)1^-'' [ml) ai, 
C = fK{nil~ml)f^^^{ml) a^, 



T" 
C" 



fn {mil 



.J f^-^ml) ai, 



)/f-^(m^) a. 



(46) 



where ai^2 — Ci,2 + C'2,i/A^c- Note that some analyses 
employ ai^2 — ^iVi which are fitted from the data and 
treated as universal fit parameters. This way of calculat- 
ing charm hadronic decay matrix elements is sometimes 
called "modified factorization" approach. The argument 
for doing this is an attempt to include unknown non- 
perturbative corrections to Eq. . While this approach 
defines a convenient model to deal with hadronic decays, 
there is no reason to believe that soft contributions are 
universal in all transitions. 

Calculations of £ and A amplitudes in factorization 
are m uch more complicated. It has been argued I Gaol . 
|2007[ ) that they can be estima ted using methods s imilar 
to those employed in B-decays ( Beneke et aLl . [T999[ ). Nu- 
merically, the calculation of the ratio of Eq. (|45l) amounts 
to 



tan^ Ocrse''''' 



(47) 



with w 1.521 and 0,. w 103° for C2/C1 w -0.5. This 
ratio will be used to estimate decay asymmetries with 
kaons later in this paper. 



V. CABIBBO FAVORED AND D+ DECAYS AND 
REFERENCE BRANCHING FRACTIONS 

The absolute branching fractions for decays of the 
ground state charmed mesons are important as they 
are used to normalize many B and D meson decays. 
For e xample, the determination of \Vcb\ from B — > 
D*£iy (Ric hman and Burchat. 1995) depends directly on 
the determination of the D branching fractions used to 
reconstruct the final state. 

To measure the absolute branching fractions we need to 
have a mechanism to determine the number of D mesons 
produced. As the cross-sections for producing D mesons 
are not directly calculable we have to count the D mesons 
in the data sample. Broadly speaking there are two meth- 
ods employed for this D counting. At threshold MARK 
III and CLEO-c have used a tagging technique described 
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FIG. 10 Basic topological amplitudes for D-meson decays. 
Top row: tree T^'' and color-suppressed C^'\ middle row: 
weak annihilation and weak exchange E^'\ bottom row: 
singlet weak exchange SE' . CSC amplitudes are usually de- 
noted by primes. 



in Sect. IIII.Al where one D meson is fully reconstructed 
and tag the existence of another D in the event. At higher 
energies the presence of a D*'^ meson can be tagged us- 
ing the slow pion in the -D*+ — > Z3°7r+ decay. The slow 
pion in this decay is often denoted tTs- This 'slow pion' 
tagging technique has been used by several experiments 
including CLEO and ALEPH to count the number of 
Z3*+ D^TT"*" decays in charm jets produced in e+e^ 
collisions. A variation of this idea has been used by AR- 
GUS, CLEO, and BABAR where D* mesons produced 



in semileptonic B decays, ^ D 



*+ 1 



are tagged by 



the presence of a slow pion and a lepton. These different 
techniques are discussed in this section. 

Before the CLEO-c measurement of the 
K^TT^TT^ branching fraction using tagging as described 
in Sect. IV. CI there wa s a statistics limited study by 
MARK III ( Adler et ai . 1988) and model dependent 
analyses. CLEO ( Balest et ali . 11993) reconstructed the 



two decay chains D*'^ D'^tt'^ 



K- 



and 



D* 



K- 



This allowed CLEO 



to measure the ratio of produced K n 
nal states, which can be expressed as 

Nk^^ _ B{D*+ -> D+7T°)BiD+ K-n+TT+)e{Kn7r) 
Nk-k ^ B{D*+ -> DOtt+)B{DO K^Tr+)e{KTr) ' 

where e{KTnr) and e{K'K) are the efficiencies for recon- 
structing the £)+ K^n^TT^ and K^tt^ final 
states, respectively, including the _D*+. To extract the 
D+ K~TT^TT~^ branching fraction CLEO used the mea- 
sured -D" K^TT^ branching fraction and the ratio 

B{D*+ D+TT°) 
B{D*+ D"tt+)' 

The deterni i nation of this ra tio is discussed by 

iBartelt et all (|l998l ): iButler et a/1 (|l992l ) and ultimately 



relies on isospin conservation. Though the errors are ex- 
pected to be small they are hard to quantify. 



A. Absolute D*^ branching fractions using slow pion tagging 

The method of tagging — > £'°7r+ decays in jets 
produced in e"''e" — > cc interactions by the presence of a 
slow pion from the D* decay is sometimes referred to as 
the HRS technique after the first experiment that used 
this method. As the Q value of the D*'^ D'^tt'^ decay 
is only about 5 MeV and the produced pion has a mo- 
mentum of only 39 MeV in the D* restframe it can at 
most contribute this amount to the transverse momen- 
tum with respect to the thrust axis. Experimentally, the 
slow pion from the D*~^ decay closely follows the origi- 
nal D* direction. Due to the soft track associated with 
this decay, the pion tends to bend out from the jet in the 
magnetic field of the trac king system. 

The HRS experiment (jAba chi et ai\ . Il988l ) used 300 
pb~^ of data collected at Ecm = 29 GeV. For candidate 
slow pions the transverse momentum, px, is calculated 
with respect to the thrust axis determined from the par- 
ticles in the opposite hemisphere with respect to the slow 
pion candidate under consideration. The choice of using 
only tracks in the opposite hemisphere for the calcula- 
tion of the thrust axis is to avoid any possible bias due 
to the decay of the D meson. In Fig. [Tl] the distribu- 
tion is shown in two ranges of the fractional momentum 
xp = 2p||/i?cm of the slow pion, where py is the com- 
ponent of the slow pion momentum that is parallel to 
the thrust axis. In the low fractional momentum range 
(0.03 < xp < 0.06) a clear excess is seen at very low 
values of the transverse momentum due to slow pions 
from D*+ — !■ D'^tt'^ decays. This excess is not present in 
the higher xp range as slow pions from D*'^ decays do 
not populate this range. The HRS collaboration use the 
excess at low to determine that they had 1584 ± 110 
£)*+ £)°7r+ decays in their sample. Next a is recon- 
structed in the Z?" — > K^tt^ channel. The 1)° candidate 



is combined with the slow pion and the mass difference 
Mx-rnis ~ Mkt; is required to be in the range 0.143 to 
0.148 GeV/c^. The yield is determined by fitting the 
Mk-k mass distribution. A total of 56 ± 9 events were ob- 
served. The efficiency for finding the K-k pair, given that 
the TTs is found, is determined to be 79% giving a branch- 
ing fraction of B{D° K-n+) = (4.5 ± 0.8 ± 0.5)%. 
The largest systematic uncertainty quoted is bias due to 
event selection criteria. This uncertainty is evaluated by 
changing the event selection criteria to remove the thrust 
and coUinearity criteria used. The analysis was limited 
by statistics. 

The same technique as pioneered abov e by the HRS 
collaboration has been us ed by A LEPH (iBarate et all . 
119971 : iDecamp aLl.ll99lh. CLEO dAkerib et a^.l ■ ll993^ ■ 
and ARGUS (|Albrecht et all |l994bD. ALEPH used a 
sample of e+e" data collected from 1991 to 1994 at LEP 
near the Z pole. CLEO and ARGUS used samples of 1.79 
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FIG. 11 The pI distribution for data from the HRS collab- 
oration. In (a) the fractional slow pion momentum in the 
range 0.03 < xp < 0.06 is shown while in (b) the fractional 
momentum range 0.07 < xp < 0.1 is shown. In the low mo- 
mentum range where we expect slow pi ons from D* decays a 
clear excess at very low pt is seen. From lAbachi et all (| 19881 ). 



fb^^ and 355 pb~^ respectively of e+e^ data collected 
near the T(4iS') resonance. 

ALEPH followed the HRS approach closely. They an- 
alyzed the data in six ranges of the slow pion momen- 
tum, from 1.0 to 4.0 GeV/c. The transverse momen- 
tum squared distributions in the six momentum bins are 
shown in Fig. [T^l A D'^ K^n^ candidate is searched 
for in events with a slow pion, and candidates where 
0.1435 < Mk^tt, - Mktt < 0.1475 GeV are accepted. 
In Table IIVI the yields and branching fractions from the 
ALEPH analysis are summarized. The results from the 
different momentum bins are combined, including corre- 
lations, to obtain the final result 



^(£1" ^ K- 



(3.90 ±0.09 ±0.12)%. 



This result includes corrections (1.9%) due to final state 
radiation. The largest systematic uncertainties come 
from the background shape in extracting the inclusive 
D* yield and the modeling of the angle between the D* 
and the jet thrust axis. 

ARGUS used the same technique to count 
D°TT+ decays. To extract the D*+ D°tt+ yield AR- 
GUS plot the distributions of \cos9\ where 9 is the angle 
between the slow pion candidate and the thrust axis of 
the jet in the opposite hemisphere. Figure [T51 shows the 
I cos 9\ distribution in two ranges of the slow pion momen- 
tum. In the momentum range 0.2 to 0.3 GeV/c a clear ex- 




FIG. 12 The pi distribution for data from the ALEPH exper- 
iment. In (a) the transverse slow pion momentum squared in 
six equal momentum bins from 1.0 to 4.0 GeV/c. The 1.0 to 
1.5 GeV/c momentum bin is the uppermost and the 3.5 to 4.0 
GeV/c bin is the lowest. The slow pion from D*^ D^-k^ 
is clearly visible in the lower momentum range. In (b) the 
transverse momentum distribution s from d ifferent sources of 
D* mesons are shown. From iBarate et al\ ((T997i ) . 



cess of events near |cos6'| = 1 is seen from D*+ D^tt^ 
decays. In the range 0.4 to 0.5 GeV/c no excess is seen 
as this is above the momentum where we have slow pions 
from D*"*" decays. From a fit to the |cos6'| distribution 
ARGUS determines a yield of 51, 327±757 £»*+ D°n+ 
decays in the sample. The systematic uncertainty on this 
yield is estimated to be 5.9% by varying the signal shape 
parameterization. ARGUS reconstructs the D° in three 
channels and determines the following branching frac- 
tions 



B{D" K-TT+TT-TT^ 



= (3.41 ±0.12 ±0.28)% 
= (6.80 ±0.27 ±0.57)% 
= (5.03 ±0.39 ±0.49)% 



The CLEO (|Akerib et al\ . Il993[) study is very similar 
to the ARGUS analysis. CLEO only studied the final 
state D" ~* K^TT+. They tagged 165, 658±1, 149 D*+ 
D^TT^ decays and measured the branching fraction 

B{D^ K-TT+) = 3.95 ± 0.08 ± 0.17%. 

This includes a correction of about 1% for the effects 
of final state radiation. The largest contribution to the 
systematic uncertainty (±4.0%) comes from the track re- 
construction efficiency for the final Ktt system. 

These measurements are limited by systematic uncer- 
tainties on the determination of the number of D*^ — > 
D'^tt'^ decays in the data sample. The yield is extracted 
by extrapolating the background into the signal region 
based on shapes determined from Monte Carlo simula- 
tions. 



B. Tagging with B° D*+rv 

Tagging semileptonic B decays with the presence 
of a lepton plus a slow pion was first used by 
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TABLE IV Event yields and branching fractions for the ALEPH study of the K~tt'^ decay in bins of the slow pion 

momentum. The first column is the momentum range, the second and third columns show the D*^ D^-k^ yield determined 
from the slow pion transverse momentum and the — > K~tv^ yields, respectively. The last column shows the D° — > K~-k^ 
branching fraction. 



Momentum Range 






B{D'> ^ K--K+) (%) 


(GeV/c) 









1.0- 
1.5- 
2.0- 
2.5- 
3.0- 
3.5- 



-1.5 
-2.0 
-2.5 
-3.0 
-3.5 
-4.0 



79, 038.2 ± 2, 021.9 ± 12, 018.0 
56, 393.2 ±1,140.4 ±921.6 
35, 303.4 ±855.8 ± 842.2 
12, 287.8 ±674.7 ± 535.1 
3, 497.4 ±499.2 ±630.4 
192.4 ± 366.8 ±401.5 



2, 472.9 ± 55.5 ± 11.0 4.400 ± 0.150 ± 1.041 
1, 558.3 ± 41.4 ± 5.4 3.990 ± 0.133 ± 0.139 



913.8 ±30.9 ±2.8 
321.5 ± 18.2 ± 1.3 
115.7 ± 10.9 ±0.7 
9.8 ±3.3 ±0.4 



3.768 ±0.157 ±0.150 
3.758 ± 0.296 ± 0.206 
5.010 ± 0.857 ±1.228 
7.44 ± 14.2 ± 19.4 
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FIG. 13 The |cos0| distribution for data from the ARGUS 
experiment. In (a) the distribution is shown for the slow pion 
momentum in the range 0.2 t o 0.3 GeV/c and in (b) for the 
range 0.4 to 0.5 GeV/c. From I Albrecht et al\ (|l994bl ). 



ARGUS (lAlbrecht et all Il994al) a nd has since been 



used by CLEO (lArtuso aU Il998l) and most recently 



BABAR ()Aubert et all l2008c[) . The BABAR analysis 
uses the largest data sample, 210 fb^^ of e+e^ data col- 
lected at the T(4S'). 

In the first study that used this technique ARGUS used 
a sample of 246 pb~^ of e+e" data collected at the T(45') 
containing 209, 000 ±9, 500 BB pairs. They obtained the 
branching fractions 

B{D'^ ^ K-7:+) = (4.5 ±0.6 ±0.4)%, 
B{D° ^ K-TT+TT+n-) = (7.9 ±1.5 ±0.9)%. 

This measurement is clearly statistics limited, ARGUS 
reconstructed a sample of 2, 693± 183± 105 £>*+ D^n+ 
candidates. 



CLEO used a sample of 3.1 fb^^ of e+e^ data collected 
at the T(45') containing 3.3 x 10^ BB events. A sample 
of 1.6 fb~^ of data collected below the T{AS) resonance 
was used for continuum subtraction. CLEO reconstructs 
44, 504 ± 360 inclusive events and 1, 165 ± 45 exclusive 
D'^ — > K^TT'^ decays and determines a branching fraction 



6(Z?° 



) = (3.81 ±0.15 ±0.16)%. 



This branching fraction does not include radiative cor- 
rections. 

BABAR used 210 fb~^ of e+e" data coUected at the 
T(45') resonance, corresponding to 230 x 10^ BB pairs, 
and 22 fb"^ collected 40 MeV below the resonance. 
The offresonance sample is used to subtract non-BB 
backgrounds. In this analysis the semileptonic B de- 
cay, B° ^ D*+e-D followed by D*+ ^ D^n+ is used. 
BABAR use the lepton in the B decay and the slow pion 
from the D* to count 5° — > D*~^£~9 decays followed by 
D*+ D^TT+. BABAR used both electrons and muons 
in the momentum range 1.4 < |p|^ < 2.3 GeV/c. For 
the soft pion candidate the momentum is in the range 
60 < IpItt^ < 190 MeV/c. As the energy release in the 
Z)*+ D^-K^ decay is very small the reconstructed slow 
pion direction is used to approximate the direction of the 
Z)*"*". The momentum magnitude of the £)*+ is parame- 
terized as a linear function of the slow pion momentum. 
Using this estimate of the momentum, the missing 
mass squared of the neutrino is approximated as 



(-Ebc 



Ed' - Elf - {pD" + P«) 



where -Ebcam is half the center-of-mass energy and the 
momentum of the B is taken to be zero. The ener- 
gies and momenta in this expression are evaluated in 
the e+e^ center-of-mass frame. For signal candidates 
it is required that the charge of the slow pion and the 
lepton are opposite. For background studies BABAR 
considers same-charge candidates. BABAR extracts the 
number of D*~^£~D decays using the missing mass 

squared, Af^, against the D* and the lepton. Besides 
the _B" D*^£^D signal events there are a few ad- 
ditional sources of events that peaks near zero in the 
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FIG. 14 The distribution of the missing mass squared, for 
(a) right sign events and (b) wrong sign events. The wrong 
sign events sh ow that the simulation of the background shape 
is good. From lAubert et al.\ l|2008d '). 



missing mass squared. BABAR includes these events 
as signal candidates, they include 1) i? — > Z3*+(n7r)£"i/ 
O'D**") where n > 1; 2) B ^ D*+D, D -> i-X; 3) 
B" ^ D*+T-9, T- l-vivr ("cascade"); 4) 5° 
D*~^h~ ( "fake-lepton" ) , where h~ is a kaon or pion that 
has been misidentified as a lepton. The distribu- 
tions are shown in Fig. [Ml A clear signal is observed for 

> — 2.0 GeV^. However, there are substantial back- 
grounds from combinatorics in BB events and in contin- 
uum production that need to be subtracted. The con- 
tinuum background is modeled using offresonance data 
and the BB combinatorial background, as well as the 
signal components, are modeled using Monte Carlo sim- 
ulations. The signal yields are extracted from fits to the 

distributions in the range from —10.0 to 2.5 GeV^. 
The data are divided into ten different lepton momentum 
ranges to reduce sensitivity to the Monte Carlo simula- 
tion. In each lepton momentum bin the continuum yields 
are fixed by scaling the off-resonance sample to the lu- 
minosity of the on-rcsonance sample; while the number 
of events from primary signal, D**, and combinatorial 
BB are independently varied. The contributions from 
cascades and fake-leptons are fixed from the simulation. 
These two contributions account for about 3% of the to- 
tal inclusive signal. 

Table |V] summarizes the event yields for the inclusive 
D*^£~v reconstruction in the column 'Inclusive'. 
BABAR finds 7V'"'=i = 2, 170, 640 ± 3, 040 ± 18, 100 ^ 
P decays followed by D*^ Z?°7r+ in their data 
sample. 

The next step in this analysis is to reconstruct the 
_D" — > K'TT^ decay. All reconstructed charged tracks in 
the event are considered except for the tracks associated 
with the lepton and slow pion candidates. Pairs of tracks 



TABLE V Event yields for the inclusive B° -» D*+rv 
reconstruction and the exclusive analysis where the 
K~Tr'^ final state is reconstructed in the BABAR analysis 
to determine the branching fraction for K~-k^ decay. 

Errors are only statistical. 



Source 


Inclusive 


(xlO") 


Exclusive 


(xlO") 


Data 


4.4124: ± 


0.0021 


4.727 ± 


0.022 


Continuum 


0.46 ± 


0021) 


0.309 ± 


0.017 


Combinatorial BB 


1.7817 ± 


0.0007 


0.819 ± 


0.005 


Peaking 






0.163 ± 


0.008 


Cabibbo suppressed 






0.055 ± 


0.001 


Signal 


2.1706 ± 


0.0030 


3.381 ± 


0.029 



with opposite charge are combined, and the track with 
the opposite charge with respect to the slow pion candi- 
date is assigned the kaon mass. The kaon candidate is 
required to satisfy loose kaon identification criteria that 
retain more than 80% of real kaons while rejecting 95% 
of pions. The kaon plus pion invariant mass is required 
to satisfy 1.82 < Mk-k < 1-91 GeV. Each D° candidate 
is combined with the slow pion and the mass difference 
AM = M{K-7:+t:+) - M{K-n+) is computed. The 
signal is looked for in the range 142.4 < AM < 149.9 
MeV. 

Besides the signal events, the exclusive sample con- 
tains: continuum, combinatorial BB, uncorrelated peak- 
ing and Cabibbo suppressed decays. As for the in- 
clusive sample, the continuum background is subtracted 
using the off-resonance sample. The combinatorial BB 
background is determined from simulated BB events, 
normalized in the AM sideband 153.5 < AM < 162.5 
MeV. The background from uncorrelated peaking 
arises from events where the and lepton comes from 
different B mesons. This background peaks in AAf but 
not in M^. This background is estimated using the side- 
band in . The backgrounds from Cabibbo suppressed 
Z?" K~ and D'^ tt^tt^ decays are subtracted 
using simulated events. 

The mass difference, AM, is shown in Fig. [151 The 
yields for this 'exclusive' sample are given in Table |Vl 
After background subtraction BABAR finds A^"''^' = 
(3.381 ± 0.029) X 10"' events, where the uncertainty is 



only statistical. The branching fraction for I?" 
is calculated using 



K- 



6(i?° 



N' 



cxcl 



where exir = (36.96 ± 0.09)% from simulation and ^ — 
1.033 ± 0.002 is the selection bias for the partial recon- 
struction. The selection bias stems from the fact that 
the reconstruction efficiency for the slow pion is larger in 
events where the D'^ K~7t~^ than in generic D decays 
with more tracks. 

BABAR has considered many sources of systematic 



uncertainties that affects the measured D 







K' 
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where and Bi are the efficiencies and branching frac- 
tions for mode i and Nj^q is the number of produced 
DD pairs. Though the yields are determined separately 
for D and D decays it is assumed that the branching 
fractions are the same. Similarly, CLEO-c reconstructs 
double tags where both D mesons are reconstructed. The 
number of double tags found is given by 



N,, 



where i and j label the D and D mode used to reconstruct 
the event and Cij is the efficiency for reconstructing the 
final state. Combining the two equations above allow us 
to solve for A^nn as 



N 



DD 



FIG. 15 The AM distribution for the reconstructed D° - 
K~tt' ^ candidates in eve nts with a B° — > D*'^i~v tag 
From lAubert et al\ (|2008d ) . 



branching fraction. The most important uncertainties 
include: selection bias (±0.35%), nonpeaking combi- 
natorial background (±0.89%), peaking combinatorial 
background (±0.34%), tracking efficiency for kaon and 
pion (±1.00%), K- identification (±0.70%), D° invari- 
ant mass selection (±0.56%), and final state radiation in 



uncertainty is estimated to be ±1.80%. BABAR obtains 
the final result 

B{D° K-'K+) = (4.007 ± 0.037 ± 0.072)%. 



C. Absolute D hadronic branching fractions using double 
tags 



CLEO-c (iDobbs et a/.l . l2007HHe et aLl . l2005D has used 
a double tag technique, where by reconstructing one D 
in the event the presence of an additional D in the event 
is tagged. By determining how often the other D meson 
can be reconstructed in the event the branching frac- 
tion for the D decays can be calculated. This type of 
analy sis was first pione e red by the Mark I II col labora- 
tion (jAdler all Il988l : iBaltrusaitis et all Il98d ). The 
CLEO-c analysis described here uses the same basic idea. 

The CLEO-c analysis determines the number of single 
tags, separately for D and D decays. 



DD 



and 



This gives us the number of produced DD events. Note 
that many systematic uncertainties cancel in the ratio of 
efficiencies. This includes for example track finding ef- 
ficiencies and particle identification that are common to 
efficiencies in the denominator and numerator. However, 
systematic uncertainties from, for example, the determi- 
nation of the yields do not cancel as they are not cor- 
related. In this analysis CLEO-c determines all the sin- 
gle tag and double tag yields in data and the efficiencies 
from Monte Carlo simulations. The branching fractions 
and DD yields are extracted from a combined fit to all 
measured data yields and efficiencies. 

This analysis uses three D^ decay modes (i^"7r+, 



K- 



and K Tr~^ it tt"*") and six decay modes 



the Z?" K'Tr+ decay (±0.50%). The total systematic (K 7r+7r+, K 



.B,Nr 



and K^TT^). The 7r° candidates are reconstructed 
in the 77 final state, and the candidates are recon- 
structed in the tt+tt" final state. Particle identification 
criteria are applied on kaons and pions (excluding pions 
in candidates). A mode dependent selection criteria 
on A_E, the candidate energy minus the beam energy, is 
applied. The precise criteria are listed in Table IVII To 
extract the signal yields fits are performed to the M-qc 
distributions for the candidates that pass the selection 
criteria. The fit is described in Sect. IIII.A.3I The fit 
is performed separately for D and D candidates in each 
mode. These fits are shown in Fig. [TBI where the D and 
D decays have been combined. In Table IVIII the signal 
efficiencies and data yields are are shown for all single 
tag modes. Many backgrounds have been considered in 
this a nalysis and are discussed in detail in iDobbs et al\ 
(I2007D . These backgrounds are summarized in Table IVlTl 
The double tag yields are determined separately for 
the 45 = 3^ + double tag modes. The same criteria 
on A_E that was applied for the single tags are applied to 
the double tags. This ensures that the systematic uncer- 
tainty from the selection in single and double tag yields 
cancels in the ratio for the signal mode. To extract the 
number of double tag candidates a two-dimensional un- 
binned maximum likelihood fit is performed in the plane 
of Mbc{D) vs. Mbc{D). This is illustrated in Fig. [TT] 
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TABLE VI Requirements on AE for D candidates in the 
CLEO-c analysis for the absolute D*^ and Z)+ branching frac- 
tions. The limits are set at appro ximately 3 standard devia- 
tions of the resolution. From Dob bs et al\ |2007). 



Mode 






Requirement (GeV) 




K-TT+ 




\AE\ < 0.0294 




K~TT + 




-0.0583 < AE < 0.0350 




K-TT + 




\AE\ < 0.0200 


D+ 


K'TT^ 


-7r+ 


\AE\ < 0.0218 


D+ ^ 






-0.0518 < AE < 0.0401 


D+ 






\AE\ < 0.0265 


D+ ^ 






-0.0455 < AE < 0.0423 


D+ ^ 






\AE\ < 0.0265 


D+ ^ 


K+K 




\AE\ < 0.0218 



The signal peaks at Mbc{D) — Mbc{D) — Md- Beam 
energy smearing affects both Mbc{D) and Mbc(-D) in a 
correlated fashion to spread the signal along the Mbc [D) 
vs. Mbg{D) diagonal. In addition, the effects of initial 
state radiation will spread the signal along the same di- 
agonal to larger values of Mbg{D) and AfBc(^)- If all 
particles produced in the e+e^ interaction are used to 
form the D and D candidate, but the particles are ei- 
ther from continuum, or from a DD event but not as- 
signed to the right D candidate (mispartitioning) the 
reconstructed Mbc{D) and Mbc{D) will lie on the di- 
agonal. There are also events in which one of the two 
D candidates are misreconstructed. These events form 
horizontal and vertical bands in Mbc(-D) vs. Mbc{D)- 

The combined double tag data with the sum of the fits 
are shown in Fig. [TH] for the D^D^ and D^D~ modes. 
There are a total of 13, 591 ± 119 D^D^ double tags and 
8, 870 ± 96 D+D^ double tags. For most of the modes 
studied in this analysis the statistical uncertainty on the 
measured branching fraction is limited by the number of 
double tags. For the modes this statistical uncer- 
tainty is ±0.88% and for the _D+ modes this is ±1.1%. 

A detailed study of systematic uncertainties has been 
performed. Some of the studies performed by CLEO-c to 
determine the systematic uncertainties are described in 
Sect. IIII.AI The systematic uncertainties considered in 
this analysis are presented in Tables IVIIII and IIXI where 
the first table list systematics common to all modes and 
the second table list systematics that are mode specific. 
The signal shape systematic uncertainty for double tags 
are taken to be ±0.2%, while for the single tags a range of 
systematic uncertainties from ±0.3%, for K~tt^, 
to ±1.3%, for D+ K~7t~^tt^tt'^ , are assigned. These 
systematic uncertainties were assigned based on trying 
alternative signal shape parameterizations in the fit. For 
the neutral D decays there is an uncertainty due to 'dou- 
ble Cabibbo suppressed interference'. The source of this 
uncertainty comes from the interference between signal 
decays and decays where both the D° and the l)° de- 
cays via doubly Cabibbo suppressed decays. The relative 
size of this interference is A w 2i?ws cos 26 where i?ws is 
the ratio of the doubly Cabibbo suppressed rate to the 



3970407-010 




(GeV/c2) 

FIG. 16 Distributions of measured Mbc{D) and Mbc (5) 
values for single tag and candidates with D and D 
candidates combined in each plot. The points are data and 
the curves are fits to the data. In each plot, the dashed curve 
shows the background contributions and the solid curve shows 
the sum of the background and signal function. The number 
of eve nts is shown on a square-root scale. From fPobbs et al\ 
(|2007l ). 



Cabibbo favored rate and 6 is the relative strong phase 
between the doubly Cabibbo suppressed amplitude and 
the Cabibbo favored amplitude. CLEO-c assigns a sys- 
tematic uncertainty of ±0.8% for this effect. This covers 
the range of allowed values of A for i?ws = 0.004 and 
incorporates the uncertainties in S. 

For the charged track reconstruction CLEO-c assigns 
±0.3% uncertainty and for charged kaons an additional 
±0.6% added in quadrature. In addition CLEO-c as- 
signs a ±1.8% uncertainty on the Kg reconstruction in 
the 7r+7r~ final state and a ±2.0% uncertainty for the 
tt" reconstruction in the 77 final state. These system- 
atic uncertainties were discussed in Sect. IIII.A.41 Kaons 
and pions, except for pions in the reconstruction of 
Kg — > Tr~^iT^ candidates, are required to satisfy parti- 
cle identification criteria. Uncertainties of ±0.25% and 
±0.3% respectively for pions and kaons are assigned for 
the particle identification. 

For the uncertainty of the modeling of the AE selec- 
tion criteria an uncertainty of ±1.0% is assigned for the 
D+ i^gTr+TfO and D+ K+K-tt+ decays, for all 
other modes an uncertainty of ±0.5% is assigned. These 
uncertainties are taken to be correlated across all modes. 

The uncertainty in the background shape for the single 
tag fits is determined by using alternative descriptions of 
the background shape. CLEO-c uses sidebands in AE 
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TABLE VII Single tag efficiencies, yields from data, and peaking background expectations. The efficiencies include the 
branching fractions for tt" ^ 77 and K'g tv^tv~ decays. The entries in the column labeled "Background" are the number of 
events in the signal peak produ ced by non-signal e vents and the associated systematic uncertainty. The quoted yields include 
these background events. From [Pobbs et a/.l ((20071 ). 

Single Tag Mode Efficiency (%) Data Yield Background 





K~7V + 






64.18 ±0.19 


25, 760 


± 


165 


96 ± 27 




K+n- 






64.90 ±0.19 


26, 258 


± 


166 


96 ± 27 


D° 


K~TV + 


tt" 




33.46 ±0.12 


50, 276 


± 


258 


114 ± 10 


D° 


K+n- 


tt" 




33.78 ±0.12 


50, 537 


± 


259 


114 ± 10 


D" 


K-TI+ 


7r+ 




45.27 ±0.16 


39, 709 


± 


216 


889 ± 135 


W 


K+n- 


7r~ 


■K+ 


45.81 ± 0.16 


39, 606 


± 


216 


889 ± 135 


D+ 


K-n^ 






54.07 ±0.18 


40, 248 


± 


208 


< 1 


D- 


K+TT- 






54.18 ±0.18 


40, 734 


± 


209 


< 1 


D+ 


K--n^ 






26.23 ±0.18 


12, 844 


± 


153 


< 1 


D- 


K+n- 






26.58 ±0.18 


12, 756 


± 


153 


< 1 


D+ 








45.98 ±0.18 


5,789 


± 


82 


81 ± 22 


D- 








46.07 ±0.18 


5,868 


± 


82 


81 ± 22 


D+ 








23.06 ±0.19 


13, 275 


± 


157 


113 ± 53 


D- 








22.93 ±0.19 


13, 126 


± 


155 


113 ± 53 


D+ 


A-«7^^ 


-7r+ 


7r~ 


31.70 ±0.24 


8,275 


± 


134 


173 ± 83 


D- 


^ A-Ott- 




7r+ 


31.81 ± 0.24 


8,285 


± 


134 


173 ± 83 


D+ 


K+K 


~7r 


■f 


45.86 ± 0.36 


3,519 


± 


73 


< 1 


D- 


K-K+TT 




45.57 ±0.35 


3, 501 


± 


73 


< 1 
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FIG. 17 Scatter plot of Mbc(D) vs. Mbc{D) for D°D° dou- 
ble tag candidates. Signal candidates are concentrated at 
Mbc(D) = Mbc{D) = Md- The signal shape and dif- 
ferent ba^kgrom^ contributions are discussed in the text. 
From iDobbs et aH (,2007,) • 



to parameterize the combinatorial background shape de- 
scribed by the ARGUS function. Using the background 
shapes fixed from the AE sidebands to fit the signal re- 
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FIG. 18 Projections of the double tag candidate mass on the 
Mbc{D) axis for (a) the nine D^D^ double tag modes and 
(b) the 36 D'^ D~ double tag modes. The points show the 
data and the curves the projection of the fit results. The 
dashed lines shows the background contributions and the solid 
line th e signal shape plus the background. From lOobbs et all 
(|2007l ). 



gion in AE CLEO-c determines mode dependent system- 
atic uncertainties from ±0.4% to ±1.5%. 

In this analysis only one candidate per mode per event 
is selected. If there are multiple candidates the candi- 
date with the smallest AE is selected for single tags. 
As this criteria does not always pick the right candidate 
and the Monte Carlo simulation might not properly sim- 
ulate the rate of multiple candidates in an event there 
is a systematic uncertainty associated with the multiple 
candidate modeling. CLEO-c uses Monte Carlo simula- 
tions to determine the probability that the wrong candi- 
date is selected in events with multiple candidates. This 
is used, together with the difference between data and 
Monte Carlo of the rate of events with multiple candi- 
dates, to estimate the uncertainty associated with the 
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multiple candidate resolution. The uncertainty for many 
modes are negligible due to a very small rate of multi- 
ple candidates. The largest uncertainty, ±0.8%, is in the 
Z?" i^r^Tr+vr" mode due to fake 7r° candidates. 

Multibody final states suffer from an uncertainty in the 
simulation of the efficiency due to imperfect modeling of 
the resonant substructure. The uncertainties associated 
with the three- or four-body final states were estimated 
by comparing the kinematic distributions in these decays 
between data and Monte Carlo simulations. Many three- 
body final stats has been studied using E )alitz plot fit s 
and are well described in the Monte Carlo ( Langd . [20011 ). 
The Dalitz plot analyses are described in Sect. IIXI 

Last, final state radiation, as discussed in Sect. IIII.DI 
was considered. CLEO-c compared the signal efficiencies 
with and without FSR included in the Monte Carlo sim- 
ulation. A systematic uncertainty of ±30% of the change 
due to not including final state radiation was assigned. 
This gives the largest uncertainty of about 0.9% in the 
D° K~TT+ mode. 

The signal yields for single and double tags and the 
efficiencies determine d from Mon te Carlo simulations are 
combined in a fit ( Sunl . [20061 ) . This fit includes both 
statistical and systematic uncertainties. The fit extracts 
the branching fractions for the nine D decay modes stud- 
ied in this analysis and the produced number of D^D^ 
and D^D~ pairs. The result of this fit is shown in Ta- 
ble [Xj The correlation coefficients for the measurement 
is shown in Table KIl The of the fit is 39.2 for 52 
degrees of freedom, corresponding to a confidence level 
of 98%. The includes systematic uncertainties. 

The CLEO-c analysis obtains the main branching frac- 
tion results 

B{D° ^ K-n+) = (3.891 ±0.035 ± 0.059 ± 0.035)%, 
B{D+ ^ K-7r+n+) = (9.15 ± 0.10 ± 0.16 ± 0.07)%, 

where the errors are statistical, systematic, and from final 
state radiation respectively. In addition the DD yields 
determined from this analysis are used to normalize many 
other CLEO-c measurements. The cross-sections for DD 
production are discussed in Sect. IIII.AI 



D. Summary of ^ K~'k+ 

The absolute branching fraction for D^ K~Tr~^ has 
been measured by many different experiments, using dif- 
ferent techniques as discussed in this Section. The differ- 
ent measurements are summarized in Table IXlIl The two 
most recent, and most precise, measurements are from 
CLEO-c and BABAR. They use very different techniques 
but finds branching fractions that are in good agreement. 
We adopt the PDG average 



) = 3.89 ±0.05. 



These measurements are now limited by systematic un- 
certainties. There are many sources of systematic uncer- 
tainties that contribute. Some of these can be improved 



with additional data. Both CLEO-c and BABAR can 
increase the data samples used in their analyses. 



E. Modes with K? or K'i in the final states 



It has commonly been assumed that T{D 
r{D ^ Kj X). However, as pointed out by Bigi and Ya- 
mamoto ( Bigi and Yamamotol . 1 199*51 ) this is not generally 
true as for many D decays there are contributions from 
Cabibbo favored and Cabibbo suppressed decays that in- 
terfere and produce different rates to final states with 
versus K^. As an example consider D^ Kg ^7r°. Con- 
tributions to these final states involve the Cabibbo fa- 
vored decay A'°7r° as well as the doubly Cabibbo 
suppressed decay D° ^_i4r°7r°. However, we don't ob- 
serve the K° and the K° but rather the and the 
Kl- As the amplitudes for D^ K°tt° and Z)° ^ K^n^ 
interfere constructively to form the Kg final state, and 
destructively to form a K^, we see a rate asymmetry be- 
tween the K^ and Kg final states. Using SU(3), and 
in particular the U-spin subgroup, one can predict the 
asymmetry in D'^ — > K 



S,L' 



r{D° 



T{D° -> KItt^) 



T{DO i^OTr") + r(DO ^ KIttO) 
2tan2 Oc = 0.109 ±0.001. 



For the corresponding charged D mode, _D+ — > Kg ^tt"*" a 
similar prediction based on SU(3) is not possible. Rather 
one has to rely on calculations based on factorization. 
Some of these predictions are discussed below. 

Experimentally these channels are challenging as they 
involve final states with a K^. CLEO-c has studied these 



modes ( He et al ],[200i). They infer the presence of a K^ 
using a missing mass technique after vetoing events with 
a Kg decaying to either a tt+tt^ or tt'^tt'^ pair if there are 
tracks or tt'^ candidates reconstructed in the event. 

In addition to the challenge with the K^l final state, 
these decays are CP eigenstates and at the -0(3770) 
where CLEO-c recorded the data for their analysis we 
need to disentangle the effects from quantum coherence 
with the rate asymmetry we are interested in here. The 
effect of the coherently produced D^D^ pairs at the 
-0(3770) was discussed in Section lTlI.A.ll 

CLEO-c has studied both D+ K%j^-k+ and D° 



K^ 



TT". First the D'^ analysis is discussed as it 
does not involve the complication of quantum coher- 
ence. The branchi n g frac tion for — > K^ir'^ is taken 
from iDobbs et al\ (|2007t) . In this analysis only the 
branching fraction for D"*" A'^tt"*" is directly measured. 
CLEO-c uses a tag technique, in which one charged D 
is fully reconstructed. Six different charge d D tags are 
used, these modes are the same as in Ref. ( Dobbs et al\ . 
l2007l ) described in Sect. IVTCl The tag D' is combined 
with a TT"*" and events consistent with a K^ are vetoed. 
An event is vetoed if an additional charged track or neu- 
tral pion, reconstructed in the tt^ — > 77 channel, was 
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TABLE VIII Systematic uncertainties and the quantities to which they are apphed in the branching fraction fit. Uncertainties 
not correlated between decay modes are given in the first section, and correlated uncertainties in the second. The symbols y and 
e denote yields and efficiencies, respectively. Yield uncertainties are additive and efficiency uncertainties are multiplicative. See 
the text for the distinction between e(Charged) and t{K'^). The detector simulation uncertainties are determined per charged 
track or per neutral pion or kaon. Uncertainties for other efficiencies are determined per D. In additi on, to the sy s temat ic 
uncertainties listed here, five more mode-dependent systematic uncertainties are listed in Table llXl From lDobbs et al\ \2QQt ). 



Source 


Uncertainty (%) 




Quantity or Decay Mode 


DT Signal Shape 

Double DCSD Interference 


0.2 
0.8 




y(All DT Modes) 
y(Neutral DT) 


Detector Simulation 


0.3 




e(Charged) Tracking 




0.6 




e{K^) Tracking 




1.8 




<K%) 




2.0 








0.25 




e(7r±) PID 




0.3 




e{K^) PID 


Lepton Veto 


0.1 




e{D° K--K+) ST 


Trigger Simulation 


0.2 




e{D° -> K-Tv+7v°) 




0.1 




e{D+ A'§7r+) 


AS] Requirement 


1.0 




K'^ tt+tt") and e{D+ K+K-n+) 




0.5 




e(All Other Modes) 



TABLE IX Mode-dependent systematic uncertainties. The systematic uncertainties for the signal shapes are correlated among 
all ST modes. The systematic uncertainties for FSR are correlated among all ST and DT modes. Other uncertainties are uncor- 
related. The background and signal shape uncertainties are uncertainties on the yields, the other uncertainties in the table are 
uncertainties on the efficiency. Yield uncertainties are additive and efficiency uncertainties are multiplicative. From lPobbs et al\ 
12007). 



Mode 




Background 


ST Signal 


FSR (%) 


Resonant 


Multiple 






Shape (%) 


Shape (%) 




Substructure (%) 


Candidates (%) 


D"^ 


K-TT + 


0.4 


0.3 


0.9 




0.0 






1.0 


0.5 


0.3 


0.3 


0.8 


D"^ 


TX^ TT^ ■K~ 


0.4 


0.7 


0.8 


1.2 


0.0 


D+ - 


K--K + -K + 


0.4 


0.3 


0.7 


0.6 


0.0 


D+ - 


K--K + -K + -K° 


1.5 


1.3 


0.3 


0.5 


0.5 


D+ - 


K%^ + 


0.4 


0.4 


0.5 




0.2 


D+ - 




1.0 


0.5 


0.1 


1.2 


0.0 


D+ ^ 


A'g TT + TV + TT- 


1.0 


0.6 


0.6 


0.5 


0.0 


D+ - 


K+K-TV+ 


1.0 


0.6 


0.3 


1.3 


0.2 



found. This veto removes about 90% of the Kg back- 
ground as well as many other backgrounds while retain- 
ing 98% efficiency for signal events. 

Figure[Tn]shows the invariant mass distribution against 
the tag D and charged pion. The signal peaks at a miss- 
ing mass square of about 0.25 GeV^ corresponding to the 
K^. From the fit to the data CLEO-c extracts a signal 
of 2, 023 ± 54 events. With 165 x 10^ charged D tags and 
an efficiency of 81.6% for finding the pion the branching 
fraction is calculated to be 



signal peak as seen in Fig. [T9l 



KItt+) = (1.460 ± 0.040 ± 0.035 ± 0.0005)%, 



where the errors are statistical, systematic, and from the 
branching fraction for £)+ Kgir^. The largest con- 
tributions to the systematic uncertainty come from the 
extra track and 7r° veto (±1.1%) and the signal peak 
width (±1.6%). The sensitivity to the peak width comes 
from the rjir'^ events just on the high side of the 



if^TT"'" bra nching fracti o n with 



K^Tr+ measured in iDobbs et all (|2007[ ) 



Combining the 

the D+ 

CLEO-c obtains the asymmetry 



R{D+) = 0.022 ± 0.016 ± 0.018. 

There is no evidence for a significant asymmetry in the 
£)+ — > Kg ^7r+ mode. Predictions for the asymmetry 
in charged D decays is more involved than for neut ral D 
decays. D.-N. Gao, based on factorization, predicts ( Gaol . 
l2007i) this asymmetry to be in the range 0.035 to 0.044, 
which is consistent with the observed asymmetry. 

For the D'^ Kg analysis the effects of the quan- 
tum coherence has to be accounted for. In addition, ex- 
perimentally this mode is more challenging as the reso- 
lution for a 7r° is worse than for a charged pion. CLEO-c 
first measures the branching fraction for — s- K^t:^ 
without using a tag. Next the 'branching frac- 
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TABLE X Fitted branching fractions and DD pair yields. For Nj^q-qo and Njj+j^- , uncertainties are statistical and systematic, 
respectively. For branching fractions and ratios, the systematic uncertainties are divided into the contribution from FSR (third 
uncertainty) and all others combined (second uncertainty). The column of fractional systematic errors combines all systematic 
errors, including FSR. The last column, Apsr, is the relative shift in the fit results when FSR is not included in the Monte 
Carlo simulations used to determine efficiencies. From iDobbs et all (|2007l ). 



Parameter 


Fitted Value 


Fractional Error 


ApsR 






Stat.(%) 


Syst.(%) 


(%) 


B{D° -* K-TT+) 


(1.031 ±0.008 ±0.013) X 10" 


0.8 


1.3 


±0.1 


(3.891 ± 0.035 ± 0.059 ± 0.035)% 


0.9 


1.8 


-3.0 


B{D° K~Tv+n'-^) 


(14.57 ± 0.12 ± 0.38 ± 0.05)% 


0.8 


2.7 


-1.1 


B{D° K-Tv+TT+n-) 


(8.30 ± 0.07 ± 0.19 ± 0.07)% 


0.9 


2.4 


-2.4 


Nd+d- 


(0.819 ± 0.008 ± 0.010) X 10® 


1.0 


1.2 


±0.1 


B{D+ K-n+TT+) 


(9.15 ±0.10 ±0.16 ±0.07)% 


1.1 


1.9 


-2.3 


B{D+ A'-tt+tt+tt") 


(5.98 ± 0.08 ± 0.16 ± 0.02)% 


1.3 


2.8 


-1.0 


B{D+ Ks 7r+) 


(1.526 ± 0.022 ± 0.037 ± 0.009)% 


1.4 


2.5 


-1.8 


B{D+ K'^tv+tt") 


(6.99 ±0.09 ±0.25 ±0.01)% 


1.3 


3.5 


-0.4 


B{D+ ~* K'^Tv+TT+n-) 


(3.122 ± 0.046 ± 0.094 ± 0.019)% 


1.5 


3.0 


-1.9 


B{D+ K+K-n+) 


(0.935 ± 0.017 ± 0.024 ± 0.003)% 


1.8 


2.6 


-1.2 


BiD" K-tv+-k'>)/B{K-ti+) 


3.744 ± 0.022 ± 0.093 ± 0.021 


0.6 


2.6 


±1.9 


B{D° K~Tv+n+n-)/B{K-n+) 


2.133 ± 0.013 ± 0.037 ± 0.002 


0.6 


1.7 


±0.5 


B{D+ -> K--K+-K+n°)/B{K-Tr+Tv+) 


0.654 ± 0.006 ± 0.018 ± 0.003 


0.9 


2.7 


±1.4 


B{D+ K°sn+)/B{K'Tv+7v+) 


0.1668 ± 0.0018 ± 0.0038 ± 0.0003 


1.1 


2.3 


±0.5 


B{D+ K°sTV + TT°)/B{K-TV + TT+) 


0.764 ± 0.007 ± 0.027 ± 0.005 


0.9 


3.5 


±2.0 


bId+ K%-K + 1T+TT-)/B(K- 11+11+) 


0.3414 ± 0.0039 ± 0.0093 ± 0.0004 


1.1 


2.7 


±0.4 


B{D+ K+K~iT+)/B{K-iv+iT+) 


0.1022 ± 0.0015 ± 0.0022 ± 0.0004 


1.5 


2.2 


±1.1 



TABLE XI The correlation m atrfac, including syst ematic uncertainties, for the branching fractions and numbers of DD events 
determined from the fit. From iDobbs et al\ l|2007l ). 





TV^oTjo Kiv Kiriv" Kiviriv 


No+a- Km Kivinr° K%ir K%iiit° Tfg tttttt KKit 


B{K-^iv+) 

B{K-iv+iv°) 

B{K-iv+iv-n+) 


1 -0.65 -0.34 -0.41 
1 0.44 0.70 
1 0.38 
1 


0.39 -0.19 0.01 -0.14 -0.09 -0.08 -0.09 
-0.22 0.52 0.23 0.28 0.15 0.30 0.35 
-0.11 0.28 0.66 0.14 0.51 0.17 0.21 
-0.09 0.51 0.29 0.28 0.17 0.37 0.34 


Nn+D- 

B{K-iv+iv+) 

B{K-iv+iv+ir°) 

B{K%ii+) 

B{K°siT+iT°) 

B{K'^;iT+iT+ir-) 

B{K+K-iT+) 




1 -0.61 -0.24 -0.48 -0.30 -0.33 -0.38 
1 0.43 0.52 0.32 0.51 0.55 
1 0.27 0.56 0.29 0.32 
1 0.55 0.72 0.31 
1 0.50 0.20 
1 0.30 
1 



tion' for 13° K^t:^ is measured in a tagged analy- 
sis where the is reconstructed in three modes. Due 
to the coherence the 'branching fraction' measured in 
the tagged analysis is B{D^ K^gT:^){\ — Cj), where 
— (?'/-z/±y)/(l±^«)s,/), as described in Sect. lIII.A.ll 
For the three tag modes C/ can now be calculated. Fi- 
nally, the 'branching fraction' for 13° J^l'^'^ is mea- 
sured using the same three tag modes, each of the tag 
modes give us B{D^ if£7r°)(l ± C/), and using the 
measured values of C/ from above the branching frac- 
tion B{D^ i^^7r°) can be determined. 

The Kg is reconstructed in the Kg tt+tt" final state. 
There is a background from D'^ tt+tt^tt^. This back- 
ground is subtracted using the Kg mass sideband. The 
signal yield in this analysis is extracted using a cut-and- 



count technique. CLEO-c looks in a 3 standard devia- 
tion window around the nominal values for the beam- 
constrained mass and AE. A sideband in AE is used to 
subtract the combinatorial backgrounds. The num ber of 
D^D^ pairs in the data sample is taken from lDobbs et all 
(I2007D . CLEO-c obtains the branching fraction 

B{D° -> K^stt") ^ (1.240 ± 0.017 ± 0.031 ± 0.047)% 

where the last error is due to the tt*^ reconstruction ef- 
ficiency. In the asymmetry R^D'^) this uncertainty will 
cancel. 

Next the 'branching fraction' for B{D° K'^tt°) is 
measured with a tag. The three tags modes used are 
5" ^ K+TT-, 5° ^ K+Tr-n°, and 5° ^ K+n-n+n-. 
The results for the tagged analysis is summarized in Ta- 
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TABLE XII Summary of measurements of the 
are used in the average by the PDG. 



K TT^ branching fraction measurements. Only the top six measurements 



Experiment 


Ref. 




B{D"^K-Tv^) (%) 


CLEO-c 


Dobbs et al. 


(2007) 


3.891 ± 0.035 ± 0.059 ± 0.035 


BABAR 


Aubert et al. 


f2008c) 


4.007 ± 0.037 ± 0.072 


CLEO IP 


Artuso et al. 


(1998) 


3.82 ±0.07 ±0.12 


ALEPH 


Barate et al. 


(1997) 


3.90 ±0.09 ±0.12 


ARGUS 


Albrecht et al. 


(1994a) 


3.41 ±0.12 ±0.28 


ALEPH 


Decamp et al. (1991) 


3.62 ±0.34 ±0.44 



CLEO-c 
CLEO II 
CLEO II 
ARGUS 
CLEO II 
HRS 

MARK III 
MARK II 
LGW 



He et al. (2005 ) 
Artuso et al. 
Coan et al. 
Albrecht et al 
Akerib et al. 
Abachj^et al. 
A dler et al. 
Schi ndler et al. (1981) 
Peruzzi et al. ( 1977) 




3.91 ±0.08 ±0.09 
3.81 ±0.15 ±0.16 
3.69 ±0.11 ±0.16 

4.5 ±0.6 ±0.4 
3.95 ±0.08 ±0.17 
4.5 ±0.8 ±0.5 
4.2 ±0.4 ±0.4 
4.1 ±0.6 
4.3 ± 1.0 



Average 



3.89 ±0.05 



"This is an avera g e of the results in lAkerib et al\ lll993f) : 
lAubert et oZI ll2008d) : [Coan et al\ l ll998t) . 



TABLE XIH The efficiency is for the reconstruction of the 
Kgir^' after the D tag has been found, the tag yield is the 
number of D tags reconstructed, the signal yield is the number 
of D° ^ K°sTT° candidates are reconstructed against the tag 
D, and the tag bias is a correction due to the fact that it is 
easier to reconstruct the tag in ev e nts w ith the signal than in 
generic D decays. From iHe et"^ (|2008l '). 



TABLE XIV The efficiency is for the reconstruction of the 
A'^tt", including the Kg veto, after the D tag has been found, 
the tag yield is the number of D tags reconstructed, the sig- 
nal yield is the number of — > A'stt" candidates are recon- 
structed against the tag D, and the tag bias is a correction due 
to the fact that it is easier to reconstruct th e tag in e vents with 
the signal than in generic D decays. From iHe et al. {200& ). 



Tag mode 



K 



TT TT Tag mode 



Efficiency (%) 31.74 31.29 29.97 

Tag yield 47,440 63,913 71,040 

Signal yield 155 203 256 

Tag bias correction (%) 1.000 1.014 1.033 

g(£>° ^ K%-K°)(l - Cf) 1.03 ± 0.09 1.00 ± 0.09 1.16 ± 0.08 



K 



Efficiency (%) 55.21 52.72 49.88 

Tag yield 47,440 64,280 71,040 

Signal yield 334.8 414.5 466.5 

Tag bias correction (%) 1.000 1.037 1.057 

B{D° KItv°){1 + Cf) 1.28 ± 0.08 1.03 ± 0.06 1.12 ± 0.06 



ble IXIIII Similarly the tagged 'branching fraction' for 
— > K^TT*^ was studied using a missing mass technique 
where the event was fully reconstructed except for the 
K^. The results are summarized in Table IXlVl 

Combining these measurements CLEO-c finds an av- 
erage asymmetry for the neutral D decays 

R{D°) = 0.108 ± 0.025 ± 0.024, 

which is in good agreement with the prediction. 

F. Final states with three kaons 

Final states with three kaons are not generally Cabibbo 
suppressed, but the smaller branching fractions for these 
decays are due to the small phase space available in these 



decays. The decay D"*" K^K^K^ is Cabibbo sup- 
pressed and is included in Sect. lVII.DI The limited phase 
space availabl e has been taken adva ntage of to measure 
the D° mass (jCawlfield et aLl . l2007t ). 



G. Summary of Cabibbo favored D" and decays 



In Table I XVII a summary 
and decays are given. 



of the Cabibbo favored 
Assuming that V{D ^ 



K^X) = T{D K^X) for modes where the final states 
with a has not been explicitly measured the Cabibbo 
favored branching fractions adds up to (50.8 ± 1.4)% for 
meson decays and (38.3 ± 1.1)% for D+ decays. The 
mode if "tt+tt^tt" is not included here. An early 

measurement by MARK III ()Adler et all . Il988f) reported 



28 



TABLE XV Summary of final states with three kaons. If there are more than one measurement we quote here the PDG 
average. 



Mode 




Rof. 






K%K%K% 
K+K-K--K+ 


BABAR CAubert et al. 2005a~) 
PDG Ave. CAmsler et al. 2008 ~) 
PDG Ave. (Amsler et al. 2008) 
FOCUS (Link et al. 2005a) 


4.72 ±0.03 ±0.15 ±0.27 
0.96 ±0.12 ±0.05 
0.221 ± 0.033 ± 0.009 
0.63 ±0.11 ±0.06 ±0.04 



4110907-002 




TABLE XVI Summary of branching fractions for Cabibbo 
favored D° and decays. 



0.15 0.2 0.25 0.3 0.35 0.4 
Missing mass squared (GeV^) 

FIG. 19 Missing mass squared distribution for all six tag 
modes for Xn^ . Events w ith extra t racks or vr" can- 

didates have been removed. From iHe et al\ (|2008l ). 



a large branching fraction of 15 ± 5%. The PDG is not 
using this result anymore in their summary and there has 
not been any newer measurements. However, CLEO- 
c has used this mode for tagg ing £)° decays in their 
studies of semileptonic decays ( Ge et al\ . l2009l ). They 
provide enough information that the branching fraction 
B{D° K-tt+tt^tt") = (7.90 ± 0.14)% can be calcu- 
lated. The error quoted only includes the statistical error 
and the uncertainty from the K^tt^ normalization 

mode. In particular experimental systematic uncertain- 
ties are not included and hence this is not included in the 
summary. But it does show that there is a substantial 
branching fraction to the K~tt^tt^tt^ final state. 



VI. CABIBBO FAVORED Ds 
BRANCHING FRACTIONS 



DECAYS AND REFERENCE 



The determination of the absolute branching fraction 
scale for Dg deca ys has been a cha llenge since the dis- 
covery of the Ds ( Chen et~al\ . Il983l ). Until recently the 
focus has been on the final state — > 07r+ , followed by 
(j) K~ . This final state is easy to reconstruct with 
small backgrounds; the (/> is a narrow resonance and the 
final state consists of all charged particles. However, this 



Mode 


Ref. Branching Fraction 


DO _ K-TV+ 


(3.89 ± 0.05)% 


DO ^ A-§7rO 


(1.22 ±0.06)% 


DO -> 


(1.00 ± 0.07)% 


DO ^ A-gTT+TT- 


(2.99 ±0.17)% 


DO A-tt+ttO 


(13.9 ±0.5)% 


DO ^ K-TV + TV + TT- 


(8.10 ±0.20)% 


DO ^ 7i-07r+7r-7rO 


(5.4 ± 0.6)% 


DO ^ A-TT+TT+Tr^TrO 


(4.2 ± 0.4)% 


DO A-gr/vrO 


(5.6 ± 1.2) X 10"^ 


DO ^ K'^TT+TT + TV-TT- 


(2.84 ±0.31) X 10"^ 


DO ^ K-TV + TV + TT + TV'TT- 


(2.2 ±0.6) X 10"" 


D+ ^ A-07r+ 


(1.45 ± 0.04)% 


D+ KIt:+ 


(1.46 ± 0.05)% 


D+ K-TV + TV + 


(9.22 ±0.21)% 


D+ K°n+n° 


(6.8 ±0.5)% 


D+ ^ K-TV + TV+TT° 


(6.00 ± 0.20)% 


D+ ^ K'^s^+TT + TV' 


(3.02 ±0.12)% 




(5.6 ±0.5) X 10-3 



final state is not as 'clean' as one would wish. There are 
non-(j) contributions, such as the /o(980), to the K~^K~ 
mass near the mass that pollutes the Df 07r+ signal. 
Of course, these decays are still real Df K'^K~tt~ de- 
cays. This is discussed further in Sect. lIX.A.lll on Dalitz 
plot analysis of K^tt'^ . 

As measurements have gotten more precise the defi- 
nition of what is measured has had to be made more 
pr ecise. One of the most r ecent measurements by CLEO- 
c (Alexand er et all . |2008| ) does not quote a Df (j)TT~^ 
branching fraction, but rather a partial branching frac- 
tion in a invariant mass region near the 0. The 
first attempts at establishing the branching fraction scale 
was based on model dependent assumption about equal 
partial widths for semileptonic decays of the D+ and . 

This Section will discuss the different approaches used 
to measure the Dg absolute branching fractions. The 
early measurements are described very briefly and the 
more recent, and precise, measurements are described in 
more detail. 
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A. Model dependent approaches 

These measurements are typically no l onger used in 
avera ges, e.g. by the particle data gro up (lAmsler et al. 
20081). Several experime nts, CL EO (lAlexander et al 



ZUU8II. several experime nts, '^L tjD llAiexander et au . 
IQQOal : iButler et a/.[ll994|). E6 87 (iFrabetti et all Il993^ 



ARGUS (|Albrecht aZ.Ul99lh . and E691 (Anjos et all 
fl990i) measured the ratio 

B{Dt ^ (/)7r+) ■ 
Using theoretical predictions for the ratio 

r(£)+ ^ c^i+ui) 

T{D+ k*°l+vt) 

and the measured £)+ and £)+ lifetimes these experi- 
ments determined the branching fraction for — > ^tt"*" . 
Comparing these results require some care as slightly 
different assumptions were made about the ratio of the 
semileptonic rates. Also, combining these measurements 
require care as there are strong systematic correlations 
between the measurements due to the common, or at 
least similar, assumptions about partial rates for the 

semileptonic decays. 

The NA14 experiment (jAlvarez et all Il990l ) used the 
Lund model to estimate the ratio of D+ to D+ pro- 
duction cross-sections, which allowed them to determine 
the D+ (j)Tr~^ branching f ractions. The CLEO col- 
laboration ( Chen et~al\ . Il989l ) used estimates of the Df 
production rate to determine the branching fraction for 
D+ -> (j)T:+. 

All of these measurements use model dependent as- 
sumptions and have associated systematic uncertainties 
that are hard to quantify. With larger data samples 
model independent measurements became possible. 



B. The branching ratio for Ds 



Att from DID* 



The first statistically significant, see Sect. I VI. Dl model 
independent measurement of the ab solute Dt br a nchin g 
fraction was performed by CLEO ( Artuso et all [l996l ). 
They used 2.5 fb^i of e+e" data collected at the T(45) 
resonance, corresponding to 2.7 x 10® BB pairs, to study 
B — > D'*D* decays. The same te chnique has been used 
by BABAR (jAubert et all l2005cl) . They have analyzed 
a sample with (123 ± 1) x 10® BB pairs. 

In these analyses the decay B D*D* is recon- 
structed in two different ways. First, the D* is fully 
reconstructed using — > £'^7 followed by Df — > (pn^ 
and the D* is partially reconstructed using the slow pion 
from the D* decay. In the second method the D* is 
fully reconstructed and the Z?*^ Df'y is only iden- 
tified through the presence of the 7. From this study 
BABAR quotes B{D+ ^ 07r+) = (4.81 ± 0.52 ± 0.38)% 
and CLEO B{D+ -> (j)T:+) = (3.59 ± 0.77 ± 0.48)%. 
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FIG. 20 The recoil mass against a D or D* 

From lAubert et al\ (|2006d '). 



More recently, BABAR (jAubert et a/.l . l2006cD has pre- 
sented results based on 210 fb~-^ of data where they use 
a tag technique in which one B meson is fully recon- 
structed. In events with one fully reconstructed B me- 
son candidate BABAR reconstructs one additional D*^*^ 
or -Dsj-j) meson. Then they look at the recoil mass 
against this reconstructed candidate. The recoil masses 
are shown in Figs. [201 andj^H 

From these data BABAR extracts B{Dsj{2460y 
D*-n") = (56±13±9)%andi3(Lisj(2460)- -> D*-j) = 
(16 ± 4 ± 3)% in addition to B{D- (t>TT+) = (4.52 ± 
0.48 ±0.68)%. BABAR combines this measurement with 
their previous measurement discussed above to obtain 
^(£17 — - 



•TT+) = (4.62 ±0.36 ±0.50)%. 



C. Study of Dt K tt"*" in continuum production 



Belle (|Abe et all [2OO7I) has used 552.3 pb"! of e^ 



data to study the process e+e Dl~^D^^ followed 

by D*^K- and D*+ D+j. The very large 

data sample allow them to study this exclusive final 
state in continuum production of mesons. The final 
state is reconstructed in two ways; either by partially 
reconstructing the Dgi or the D*. Belle obtains the 
preliminary branching fraction B{Df — > K'^ K~Tr~^) = 
(4.0 ± 0.4 ± 0.4)% which is of comparable statistical pre- 
cision to the other methods discussed above. 
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FIG. 21 The recoil mass against a Ds 

From lAubert et al\ (|2006j ). 



or Dt 



D. Absolute branching fractions for hadronic Ds decays 
using double tags 



CLEO-c (jAlexander et all l2008f) has determined the 
absolute hadronic branching fractions for Ds meson de- 
cays using a double tag technique similar to what was 
done for the D hadronic branching fractions. The same 
technique wa s used by MARK III (jAdler et all Il990bl ) 
and BES ( Bai et all [1995). These initial studies were 
limited by statistics; MARK III observed no events and 
placed an upper hmit while BES observed two events and 
reported a branching fraction of 



-) ^ (3.9 



+5.1+1.8\ 
-1.9-1. 17 



recorded at 



4.03 



The BES analysis used 22.3 pb^ 
GeV. 

The CLEO-c analysis used a sample of 298 pb^^ of 
e"^e~ collision data recorded at a center-of-mass energy 
of 4170 MeV. At this energy Dg mesons are produced, 
predominantly, as DfD'*^ or DjD'*^ pairs. The eight 
hadronic decays considered in this analysis by CLEO- 
c are Df K°^K+ . Df Xg/s:-7r+7r+, D+ -> 



K+K-n+, D+ K+TT-7r+, D+ K+K-tt+tt^', 

Df 7r+7r-7r+, Df TjTr+ , and Df ?7'7r+. The 
analysis proceeds similar to the D hadronic branching 
fraction analysis described in Sect IV. CI Yields and effi- 
ciencies for single tags (separately for Df and D~) and 
double tags are extracted. The 7r° or 7 from the D* de- 
cay is not reconstructed in this analysis. The yields, in 
terms of the efficiencies, branching fractions, and data 



sample size are given by 



y^] 



where i indicates a Df and j indicate a D^ . In this anal- 
ysis a total of 16 single tags and 64 double tags are used . 
The event selection is detailed in lAlexander et aZI ( 2008f ). 
A Ds candidate is referred to as "indirect" if it comes 
from the decay of the D* in the e+e^ DgD* interac- 
tion. Otherwise the D^ is said to be "direct" . The Dg 
candidates are identified based on their momenta and in- 
variant mass. The direct Ds has a fixed momentum in the 
e'^e~ restframe, whereas the indirect Ds has momenta in 
a range due to the extra boost from the I?* Z3s(7, 7r°) 
decay. The recoil mass Mroc is defined by 



Ml 



- (Po - Pd, 



where (£^o:Po) is the e'^e~ center-of-mass four-vector, 
Po^ is the measured Ds momentum and Mo^ is the nom- 
inal Ds mass. For direct Ds candidates Mj-oc peaks at the 
Z?* mass of 2.112 GeV, while for indirect Ds candidates 
Mrec is spread about evenly over ±60 MeV around this 
peak. CLEO-c requires that Ds candidates in a double 
tag, and for most single tags, satisfies M^^c > 2.051 GeV. 
For the three single tag modes, K~ K'^tt'^tt'^ , tt+tt^tt^, 
and K'^TT'^TT~, with more substantial backgrounds it is 
required that M^^c is greater than (2.099, 2.101, 2.099) 
GeV, respectively. Note that this cut eliminates events 
from e+e" DfD~ as these events peak at Mj-^c = 
Mo^- A number of vetoes are applied to reject fake can- 
didates, primarily from D* D* events. 

The single tag signal yields are extracted from the Ds 
invariant mass distributions. The single tag event yields 
in data are shown in Fig. [221 At most one single tag 
candidate per mode and charge are accepted per event. 
If more than one candidate pass the selection criteria the 
candidate with the value of Mroc closest to M^. is se- 
lected. The data is fit to a signal shape and a background 
shape. The signal shape is determined from Monte Carlo 
simulations, but the Ds mass is allowed to float in the 
fit. The background is modeled with a linear function 



in all modes except K^K 



and 



where a 



quadratic form is used instead. The background shape is 
constrained to be the same in the two charge conjugate 
modes. 

The double tag yields are extracted by a cut-and-count 
procedure in the plot of the invariant mass of the Df vs. 
D^ . All double tag candidates are shown in Fig. (221 
At most one double tag candidate is allowed per event. 
If there are more than one candidate the combination 
with the average mass M = [M {Df ) + M{Dj))/2 closest 
to the Mo^ is kept. The combinatorial background has 
structure in M, but is more uniform in AM = M{Df) — 
M (Dj). The signal region is defined by \M-Md, \ < 12 
MeV and |AA/| < 30 MeV and the sideband region is 
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FIG. 22 Single tag yields for Da modes used in the 
CLEO-c analysis. Charge conjugate modes are combined. 
From [Alexander et ai, (|2008. ') . 
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FIG. 23 Double tag yields for Ds modes used in the CLEO-c 
analysis. The signal region is indicated by the rectangle in 
the center and the two sideband regions are the diagonally 
offset rectangles. There are 1089 double tag candidates in the 
signal region and 339 candid ates in the background region. 
From [Alexander" et al\ l|2008l '). 



defined by \M - Af^J < 12 MeV and 50 < |AM| < 140 
MeV. In this analysis the individual double tag yields 
and efficiencies are determined. The signal and sideband 
regions are shown in Fig. [23] 

All yields and efficiencies are combined in a likelihood 
fit to extract the Ds branching fractions. The branching 
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FIG. 24 The K" invariant mass near the 4> resonance in 
K~ K^n^ events from the CLEO-c double tag analy- 
sis. The single tag fit procedure used in the CLEO-c analysis 
is applied to extract the yield in each M{K~ K'^) bin, hence 
backgrounds are subtracted and the yields shown are for the 
D'^ K~ K'^-k'^ signal. The 4> resonance is clear above an 
additional broad component. Indicated in the plot are the dif- 
ferent mass windows consi dered by CLEO-c for t heir partial 
branching fractions. From [Alexander et al\ l|2008l ) 



fraction results from this fit is presented in Table IXVIII 
In addition to the branching fractions, CLEO-c deter- 
mines the number of D^D* pairs produced in their data 
sample to be Nd^d-^ = (2.93 ± 0.14 ± 0.06) x 10^. Com- 
bined with the luminosity, Cint — (298 ± 3) pb^^, they 
obtain the cross-section crDsD'.iEcm — 4.17 GeV) = 
(0.983 ± 0.046 ± 0.021 ± 0.010) nb, where the last sys- 
tematic is due to the uncertainty in the luminosity. 



CLEO-c does not quote a 



A7r+ branching frac- 



tion. The reason for this is that at the precision of this 
measurement the branching fraction for Df — > (f)Tr~^ is 
not a well defined quantity. Figure [Ml shows the K~^K~ 
invariant mass near the (j) resonance. The combination 
of the relatively broad 4> resonance and interference with 
other resonances, such as the /o(980), requires a com- 
plete amplitude analysis to determine the different con- 
tributions. Instead, CLEO-c provides partial branching 
fractions in different mass windows around the (jj reso- 
nance. These partial branching fractions, given in four 
mass windows centered at the <j) mass are pre- 
sented in Table [2£VIII| 

Systematic uncertainties from tracking efficiencies, tt'^ 
and reconstruction, and particle identification are 
common in this analysis to those from the analysis of 
the and I?+ absolute branching fractions discussed 
in Sect. IV. CI In additions, for modes containing an 77 
in the final state an uncertainty of ±4.0% is applied per 
rj. Other large systematic uncertainties in this analysis 
includes the uncertainties from the signal lineshape and 
the background parameterization in the fits for the yields. 
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TABLE XVII Branching fractions for Ds decays determined in the CLEO-c analysis. 



Mode Branching Fraction g (%) ^ A'+JC-7r+) Acp {%) 

B{D+ K'^K+) 1.49 ±0.07 ±0.05 0.270 ± 0.009 ± 0.008 ±4.9 ±2.1 ±0.9 



B{Dt ^ K+K--K+) 5.50 ±0.23 ±0.16 1 ±0.3 ±1.1 ±0.8 

B{Dt ^ K+K'Tv+n'-'') 5.65 ±0.29 ±0.40 1.03 ± 0.05 ± 0.08 -5.9 ±4.2 ±1.2 

B{D+ ^ K'^K'tv+tt+) 1.64 ±0.10 ±0.07 0.298 ± 0.014 ± 0.011 -0.7 ±3.6 ±1.1 

B{Dt ^ Tv+n~Tv+) 1.11 ±0.07 ±0.04 0.202 ± 0.011 ± 0.009 ±2.0 ± 4.6 ± 0.7 

B{Dt^TT+'n) 1.58 ±0.11 ±0.18 0.288 ±0.018 ±0.033 -8.2 ± 5.2 ± 0.8 

B{D+^TT+T]') 3.77 ±0.25 ±0.30 0.69 ± 0.04 ± 0.06 -5.5 ±3.7 ±1.2 

B(D+ K+TV+TT-) 0.69 ±0.05 ±0.03 0.125 ± 0.009 ± 0.005 +11.2 ± 7.0 ± 0.9 



TABLE XVIII Partial branching fractions in the mode — > 
K^K~tt'^ for events with a K^K~ invariant mass within 
AM MeV of the (p, \mif-ii+ - m^\ < AM. 



AM Partial Branching Fraction (%) 


5 


1.69 ±0.08 ±0.06 


10 


1.99 ±0.10 ±0.05 


15 


2.14 ± 0.10 ± 0.05 


20 


2.24 ± 0.11 ± 0.06 



These uncertainties are explored by using alternative fits. 

The CLEO-c analysis provides the to-date best de- 
termination of the hadronic branching fractions for Ds 
mesons. This analysis is statistics limited; the statisti- 
cal uncertainty in the Df — > K^K^n^ mode is 4.2% 
and the systematic uncertainty about 3%. The largest 
systematic uncertainties come from the yield extraction. 
Both the statistical and systematic uncertainties would 
improve with additional data. This analysis was based 
on 298 pb~\ CLEO-c has recorded a total of 589 pb^^ 
of data at this energy. 



TABLE XIX Summary of branching fractions for Cabibbo 
favored decays. 



Mode 




Branching Fraction 


Dt 






(1.49 ±0.09)% 


Dt 






(5.50 ±0.28)% 


D+ 






(5.6 ±0.5)% 


Dt 




K%K+Ti+Ti- 


(9.6 ± 1.3) X 10"^ 


Dt 






(1.64 ±0.12)% 


Dt 




K+K-TT+TT+n- 


(8.8 ± 1.6) X 10"^ 


Dt 




KlK%-K+-K+ir' 


(8.4 ±3.5) X 10"'' 


Dt 


^ 




(1.11 ±0.08)% 


Dt 




-K^ TV'^ 'K~ Tv'^ 


< 14% 


Dt 




■K^ 'K'^ ■k'^ ■K~ ■K~ 


(8.0 ±0.9) X 10"^ 


Dt 




TT ^ TT TT TT ^ TT ^ TT 


(4.9 ±3.2)% 


Dt 






(1.58 ±0.21)% 


Dt 






(2.5 ±0.9) X 10"^ 


Dt 




r)p+ 


(13.0 ±2.2)% 


Dt 






(3.8 ± 0.4)% 


Dt 


J, 




(12.2 ±2.0)% 



VII. CABIBBO SUPPRESSED DECAYS OF D", D+ , AND 
Dt MESONS 

A. Theoretical issues 



E. Summary of Cabibbo favored Dt decays 

The previous sections discussed the key measurements 
that established the absolute branching fraction scale for 
Dt meson decays. These measurements have evolved 
from model dependent determinations, e.g., making use 
of equal semileptonic widths as for the D^ decay, to 
model independent measurements using tagging tech- 
niques. Also as the measurements have become more 
precise we need to be more precise about what is mea- 
sured. For example, the often-used normalization mode 
Dt 07r+ suffers from a contamination from the Df — > 
/o(980)7r+ under the 4>tt^ signal. 

The results for the Cabibbo favored modes are sum- 
marized in Table IXiXl 



Studies of hadronic singly Cabibbo-suppressed decays 
of charmed mesons are important for several reasons. 
First, these decays hold the potential for future obser- 
vati on of direct (i.e. not associated with D'^D mix- 
ing (|Petrovl . [200l ) CP violation in the Z^-system. In the 
Standard Model, this is due to the fact that the final state 
particles contain at least one pair of a quark and anti- 
quark of the same fiavor, making possible a contribution 
from penguin-type amplitudes. Those amplitudes pro- 
vide an access to the third generation of quarks (&-quarks 
in the loops), needed for observation oi CP violation in 
the S tandard Model ( Bianco et aLl . 120031 : iBuccella et all . 
Il993f ). New Physics can also make an entrance through 
those transitions, affect ing both the amplitude s and CP- 
violating asymmetries ( Grossman et all l2007t ). Second, 



it offers new ground for studying strong dynamics in 
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hadronic decays, in particular, the issue of flavor SU{S)f 
breaking in D-decays. For example, one of the famous 
failures of the applications of SU{3)f symmetry involves 
the prediction that the decay rates for D'^ K'^K~ 
and — > Tr~^iT~ are equal. In reality, the first rate 
is about three times larger than the second one. Other 
puzzles include the fact that the rates for decays like 
Z?"*" — > K*~^K are so much enhanced by strong dy- 
namics that their values appear to be as large as the 
ones of Cabibbo-favored decays. One popular expla- 
nation for s uch phenomena include resonant final state 
interactions ( Chau and Cheng| . ll989t[Kamal and Vermal . 



|1987) that affect not onl y J-decays, but also D^D" 
], 119991: iGolowich and Petrol Il998^ 



mixing (,Falk et al. 
Ther e are also oth er explanations ( Chau and Cheng . 
Il992l : ISavagd . Il99l[ ). In order to study those phenom- 
ena it is convenient to select a base formalism for studies 
of hadronic transitions. 

It is convenient to use the topological dia- 
gram approach to predict unknown branching ra- 
tios for singly-Cabi bbo-suppresse d decay s. The 

re- 
and 



Chiang all (l2003f) an d 



analysis, d one in 

peated in iBhattacharya and Rosneil (|2008bl ) 
iBhattacharva and Rosnen ( 2008al ) with updated ex- 
perimental data, is displayed in Tables IXXI through 
IXXIIII 



SA' . The amplitudes C and E extracted from Cabibbo- 
favored charm decays imply values of C — X'C and 
E' = \' E which may be used in constructing amplitudes 
for singiy-Cabibbo-suppressed D'^ decays involving rj and 

i- 



2. D ^ PV transitions 

A similar technique can be applied to describe 
D PV transitions. In this case, similar 

topological amplitudes are denoted by a subscript 
"V". We present the most r ecent result in Table 
IXXIII ([Bhattacharva and Rosneil . [2008aD . 



B. Cabibbo suppressed D° and decays 

Experimentally, Cabibbo suppressed or doubly 
Cabibbo suppressed decays of £>" or £)+ mesons are 
almost always measured relative to a Cabibbo favored 
normalization mode. This includes most CLEO-c anal- 
yses as the branching fractions for Cabibbo suppressed 
modes are typically suppressed by iKd/KsP ~ 0.05 and 
the statistics in these modes using a double tag analysis 
would be limited. In some cases, e .g. the C LEO-c 
analysis of KK final states ( Bonvicinil . |2008), 



\. D ^ PP transitions 



CLEO has normalized against the number of produced 
DD events and measured directly the branching fraction. 



Topological diagram approach to singiy-Cabibbo- 
suppressed transitions can make use of the information 
obtained from the fits of CF decays discussed above. In 
particular, the ratio of primed (SCS) to unprimed (CF) 
amplitudes is fixed, it is just A' = ian Oq = 0.23. Ta- 



ble IXXI (taken from Bhattacharya and Rosnerl ( 2008131 )) 
presents the most recent compilation of the branching 
ratios, amplitudes, and representations in terms of re- 
duced amplitudes for singiy-Cabibbo-suppressed (SCS) 
charm decays involving pions and kaons. The extracted 
topological amplitudes, in units of 10^^ GeV, are 



T' 
C 
E' 
A' 



6.44 ; 

-4.15- 2.25i 
-1.76 + 3.48i 
0.55 - 1.14i . 



(48) 



The deviations from flavor SU(3) in Table IXXI are dis- 
cussed below. ^ 

Note that the decay _D'^ K^K is forbidden by 
SU{?>)f- Estimates of S'[/(3)F-breaking effects lead to 

predictions for B{D^ K^K ) that are consistent with 
experimental observations, but are by no means reliable 
(iDai et all Il999l : lEeg et all l200ll: iLipkinl . Il980l: iPhaml . 
Il987l ). We shall discuss those below. 

Final states with rj and 77' require additional consider- 
ations. In particular, new topological amplitudes, flavor- 
singlet singlet-exchange SE' and singlet-annihilation 



1. Two-body decays of D*' and D"*" 

There is a substantial amount of data on the two- 
body decays of and D+. The first measurements 
of Cabibbo suppressed 13° decays were for 13° 
K~ K' ^ and D° .^""""^ d one by the Mark II experi- 
ment ( Abrams et aLl . [l979bl ). Since the first observation 
of these modes they have been measured by many exper- 
iments with increased precision. In these measurements 
the K^K^ and 1?° 7r^7r+ branching fractions 

are measured relative to the 13° K~7r~^ yield. Ex- 
periments operating above the cc threshold tag the D'^ 
by looking at the D^-D*~^ mass difference in the decay 
D*+ D°n+. 

The results for the R - K+ an d -> TT-n+ 

decays are summarized in Table IXXIIII The most p re- 
cise measurement is that of CDF (|Acosta et aLl . [20051 ). in 
the — > K^K^ they reconstruct about 16,000 signal 
candidates. 

As can be seen from Table IXXIIII the rate for D — > 
KK is larger than the rate for D ttW by a factor 
of three. In the SU{3)f (or in the [/-spin) symmetry 
limit, those rates should be the same. SU{3)f is, in 
general, expected to work to 30%, so this is a rather 
severe violation of this symmetry. 

While the one popular explanation for this puzzle in- 
volves final state interactions (i.e. a presence of a reso- 
nance that couples stronger to K'^K~ compared to 7r"'"7r~ 
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TABLE XX Branching ratios, amplitudes, decomposition in terms of reduced amplitudes, and predicted branching ra- 
tios for singly-Cabibbo-supprcsscd charm decays involving pions and kaons. Predictions for the branching ratios are 
from (jBhattacharva and Rosner . 2008b|). 



Meson 


Decay 


B 


P* 


1-41 


Rep. 


Predicted 




mode 


{10-') 


(MeV) (10"'' GcV) 




B {W-'^) 






1.40±0.02 


921.9 


4.61±0.03 


-{T' + E') 


2.23 




TT TT 


0.80±0.08 


922.6 


3.49±0.17 


-iC'-E')/V2 


1.27 




K+K- 


3.93±0.07 


791.0 


8.35±0.08 


(T' + E') 


1.92 






0.37±0.06 


788.5 


2.57±0.35 








D+ 


+ 
TT TT 


1.24±0.07 


924.7 


2.73±0.08 


-{T' + C')/V2 


0.87 




K+K° 


6.17±0.20 


792.6 


6.58±0.11 


T' ~A' 


5.12 


Df 




2.44±0.30 


915.7 


5.84±0.36 


~(T' -A') 


2.56 




7V°K + 


0.75±0.28 


917.1 


3.24±0.60 


^iC' + A')/V2 


0.87 



TABLE XXI Real a nd imaginary parts of amplitudes for SCS charm decays involving rj and 77', in units of 10 ^ GeV from 
(|Chiang et am2003l ). 



Amplitude 



Expression 



Re 



Im I A 



1O 



3\/2 



n"7j) 2E' + SE' 0.63 

7r"ry') ^{C + E') + SE' -2.95 

^777?) C' + SE' -4.14 

-2.10 
-0.75 
2.92 
1.85 

^AjPt ^ v'K+) \{2T' + C' + 3A') + SA' 2.59 



C' + SE' 

^A{D° -> r^r^') \{C' + &E') + SE' 
^/iA{D+ TT+Tj) T' + 2C" + 2A' + SA' 
-^A{D+ ^n+ri') ^{T' -C + 2A')+SA' 
VSAiDt 'nK+) -{T' + 2C') + SA' 



9.21 7.79 ±0.54 
0.62 3.54 ±0.35 
-2.25 5.91 ±0.34 

2.66 3.48 ±0.38 
-6.77 8.21±0.26 
-0.01 3.72±0.15 

4.50 8.05±0.88 
-1.41 3.43±0.57 



state), it might be t emptiiiK to try to under stand the is- 
sue in factorization ( Chau and Chend . ll992l ). Neglecting 
for a moment the annihilation diagram contribution, 



A 



KK 



Ik m% 



(49) 



With the recent lattice evaluations /k/Ztt = 1.218 ± 
0.002j^n noi from a rec ent lattice QCD ca lculation with 



domain-wall fcrmions (jBeane et all |2007| ). assummg a 



modified pole dominance for the form- factors F (mj^) 
and F^''^ (m j), and extrac t ing th em from semileptonic D- 
decays ( see jArtuso et all 12008^) for a recent review and 
( Bessonl . l2009l ) for recent determination of parameters). 



we get 



A — ■ 



1 .32^4^^ . 



(50) 



In other words, factorization predicts about 30% 
breaking of SU ( 3)f in spectator amplitudes (c.f. 
( Chau and Chend . ll992h V Clearly, this is not sufficient 
for the resolution of the puzzle. The presence of final 
state interaction (FSI)-enhanced exchange amplitude is 
crucial for the explanation of this phenomenon. 

A number of other two body final states to pseudo 
scalars and have been studied. These decays are summa- 
rized in Table |2D£IV| 



The most complete study of D mesons decays to fi- 
nal states_conteiinngj7 and 77' mesons is done by CLEO- 
c ( Artuso et all |2008| ). The excellent electromagnetic 



calorimeter and the clean environment near threshold 
combined with the large data sample collected at the 
■0(3770) has allowed CLEO-c to measure many modes 
not previously seen. 

This analysis uses 281 pb^^ of data collected at the 
-0(3770) resonance. In this study CLEO-c makes use 
of single tags; the modes studied here have sufficiently 
small branching fractions that using D tagging is not 
useful. The tt*^ and rj mesons are reconstructed in the 
77 final state. In addition, for modes with two r; mesons 



in the final state {rjr] and 7777') the rj 



— > TT ' TT TT 



chan- 



nel is used to reconstruct rj mesons. The 77' is recon- 
structed in the channel 77' t^tt^tt^ . It is required that 



402 < M^^+^- -Mr,< 418 MeV/c2. 

The yields are extracted by fitting the Mbc distri- 
butions after selecting events consistent with AE = 0. 
In Figs. the observed signals are shown. The 

significance for all modes are over 4ct except for the 
D'^ — > rfTT^TT^ mode where the significance is estimated 
to be 3.2(7. The observed yields and branching fractions 
are summarized in Tab. IXXVI This data makes it pos- 
sible to constrain new singlet exchange SE' amplitudes 
introduced in Sect. IVII.A.ll In order to do that, one can 
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TABLE XXII Branching ratios and invariant ampl itudes for singly -Cabibbo-suppressed decays of charmed mesons to one 
pseudoscalar and one vector meson (from ( Bhattacharva and Rosnerl . [2008a l 1 . 



Meson Decay 
mode 



Representation 



B fAmsler et al.. 2008) p* 

(%) (MeV) 



n'^ p 

K+ K*- 
K- K*+ 

^0 1 



VP 
rjuj 

V<t> 
•q 'uo 



~{Tv' +Ep') 
-{Tp' + Ev') 
\ {Ep ' + Ev ' — Cp ' — C'v ') 

Tv' +Ep' 

Tp ' + Ev' 

Ey ' — Ep ' 
Ep ' — Ey ' 

2{Ep + Ev ^Cp + Cv ) 
-^{2Cv ' — Cp ' — Ep ' — Ev ') 
-^(2Cv' + Cp' + Ep' + Ev') 

-^{Cp ' — Ep ' — Ev ') 
^{Ep' + Ev' + Cp' + Cv') 
(Ep ' + Ev ' + Cp ' — Cv ') 



0.497±0.023 
0.980±0.040 
0.373±0.022 
0.153±0.015 
0.441±0.021 

< 0.18 

< 0.09 
Q.124±0.012 



7^ 



K*°K+ 

7V°p+ 



2v/3 



1-41 

(io-«) 



763.8 
763.8 
764.2 
609.8 
609.8 
605.3 
605.3 
644.7 
761.2 
652.0 
488.8 
648.1 
342.5 
333.5 



1.25±0.03 
1.76±0.04 
1.08±0.03 
0.97±0.05 
1.65±0.04 



0.81±0.04 



A. 



Cp'~Tv') 



VP 
V'P+ 



-,iAp' + Av' + Cp' + Tv') 
Cp' 

[Tv '-Av') 
^{Av' -Ap' -Cv' -Tp') 
^{Av' + Ap' + 2Cv' +Tp') 

75 (^-' 



Dt TT+K*'^ 

n°K*+ 
riK*+ 

K%+ 
K+ p° 
K+io 

K+ (j> 



~Av' 
(Tp' - 



-Ap' 
Ap') 



Tp') 



0.082±0.015 

< 0.034 
0.620±0.070 
0.435±0.048 

< 0.7 

< 0.5 
3.18±1.38 



767 
764 
647 
611 
767 
656 
349 
612 



0.32±0.03 

1.13±0.06 
1.03±0.06 



2.78±0.60 



-Tv') 

.V ' + Av' 



{Av' 

MTp' + 2Cv' +Ap' -Av') 
^{2Tp' + Cv' + 2Ap' + Av') 
{Ap' -Tp') 
-^^[Cp' + Ap') 
-^^[Cp'-Ap') 
Tv' + Cp' + Av' 



0.225±0.039 



0.27±0.05 
< 0.057 



773 
775 
661 
337 
743 
745 
741 
607 



0.79±0.07 



0.92±0.09 



TABLE XXIII Measurements of D° R- K+ and D° TV n . The branching fractions have been recalculated using 
B{D° K-n+) = (3.89 ± 0.05)%. 



Experiment g(73" ^ K' K+) (10~^) BjD" 7r~7r+) (10-^) 



CLEO-c CBonvicini, 2008: Rubin et al. 2006) 


4.08 ±0.08 ±0.09 


1.41 


±0.04 


±0.03 


BES II (Ablikim et al. 2005) 


4.75 ±0.43 ±0.17 








CDF (Acosta et al. 2005) 


3.859 ± 0.043 ± 0.069 


1.40 


±0.02 


±0.03 


FOCUS (Link et al. 2003) 


3.863 ± 0.054 ± 0.074 


1.37 


±0.05 


±0.03 


CLEO II fCsorna et al. 2002) 


4.05 ±0.13 ±0.13 


1.36 


±0.06 


±0.07 


E791 fAitala et al. 1998a) 


4.24 ±0.12 ±0.13 


1.56 


±0.08 


±0.12 


CLEO II fAsner et al.. 1996) 


4.51 ±0.27 ±0.28 








E687 fFrabctti et al. 1994a) 


4.24 ±0.27 ±0.35 








E691 (Anios et al. 1991) 


4.16 ±0.39 ±0.39 








CLEO (Alexander et al. 1990b) 


4.55 ±0.39 ±0.28 









Average 3.98 ± 0.07 1.40 ±0.03 
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TABLE XXIV Measurements of D° and D+ decays to 
Cabibbo suppres sed, non-stranee, tw o-body final states. The 
averages are from lAmsler et aZl l|2008f ). the decay D° — > ■k'^tt~ 
is discussed in Sect. IVII.B.ll 



Mode 




B (10-^) 






1.40 ±0.02 


D"^ 




0.80 ±0.08 






0.57 ±0.14 


D"^ 




< 0.26 


D+ ^ 




1.24 ± 0.07 


D+ ^ 




3.39 ±0.29 


D+ ^ 




< 0.34 



TABLE XXV Measurements of D meson decays to final 
states with r) and r)' mesons. 



Mode 



B{D+ 
B{D+ 
B{D° 
B{D° ^ 
B{D'' ^ nv) 
B{D'' W) 
B{D^' rin+T7-) 
B{D+ 777r+7r°) 
B{D° T]'Tv+n-) 

B{D+ VTT+Tr") 



Yield Branching Fraction (10 



1033 ± 42 
352 ± 20 
156 ± 24 

50 ±9 
255 ± 22 

46 ±9 
257 ± 32 
149 ± 34 

21 ±8 

33 ±9 



34.3 ± 1.4 ± 1.7 
44.2 ±2.5 ±2.9 

6.4 ± 1.0 ±0.4 
8.1 ± 1.5 ±0.6 

16.7 ± 1.4 ± 1.3 

12.6 ±2.5 ±1.1 
10.9 ± 1.3 ±0.9 

13.8 ±3.1 ± 1.6 

4.5 ±1.6 ±0.5 

15.7 ±4.4 ±2.5 



rewrite four equations for D'^ decay amplitudes to the 
final states with 77*^'-' : 

- V6AiD" ^ riT:°) = 2E'-C' + SE', 



1630108-002 



2V2 
3V2 



i (E' + C) + SE', 

c ± 5"^;', (51) 

1 
7 



(C + QE') + SE'. 



It is interesting to note that the right-hand side of each of 
Eqs. (l51|l determines a vector in a complex plane. Since 
both amplitudes and phases of C and E' are known from 
Eq. (j48p , these four equations contain a common complex 
off-set, SE' . Since only the magnitudes of the right-hand 
sides of these equations are known, they each define a cir- 
cle in the complex plane with the radius given by that 
magnitude. Plotting them on the same graph then de- 
termines SE' . 

This is done in Fig. [^H Notice that all circles inter- 
sect in two points, which determine two possible solu- 
tions for SE'. The smaller values for SE' = (-0.7 ± 
0.4) X 10-^ GeV ± i(-1.0 ± 0.6) x 10"^ GeV are the- 
oretically preferable, as SE' is an Okubo-Zwcig- Iizuka 
( O ZI)-suppr cssed amplitude (jlizukal . .1966:.Okubol . ll977l: 
IZweiel[r96l . 



2. Multi-body decays with kaons and pions 

Multibody decays of and £)+ mesons were also ex- 
tensively studied. While theoretical studies of those tran- 
sitions are limited, some of those decays can be used in 



the D alitz-plot analyses of D'^D mixing ( Artuso et al\ . 
120081 ). The measurements of branching fractions to fi- 
nal states with three of more pions, including final states 
with r] and uj mesons could be found in Table IXXVII In 
Section HXl Dalitz plot analysis of three-body final states 
are discussed. 

In addition to the D° KK decays discussed above, 
many other Cabibbo suppressed decays with two kaon in 




1.86 1.E 

M,, (GeV/c2) 



FIG 



25 



Yields for a) 
D" ^ ?77r°, d) £1° ^ 
From I Artuso et al\ (fm 



D+ 

,0 



b) 



D+ 

, and D' 



r] TV 

,0 



TABLE XXVI Measurements of D" and D+ decays to 
Cabibbo suppressed final states with three or more pions in 
the final states. Final states with 77 and lo mesons are also 
included. Limits are given at 90% C.L. 



Mode 






B (10-^) 




TT^ 




14.1 ±0.6 


D"^ 






7.44 ±0.21 






- + - 

n TV n n 


4.2 ±0.5 




n'' 


n n 


< 0.35 








(1.00 ±0.09) 




TV^ 




4.2 ± 1.2 


D+ ^ 


7V' 




3.21 ±0.19 


D+ ^ 


TV' 




4.6 ±0.4 


D+ ^ 


7V' 




1.14 ±0.08 


D+ ^ 


TV' 




1.63 ±0.16 
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(a) DO - 

L Jl n n 


-> r|7i+7t" 







(c) DO r\'n*n- 



1.83 



1.86 




1.89 1.83 



1.86 

M,, (GeV/c2) 



1.89 



FIG. 26 Yields for a) D° r]-!v+TT~, b) D+ t^tt+tt", c) 
p° - » Vtt+tt", and d) D+ ^ ?7'7r+7r". From lArtuso et al\ 
(|2008l ). 



the final states have been studied. Dalitz plot analysis 
has been performed on some three-body final states as 
discussed in Sect. IIXI The final states with two kaons in 
the final state are summarized in Table IXXVIII 

TABLE XXVII Measurements of D° and D+ decays to 
Cabibbo suppressed final states with two kaons. 



Mode 




B (10-^) 








3.5 ±0.5 


D° 






2.7 ±0.5 


D° 






3.29 ±0.14 


D° 






< 0.59 


D'> 




K+K-Tv+Tv- 


2.43 ±0.012 






KgKgTV TY 


1.30 ±0.24 


D° 




KI;K-TT + TV + TV- 


< 0.15 


D° 




K+K-Tv+Tv-n" 


3.1 ±2.0 


D+ 




K+K-TT+ 


9.63 ±0.31 


D+ 






1.69 ±0.18 


D+ 


J 




2.32 ±0.18 


D+ 




K+K-TT + TV + TT- 


2.3 ± 1.2 



C. Cabibbo suppressed Ds decays 



1630108-005 




5 10 

Re (A) (10-^GeV) 



FIG. 27 Graphical representation of Eqs. (|5ip used to de- 
termine SE' from t he amplitude analysis for D rj^''^ r){-K^) 
lArtuso et al\ (|2008l ). Circles represent absolute values of the 
decay amplitudes. The intersection points provide two possi- 
ble solutions for SE' (see text). 



The Cabibbo suppressed Ds decays are final states 
with one or three kaons. The measured decays are listed 
in Table IXXVIIII This table also includes the doubly 
Cabibbo suppressed decay Z?+ K~^K^it~. CLEO- 
c ( A dams et al, ,200711 has performed a systematic study 
of two-body Ds decays. 



D. Doubly Cabibbo suppressed decays 

The doubly Cabibbo suppressed decays have two 
Cabibbo suppressed weak couplings. Naively, the rate 
for the doubly Cabibbo suppressed decays are supressed 
by a factor of iduti^Oc ~ 2.8 x 10~^. The first observa- 
tion of a doubly Cabibbo su ppressed decay was in the 
decay channel D° K+n^ (jcinabro et'aU . \l994) . Ex- 
perimentally, the fiavor, Z?" or D°, of the initial state 
is tagged by the charged of the slow pion in the decay 
of a D*'^ D'^tt^. The simplest measurements observe 
the time integrated rate of decays and do not sep- 
arated direct decay contributions from mixing, where a 
D'^ oscillated to a l)° and decayed via a Cabibbo favored 
decays. 

The doubly Cabibbo suppressed decays that have 
been studied are summarized in Tabic IXXIXI The three 
most pre cise measurements of the K^tt~ decay 

20081). BAB AR (jAubert et al\ . 



by CD F (jAaltonen al.,.,^^^^,. ^.^^^ ^^^ 
l2007b[ ). and Belle (|Zhang et all l2006t) obtained branch 
ing rations with respect to K~Tr~^ of (4.15±0.10) x 

10-3, (3.53±0.08±0.04) x lO^^, and (3.77±0.08±0.05) x 
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TABLE XXVIII Cabibbo suppressed D+ decays. 



Mode 


Dt^ 


K+n" 


Dt^ 








Dt^ 


K+r/ 


Dt^ 


K+Tv+n- 


Dt^ 


ifgyr+Tr+TT 




K+K+K- 




K+K+TT- 



ReL 

(Adams et al. 2007) 
(Adams et al. 2007 : Link et al. 2008) 
(Adams et al. 2007) 
(Adams et al. ^01) 
(Alexander et al, 2008) 
(Xink et al, 2008) ' 
(Link et al. 2002) 
(Link et al. 2005b) 



B/10~^ 
0.82 ±0.22 
1.25 ±0.15 
1.41 ± 0.31 
1.6 ±0.5 
6.9 ±0.5 
3.1 ± 1.1 
0.49 ±0.17 
0.29± 0.11 



10^'^ respectively. The agreement between these mea- 
surements is not very good, the PDG apphes a scale fac- 
tor of 3.3 for the error on their average to obtain the aver- 
age ratio of branching fractions to be (3.80±0.18) x 10^'^. 



of branching fractions B{D+ K+K+K-)/B{D+ 
K-TT+1T+) = (9.49 ±2.17 ±0.22) x 10"''. This gives the 
branching fraction B{D+ K+K+R-) = (8.7 ± 2.0) x 

10-^ 



TABLE XXIX Doubly Cabibbo suppressed D" decays. The 
first column (B) shows the branching fraction for the decay 
and the second column (R) shows the ratio of the branching 
fraction with respect to the corresponding Cabibbo favored 
decay. 



Mode 




B (10--*) 


R (10-') 




K+Tv- 


1.48 ±0.07 
3.05 ±0.17 


3.80 ±0.18 
2.20 ±0.10 

O r)O + 0.25 

'J-^<J_0.22 



The decay D° 



K~^T T 7r° w as first observed by 
CLEO (|Brandenburg et al\ . l200lh . The PDG average 
is domin ated by the more rec ent measure ments from 
BAB AR (|Aubert et all l2006bD and Belle (|Tian et all . 
I2005D . 

The first significant D° K+tt~tt~^tt~ observation 
was made by CLEO (jPvtman et all l200lh . The most 
recent an d precise rr i easure ment of this decay was done 
by Belle (|Tian et all l2005l). 



Both CLEO-c 



BABAR (jAubert et all 



IPvtman et all l2006h and 
has studied the doubly 
Cabibbo suppressed decay D+ K+n^ . CLEO-c has 
reconstructed candidates in a 281 pb^^ sample of e+e~ 
data recorded at the -0(3770). BABAR has used a sample 
of 124 fb-i recorded at the T(4S'). CLEO-c and BABAR 
finds branching fractions in good agreement with each 
other, B{D+ K+t:^) = (2.24±0.36±0.15±0.08) x lO^'^ 
and B{D+ K+ir^) = (2.52±0.46±0.24±0.08) x lO""*, 
respectively. The average branching fraction obtained is 
(2.37 ±0.32) X lO-'^. 

The fin al state D+ — > K~^tt+ tt~ has been studied 
by E 687 (iFrabetti et all ll995bD . E791 (|Aitala all 
I1997D . and FOCUS (iLink et aLl . l2004bD . The average 
branching fraction from these measurements is B{D+ — > 
K+TT+TT-) = (6.2 ±0.7) X 10-*. 

The de cay D+ — » K~^K ~^K~ has been observed by 
FOCUS (jLink et all 12002^ . They measure the ratio 



VIII. FINAL STATE INTERACTIONS AND AMPLITUDE 
ANALYSIS 

One of the simplest ways to analyze decays of D- 
mesons is to employ the flavor flow diagram technique 
described earlier. One potential problem with the appli- 
cation of this technique^ to charm decays involves assign- 
ment of quark amplitudes (T, A, etc.) to a particular 
decay. The root of the problem involves inelastic final 
state interactions (FSI). 



A. Hadronic decays into meson states 

Historically, the issue came up with decays of the type 
4)K^ , which have been claimed to originate entirely 
from quark exchange amplitudes. Thus, in the topologi- 
cal SU{Z) or flavor-flow analysis of this transition only an 
exchange amplitude E should be assigned to this decay. 
However, FSI contribution of the type 
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cf,K° 



(52) 



could proceed through the color-suppressed internal W- 
emission diagram C followed by strong-interaction rescat- 
teving r]^'^ K'^* (pK^ . This contributi on is not op t ional , 
but is, in fact, req uired by unitarity ( Donoghug . Il986l : 
iFaifer all Hool)- While in the example above par- 
tial cancelation o ccurs between the intermediate rjK^* 
and rj'K^* states ( Lipkinl . 119871 ). this cancelation is not 
generic. If large, the contributions of this type could be 
important in the top ological flavor -flow amplitude anal- 
ysis of charm decays (IChend . [20031) . 



^ Similar problem could affect charm decay analysis using factor- 
ization approximation. 
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Eq. I|53p . dominate the branching ratio for D° KsKs- 
There is, however, no reUab le way to estimate those (see, 
however, ( Eeg et all . l200lh '). Another way would be to 
accept that this, and similar branching' ra tios a r e dom- 
inate d by final state interactions (Li pkin| , llQSOl : iPhaml . 
Il987t ). A simple two-channel model estimates give 



FIG. 28 The two quark diagrams that contribute to the 
decay D° ^ T^gi^g. Since Vcd = —Vus the two amplitudes 
represented by these diagrams largely cancel. In the limit that 
the d and s quark masses where the same the cancellation 
would have been exact. 



One way to study the importance of inelastic FSI con- 
tributions in charm decays is to seek guidance from ex- 
perimental studies of " annihilation" decays, i.e. decays 
whose contribution is dominated by weak annihilation or 
exchange amplitudes in the topological flavor-flow anal- 
ysis. 

Another related decay mode that is interesting from 
this perspective is D° KsKs- Naively, there are two 
W exchange diagrams that contribute to this final state 
as illustrated in Fig. [281 Since Vcd — —Vus, these am- 
plitudes interfere destructively, so in the flavor SU{3)f 
limit the branching ratio for this process is zero. Thus, 
in addition to being the "pure annihilation" decay, the 
rate of D'^ KsKs transition explicitly probes SU{3) p- 
breaking corrections. It should be rather small. 

Interestingly enough, a naive calculation of this decay 
rate in factorization gives exactly zero. 



r(i:>° 



r(i?°^X+A'-)tan2 -((5„-<5i) 



^ KsKs) 



1 



foPD ■ {pk'> -Pk°) = ^1 (53) 



so Sfact(-D" — > KsKs) — 0. As we discuss later in this 
section, experimental analyses of this transition, how- 
ever, clearly yield a non-zero result. 

The ratio of branching fractions 
B{K^K°c}/B(K^tt±jtZ) has been measured b^ 



£>(AgA y)/D(Ag7r' TT ) nas been mcasure a Dy 
CLEO (Alexander eToZI. Il990bl), E 687 (Frabct ti et all . 
Il994dl. CLEO II (lAsner et all I1996D. and FO- 



CUS (iLink et al\ . 1200531) . CLEO-c ( iBonvTchdI . I2008D 



has studied this decay using a single tag technique 
and normalized to the number of D^D^ events pro- 
duced. These measurements are summarized in 
Table IXXXI Measurements of the branching ra- 
tios B{K^K^)/B{K^TT^'K~) has been rescaled using 
BiKl-K+TT-) = (2.99 ±0.17)% (|A^sieLiIaZl,|200i). 

The most recent, and most precise, measurement from 
CLEO-c gives the smallest central value. Given the large 
uncertainties in the earlier measurements there is no 
strong inconsistency between the different measurements. 
This clearly points to shortcomings of factorization cal- 
culation outlined above. 

One way to understand this branching ratio would 
be to assume that non-factorizable pieces, dropped in 



1 



(54) 

where 6q and 5i are the phase shifts for / = and 1=1 
amplitudes. Estim ates with other models of FSI give 
comparable results ( Dai et all . Il999h . While these esti- 
mates are by no means reliable, they serve as an indica- 
tion of importance of FSI in charm hadronic decays. 



B. Baryonic decay Dt — » p^n 

Final states with baryons are not possible for the 
and . The lightest neutral final state pp has a mass of 
1876.54 MeV/c^ and is just above the D° and D+ mass. 
However, the is kinematically allowed to decay to 
p+n. This decay is also quite interesting because the 
flavors of all valence quarks that constitute the initial 
state (cs) differ from the flavors of the flnal-state quarks 
composing the p^n pair. Thus, it is quite tempting to 
declare that the transition Dt p'^n proceeds o nly via 
the weak annihilation graph ( Chen et al\ . l2008t iPhaml . 
il980a..b). 

A factorization ansatz can be employed in order to es- 
timat e the branching ratio for this process ( Chen et al\ . 
l2008h . It must be emphasized again that contrary to 



hadronic B-decays, simple factorization has not been 
proven in charm transitions, especially as applied to an- 
nihilation amplitudes. Nevertheless, a factorized decay 
amplitude is 

MDt Pn) = -j=VcsV*^aifD,PD, {pn\u-ft, (1 - 75) d\0), 

(55) 

where po^ — Pp + Pn is the four-momentum of a Ds- 
meson. The matrix element between the vacuum and 
the final state can be parametrized. First, let us note 
that vector current conservation implies that 

p'^^{pn\u^f^ (1 -75)d|0) = {mp + mn){pn\u'j5d\0), (56) 

so the decay amplitude can be parameterized as 

A{Dt^pn) = ^VcsV:^mfD^ 



V2 



X ( 2mNgl"' + ^^^33" ) Up75«n(57) 



where " are the formfactors parameterizing the 
baryon current, and mjv is the nucleon's mass. The 
two formfactors g^^ and (jf" can be related to each 
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TABLE XXX The observed branching fractions for D" 
branching fraction K%ti'^-k~ when used. 



K%K%. The errors are statistical, systematic, and from normalization 



Experiment 



events BiD" A'gA"g)/10" 



CLEO-c (Bonvicini, 2008) 
FOCUS (Link et al. 2005a) 
CLEO II (Asner et al. J^96) 
E687 (Frabctti et al, 1994c) 
CLEO (Alexander et al, 1990b) 



68 ±15 1.46 ±0.32 ±0.09 
79 ± 17 4.31 ± 0.96 ± 0.48 ± 0.24 
26 3.02 ±0.66 ±0.48 ±0.17 
20 ±7 11.7 ±3.9 ±3.9 ±0.7 



6.3 



+3.3 



±0.6 ±0.4 



Average 



1.93 ±0.30 



other (|Chen et all 120081: iPhaml . Il980al lbl) , 



^^5f(pk), (58) 



so that the decay amphtude takes the form, 



V2 



( — ^ ) 3l" Upj5Vn. (59) 



This amplitude leads to the estimate of the decay branch- 
ing r atio B(Dt pn) in the factorization approxima- 
tion (|Chen et a/.l . l2008l) . 



10- 



(60) 



The theoretical error quoted in Eq. ([60|1 is entirely 
due to t he uncertainty in the form-factor value of 
(|Chen et aLl . l2008l ). which was obtained by ex- 
trapolation of the nucleon data with a particularly as- 
sumed shape of g'^-dependence. This estimate gives a 
rather small branching ratio, which nevertheless can be 
tested experimentally. CLEO-c has studied this final 
state (|Athar et a/.l . l2008l ). 

As (anti-)neutrons are hard to reconstruct, CLEO-c 
uses a missing mass technique to identify this signal. All 
particles in the event, except for the (anti-)neutron, is 
reconstructed and the signal is extracted by looking in 
the missing mass distribution of the events, which for 
signal will peak at the neutron mass. 

CLEO-c uses 325 pb~^of e+e" annihilation data col- 
lected at a center-of-mass energy of 4170 MeV. At this en- 
ergy pairs of D+D*~ and D~Dl+ are produced. CLEO-c 
uses 8 tag modes {K~^K~Tr~ , KgK~, r]TT~ , r]'TT~, (j)p^ , 
7r~7r"'"7r~, K*~ K*^ , and rip~) to first reconstruct a 
candidate. It is required that this Dg candidate has a 
reconstructed invariant mass which is within 2.5cr of the 
known Dg mass. Next, this candidate is combined with 
a photon. The recoil mass squared against the Dj + 7 is 
calculated and required to be consistent with the mass of 
the Ds- Note here that it does not matter if the photon 
came from the D* that is the parent of the D~ or from 
the parent of the other Dg in the event. This missing 



mass squared distribution is fit to determine the number 
of tags, CLEO-c reports finding 16,995 Ds tags. This 
yield will be used as the denominator in the branching 
fraction calculation. 

CLEO-c then search for the proton candidate in the 
momentum range from 150 to 550 MeV/c. In this mo- 
mentum range CLEO-c uses dE/dx to identify the pro- 
ton, 550 MeV/c is below Cherenkov threshold. Kine- 
matic fits are performed to the , photon, and proton 
candidate. As it is not known if the photon came from the 
D*~ Dj-f decay or the Dfj decay, two differ- 

ent fits are performed. First the photon is added to the 
Dj to form a D^^ . The I?*^ momentum is constrained 
to the known momentum from the e^e^ —^ D*~ D'^ reac- 
tion and the D*~-Dj mass difference. For the other hy- 
pothesis the Dj momentum is constrained to the known 
momentum from the e"'"e~ — > D*~ Df reaction and the 
proton is combined with the missing momentum of the 
event to make a Df candidate, add the photon and con- 
strain the Dl'^-Df mass difference. The combination 
with the lowest is selected. Based on Monte Carlo it 
is estimated that the right combination is selected 95% 
of the time. Applying these kinematic constraints im- 
prove the resolution on the missing mass by a factor of 
two. In addition cuts are applied on the to reject 
combinatorial background not consistent with the signal. 
If there are multiple photon candidates in an event the 
combination with the lowest overall is selected. 

In Fig. [blithe distribution of the recoil mass against the 
proton shown. There are 13 candidate events consistent 
with the — s- pn signal. From this yield, the number 
of tags, and the efficiency for reconstructing the proton 
CLEO-c determines the branching fraction 



BiDt -> pn\ 



(1.30 ±0.36 



+°;-) X 10-^ 



(61) 



This result shows quite unambiguously that the 
factorization-ansatz estimate of Eq. (pO)) fails by more 
than three orders of magnitude! This could be because of 
the following two reasons. First, the use of a factorization 
ansatz could be completely misleading for the descrip- 
tion of Df ph. This could be due to the fact that the 
charm quark is too light for the factorization approach 
to be reliable. In fact, since the mass of the Dg lies right 
in the middle of the region populated by highly excited 
light quark resonances, it is possible that the presence of 
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FIG. 29 The missing mass distribution for all D j" — > pn can- 



didates. CLEO-c sees 13 signal candidates. From lAthar et all 
(f2008l ). 



tudes and phases of the different components that con- 
tribute to these final states. It is seen that most three- 
body final states are dominated by pseudo two-body de- 
cays. 

There is an enormous number of appUcations of three- 
body decays of D-mesons. Indeed one of the most im- 
portant ones involves proper determination of branch- 
ing fractions of quasi- two- body decays, such as D ^ pir. 
Also, the possibility of determination of all relative decay 
amplitudes and phases in the Dalitz analysis of D° decays 
allows for novel studies of D'^D mixing and searches of 
CP violation in the charm system. Finally, Dalitz anal- 
yses of £)-decays offer unique ways to study formation 
of light-quark resonances (such as a and n) that are not 
reachable in direct e"'"e~-annihilation experiments. 

In a Dalitz plot analysis the dynamics of a decay 
is investigated by analyzing the kinematic distributions 
by plotting the data such that the event dei isity is 
prop ortional to the matrix element squared ( Dalit j . 
I1953I ). For the three-body decay D abc where a, 
6, ar id c are pseud o-scalars the decay rate can be writ- 
ten (jAmsler et all [2OO81) 



nearby states could significantly affect the decay. In ad- 
dition, the decay happens almost at the threshold for pn 
production, with no large energy release - something that 
factorization-based approaches usually require. Second, 
there could be other decay mechanisms that contribute 
to this transition besides annihilation. For example, in- 
elastic rescattering discussed above could be responsible 
for the bulk of the result. An example of this mecha- 
nism would be a tree-level transition Df — > 77^'^+ with 
subseque nt rescat tering ri^'^n~^ — > pn. It has been ar- 
gued (Ch en et g/.l .r2008) that this mechanism can provide 
a contribution that is consistent with the experimentally- 
measured branching ratio. More work is definitely needed 
for complete theoretical understanding of this and related 
processes. 



1 



32(27r)3Af3 



\Mf dmlf, dml^, 



(62) 



where M is the decay matrix element and raf^ = 

{pi + Pj)^ is the invariant mass squared of particles i and 
j. Note that for Al = constant, the Dalitz plot in vari- 
ables (m^j, m^^) of Eq. ((62|) represents a homogeneously- 
filled shape. Any apparent structures would then repre- 
sent resonant interactions of the final state particles. 



1. Formalism for Dalitz plot fits 

In general, the amplitude for the process D Rc, 
R ^ ah where R is an intermediate resonance and a, 6, 
and c are particles of arbitrary spin, can be written 



IX. DALITZ DECAYS OF D MESONS 

In this Section multibody decays of D mesons are dis- 
cussed. The most extensive studies of multibody decays 
are the Dalitz plot studies performed in three-body de- 
cays. We give an overview of the analysis techniques 
used, and discuss some of the final states that have been 
investigated. A few four-body final states have also been 
investigated and they are discussed next. The last topic 
is the study of inclusive distributions in D decays. 

A. Three-body Dalitz plot analyses 

Many hadronic three-body final states of and 
meson decays have been studied using a Dalitz plot 
analysis in which the resonant substructure has been 
probed. From these analyses we learn about the ampli- 



Mb.{L, niab, ) = J2{ab\Rx)TBimab){cRx\D) (63) 

A 

where L is the spin of resonance i?, and the sum is over 
the helicity states A of R. It is customary to break the 
amplitude of Eq. (j63p into three parts, 

Mr{L, niab, mhc) - ZiL, p, q) (|p|) (|q|) Tfl(m,b), 

(64) 

where Z depends on the spin of resonance R and 
describes the angular distribution of the decay prod- 
ucts. If all final state particles are spin-0, which 
is the case for all of the decays described here (see 
Eq. (HH)), it reduces to Legendre's polynomials. The 
Sl's are the spin-dependent Blatt-Weisskopf penetra- 
tion functions that incorporate effects due to finite- 
size of the final-state hadrons, and Tr is a function 
that describes dynamics of the final-state mesons that 
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incorporate a prescription on how to treat the final- 
state resonances R. The momenta p and q of c and 
a, respectively, are defined in the R rest frame (e.g. 
|q| = - ("^Q + "ib)2)(m^ - (ma - mbY)/2mR. 

The main difference between various analyses of Dalitz 
plots is related to the chosen model for Tr. 

The most common description of Dalitz plots in three- 
body decays is the so called isobar model. In this model 
amplitudes are added coherently for each resonance. A 
nonresonant contribution, which describes a direct decay 
of the D into a 3-body final state, is usually added as 
a coherent contribution uniformly distributed across the 
Dalitz plot, making the total amplitude 

M=M NR ,mab,mbc). (65) 

R 

In the isobar model each resonance is described by a 
Breit-Wigner lineshape. 



TR[mab) = [m\ - TO^fc ^ irriRTabiq)] 



(66) 



Here Tab{q) describes a momentum-dependent width of 
the resonance R, which generalizes narrow- width approx- 
imation, 



rafc(<7) = Yr 



{^T"{^)B'd,M'. (67) 



Resonant fractions, or fit fractions, are defined, for each 
resonance i?, as 



fR 



S\Mr? 



NR 



ErMuI^-^ 



(68) 



where the integration above is over the whole Dalitz plot. 
The sum of fractions, so defined, is not required to be 
unity. One must remember that isobar model is break- 
ing the unitarity and is partly the result of interference 
terms, missing from the denominator, and partly due to 
kinematic limits imposed on the integrals. 

The K -matrix model is used when a proper description 
of a Dalitz plot dominated by broad scalar resonances is 
needed. The if-matrix formalism is, by construction, 
unitary. It follows from a specific parameterization of 
the scattering matrix, 

% = % + 2iT,j = % + 2i{p,y/^ f,f {pfY'^, (69) 

where Tj/ is a Lorentz-invariant scattering amplitude and 
Pi = 2qi/mi are the diagonal elements of the (diagonal) 
phase space matrix. Here qt = rrii^Jl — Ami I * 
breakup momentum for decay channel i. 

The ii'-matrix represents a particular parameterization 
of f. 



T =[I - iKp 



K. 



(70) 



The final-state resonances appear in the iiT-matrix as 
a sum of poles. A particular parameterization of the 



JC-matrix can be chosen, which incorporates data from 
scattering experiments. One useful parameterization of 
thefsT-matrix can be found in ( Anisovich and Sarantsevl 
I2OO3I ). A good descript ion of iiT-matrix formalism can be 
fou nd in (IChung et aZ. . [19951) . See also D. Asner's review 
in (|Amsler et a/.l . l2008t) . 

In addition to the isobar model and the if-matrix 
models presented above, several experiments has use the 
Model-Independent Partial Wave Analysis (MIPWA). 
This approach was firs t used by the E791 collabora- 
tion ( Aitala et al\ . l2006l ). Instead of trying to describe 
the S'-wave as a sum of broad Breit-Wigners resonances, 
which often leads to unitarity violation when they over- 
lap, or using the /^-matrix parameterization this method 
parameterizes the amplitude and phase by dividing the 
7r''"7r~ mass spectrum into descrete slices. The ampli- 
tude a nd phase are interp olated using a Relaxed Cubic 
Spline (jKolbig and Lippsl . Il990t) . 



2. Experimental considerations 

When analyzing data using a Dalitz plot analysis there 
are several experimental effects to consider. The recon- 
struction efficiency for the D candidates is not uniform 
across the Dalitz plot. The momentum spectrum of the 
observed particles will depend on the position in the 
Dalitz plot and affect the efficiency for finding and recon- 
structing the particles. The effect of efficiency variations 
across the Dalitz plot is typically incorporated using a 
Monte Carlo simulation and parameterization of the ef- 
ficiency as a function of the Dalitz plot variables. 

The finite detector resolution is usually neglected as 
the resonances studied are mostly broad compared to the 
detector resolution. There are a few exceptions such as 
(f) K^K^ and lo tt+tt". In these cases the reso- 
lution function has to be convolved with the truth level 
probability distribution. A related effect is resolution ef- 
fects near the phase-space boundary in the Dalitz plot. 
To avoid smearing near the phase-space boundary the fi- 
nal state particles momenta can be recalculated using a 
constraint to the D mass. This forces the phase-space 
boundary to be strictly respected. 

Experimentally we also have to consider backgrounds 
that pass the event selection criteria. The backgrounds 
can be classified into different categories. Combinato- 
rial backgrounds where the selected particles do not all 
come from the decay of a I?. This background may con- 
tain resonances, such as a K* or p. We also have back- 
grounds where all candidates come from a D decay but 
are not signal. These backgrounds include final states 
with identical particles, e.g. ^5""° contributing 

to D° —^ tt+tt^tt" or a decay incorrectly identified as 
a D°, or misidentificd particles such as 7r^7r+7r+ 
reconstructed as I?+ K^tt^tt^. 

In the following Sections different Dalitz plot analy- 
ses will be discussed. As in general it is impossible to 
average the results of different analysis the most recent. 
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FIG. 30 The left plot shows the D" i^^Tr+vr" recon- 
structed mass distribution. The plot to the right shows 
the M^(7r+7r") versus M^(Ji:~7r+ ) Dalitz plot for t he 7,070 
D^' K-TT+Tv° candidates. From lKopp et al\ l|200ll '). 



or precise, results are discussed in more detail for each 
mode. 
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FIG. 31 The D° R-n+Tv" Dalitz fit. From IKopd et al\ 
(|200ll ). 



3. 75° 







The decay Z?" — > ii' 7r+7r° has been stud 
ied by the tagged photo n spectrometer at Fermi 



lab (ISummers eraLl. Il984h. MA RK II I (lAdler et al 



19871) ■ E691 (lAnios et all Il993l) E 687 ( Frabetti et aL . 
1994bl ). and CLEO II (Kopp eiolJ, |200l|). The first of 



these analyses was a simplified Dalitz analysis that did 
not include the interference. The data was fit to an inco- 
herent sum of K~ , J^*°7r°, isT* ~7r+, and nonresonant 
decays. The latest analysis by CLEO II has about a 
factor of 10 higher statistics than any of the earlier mea- 
surements. 

The analysis by CLEO II used 4.7 fb^^ of e+e" colli- 
sion data collected at y^s — 10.6 GeV. The D° candidate 
is required to come from a D*^ D^tt~^ decay. The 
candidate is required to form a _D*+ candidate which 
satisfies 144.9 < M{D*+) - M{D°) < 145.9 MeV. The 
invariant mass distribution of the K~'k^tt^ candidates 
and the 7,070 event selected for the Dalitz plot analy- 
sis are shown in Fig. [301 This sample has a purity of 
96.7 ± 1.1%. The large p+, and K*' resonances 

and their interference is easily seen in this plot. The re- 
sults of the Dalitz plot fit are summarized in Table lXXXIl 
The (0(770)+ resonance dominates the Dalitz plot with a 
fit fraction of about 78.8%. 







4. n 



This decay is of interest for the extraction of the 
CKM an gle 7 in the de c ays B "^ D'^*'^K^ and 
D^*^K'^_ (|Atwood e^"aH . l200lh . When the decay of the 
or in these decays is to a common final states, such 
as -fCgTr+TT", the two decays above interfere and this al- 



low us to measure the CKM angle 7 To extract 7 from 
this analysis a good understanding of the 
Dalitz plot is required. 

This final state has been investigated by many experi 
ment s. The first studies w ere pe rfo rmed by (lAdler et al. 



19871: lAlbrecht et all Il993l: lAnios et all Il995 : 



Frabetti aZ.I . Tl 99211 1994bl ). CLEO was the first exper- 



iment to include doubly Cabibbo suppressed de cays in 
the Dalitz plot analysis ( Muramatsu et a/] . l2002h of this 
decay. They used 10 resonances in their fit: Kgp°, Kgiv, 
ifO/o(980), ZCg/2(1270), ZCg/o(1370), K*{892)-7t+, 
Zs:^(1430)-7r+, ii:2*(1430)-7r+, A'*(1680)~7r+, and the 
Cabibbo suppressed mode A'*(892)+7r^. CLEO found a 
very small fit fraction for the nonresonant contribution 
of (0.9 ± 0-4:'to"lol)%- They also determined that 
the phase difference between the Cabibbo allowed 
Ar*(892)~7r"'" and the doubly Cabibbo suppressed decay 
Ar*(892)+7r~ is consistent with 180° as expected from the 
Cabibbo factors. The significance of the Ar*(892)"'"7r~ 
resonance is 5.5 standard deviations in the study by 
CLEO. 

Bot h BABAR (lAubert et al. 



2005bL l2008bD and 



Belle jAhe et aZ.I . [20081 : iPoluektov ei aZ.l . l2006l ) have stud 



ied this decay with samples well over an order of mag- 
nitude larger than CLEO in t heir program t o determine 
the CKM angle 7. BABAR (|Aubert et all . [2008^ has 
used a data sample of 351 fb~^ collected at the T(4S') 
to study the D'^ AT^tt+tt" Dalitz plot. They recon- 
struct 487,000 D*+ D°n+, AT^tt+tt- decays 
with a purity of 97.7%. The Dalitz plot is fit to a sum of 
eight different P and D wave resonances. They use three 
Cabibbo favored resonances Ar*(892)~, Ar*(1680)~, and 
A'2(1430)~; two doubly Cabibbo suppressed resonances 
a:* (892)+ and A:2*(1430)+, and three CP eigenstates 
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TABLE XXXI Dalitz plot parameters from CLEO II analysis of D° K'ti+-k° (IKopp et aZ ll200ir). 



Mode 



Fit fraction 



Phase (deg) 



p{77Q) + K- 

K*{S92)°n° 
pil700)+K- 

A'o*(1430)-7r+ 
/S'*(1680)-7r+ 
Nonresonant 



0.788 ± 0.019 ± 0.013 ± 0.046 

0.161 ±0.007 iO.OOTt'oSSo^ 
0.127 ± 0.009 ± 0.005 ± 0.015 
0.057 ± 0.008 ± 0.007 ± 0.006 

0.041 ± 0.006 ± O.007^oo35 
0.033 ± 0.006 ± 0.007 ± 0.012 
0.013 ± 0.003 ± 0.003 ± 0.003 
0.075 ± 0.009 ± 0.006: 



-+0.056 
-0.009 



0.0 (fixed) 
163 ±2.3 ±3.1 ±4.3 
-0.2 ±3.3 ±2.2 ±7.0 
171 ± 6 ± 5tl5 
166±5±4.6± 12 
55.5 ± 5.8 ± S.Sttf 



103±8±7±14 
31±4±5.5:! 



;+14 
-3.7 



— 3 
o 

> 




m? (GeV^/c") 



a) 



2 3 

m? (GevV) 




TABLE XXXII Dalitz plot parameters from BABAR anal- 
ysis of D° K'i TT^TT . The errors for the amplitudes and 
phases include only the statistical errors. The fit fractions 
quoted include also the systematic uncertainties. Upper lim- 
its on fit fractions are quoted at 95% confidence level. 



Resonance 


Amplitude 


Phase (deg) 


Fit fraction 


J^*(892)- 


1.740 ±0.010 


139.0 ±0.3 


55.7 ±2.8 


^!:o*(1430)- 


8.2 ±0.7 


153 ±8 


10.2 ± 1.5 


7^2* (1430) - 


1.410 ±0.022 


138.4 ± 1.0 


2.2 ± 1.6 


A'*(1680)~ 


1.46 ±0.10 


-174 ±4 


0.7 ± 1.9 


A"* (892)+ 


0.158 ±0.003 


-42.7 ± 1.2 


0.46 ±0.23 


K^{U30)+ 


0.32 ±0.06 


143 ± 11 


< 0.05 


^^1(1430)+ 


0.091 ±0.016 


85± 11 


< 0.12 


p(770)" 


1 (fixed) 


(fixed) 


21.0 ± 1.6 


c^(782) 


0.0527 ± 0.0007 


126.5 ±0.9 


0.9 ± 1.0 


/2(1270) 


0.606 ± 0.026 


157.4 ±2.2 


0.6 ±0.7 


TTTT S-wave 






11.9 ±2.6 



FIG. 32 BABAR L>° - > K'^-k+tt~ Dafitz plot analysii 
From lAubert et al\ l|2008bl 'l. 



/9(770)°, w(782), and /2(1270). The K matrix formalism 
with the P-vector approximation is used to describe the 
contribution to the amphtude from the tt+tt" 5'-wave. 
The Kn S'-wave includes the K^{USO)- and is:^(1430)+ 
resonances and a nonresonant component. The data and 
the fit projections are shown in Fig. [311 The result of the 
fit is shown in Table |2£2Q£lD 

Belle (|Poluektov et aLl . l2006[ ) has used a 140 fb"! — 



pie collected at the T(4S') to study the D" K^n+i 
Dalitz plot. They select a sample of 104,204 events for 
their analysis with an estimated purity of 96.9%. They 
fit their data to a sum of 15 resonances plus a nonres- 
onant amplitude. The data and projections of their fit 
are shown in Fig. [331 The result of their fit is summa- 
rized in Table IXXXIIll For the two a resonances that 



are included in the fit Belle obtained Mo- 



539 ± 9 



MeV, r,, 



453 ± 16 MeV, = 1048 ± 7 MeV, and 



related with the nonresonant componen t . Bel le has also 
reported a preliminary study ( Abe et a/] . l2008l ) using 605 
fb~^ of data to study this Dalitz plot. 

At this point the uncertainties in 7 are limited by 
statistics. Contributions to the uncertainty on 7 from 
these measurements are not limited by the Dalitz plot 
uncertainty. But with increased statistics the 7 measure- 
ment should improve and a better understanding of the 
Dalitz plot will be needed. At threshold, CLEO-c or in 
the near future BES III, can perform a tagged Dalitz plot 
analysis against CP eigenstates. These analyses will al- 
low a direct determination of the phase required for the 
extraction of 7. This avoids the uncertainties from mod- 
eling of the Dalitz plot. 



5. D" ^ n--K+7v" 

The Dalitz plot of — > 7r~7r+7r° has been studied 
by BABAR as a m eans to extract inform ation about the 
CKM parameter 7 (|Aubert et aLl . l2007cl) similar to what 



was done with D° K^tt^tt . CLEO has also stud- 



Fcra = 109± 11 MeV. The wide (71 resonance is highly cor- ied this decay ( Muramatsu et all . |2002| ). BABAR re- 
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FIG. 34 BABAR D" tt+tt^tt" Dalitz plot analysis. 



FIG. 33 Belle D° ^ JCgTr+yr" Dalitz plot analysis, (a) the 
m^(is:g7r+), (b) the m'^(K'^-K-), and (c) the m^{Tv-Tv+) dis- 
tributions are shown and in (d) the Dalitz plot distribution. 
The points with error bars show the data and th e smooth 
curve is the result of the fit. From iPoluektov et al\ (|2006l ). 



TABLE XXXIII 

of D° ifgyr+TT 
statistical. 



D alitz plot parameters fr om Belle analysis 
(jPoluektov et all 120061 1 . Errors are only 



Resonance 


Amplitude 


Phase (deg) 


if* (892)+ 


i.656 ± 0.012 


137.6 ±0.6 


^"(892)- 


0.149 ±0.007 


325.2 ±2.2 


7^0* (1430)+ 


1.96 ±0.04 


357.3 ± 1.5 


7^0* (1430) - 


0.30 ± 0.05 


128 ±8 


A-2*(1430)+ 


1.32 ±0.03 


313.5 ±1.8 


A-2*(i430)- 


0.21 ±0.03 


281 ±9 


^"(1680)+ 


2.56 ±0.22 


70 ±6 


7<'*(i680r 


1.02 ±0.22 


103 ± 11 


p(770)° 


1 (fixed) 


(fixed) 


i.j(782) 


(33.0 ± 1.3) X 10~^ 


114.3 ±2.3 


/o(980) 


0.405 ± 0.008 


212.9 ±2.3 


/o(1370) 


0.82 ±0.10 


308 ±8 


/2(1270) 


1.35 ±0.06 


352 ±3 




1.66 ±0.11 


218 ±4 




0.31 ±0.05 


236 ± 11 


nonresonant 


6.1 ±0.3 


146 ±3 



constructs 44, 780 ± 250 signal events over a background 
of 830 ± 70 events. The Dalitz plot of these events is 
shown in Fig. [Mj The three p bands are clearly vis- 
ible with a strong destructive interference. BABAR 
used 15 resonances plus a nonresonant contribution to 
fit the data. The resuhs of the fit are summarized in 
Table IXXXIVl The p{770) resonances are clearly the 



strongest features on the Dahtz plot, with fit fractions 
adding to (128.6 ± 1.6)%. The p(1700) resonances con- 
tribute with fit fractions of 3 to 5% each, much smaller 
than the dominant contributions. The remaining am- 
plitudes, including nonresonant, is much smaller. The 
large, destructively interfering, pn amp litudes are sug - 
gestive of an / = dominated final state ( Zemachl . ll965( ). 
This is consistent with the observation that 37r° is 

strongly suppressed. 



6. D° K+K'n'^ 



CLEO 
BABAR lAubert et al. 



(iBesson et all 



2007a[ ) 



[looi) and 
have both studied 
the Dalitz plot of this decay. The BABAR analysis used 
358 fb"-'^ of e~^e~ collision data collected near the T(4iS') 
resonance. A sample with a high purity of about 98.1% 
was selected for this study containing 11,278 ± 110 
£)*+ D^TT~^ tagged candidates. The Dalitz plot and 
the best isobar fit is shown in Fig. [351 The isobar 
model allows for several different solutions that each 
give a similarly good description of the data. At low 
invariant mass an 5-wavc K^K^ contribution 
is needed, but the fit can not distinguish between 
an ao(980) and a /o(980). Similarly, at intermediate 
K^K^ invariant mass either a /2(1525) or an /o with 
a similar mass works. In the study of this Dalitz plot 
the relative amplitude and phase of the amplitudes for 
£>" K*~K+ to K*+K~ can be measured. 

Defining and by 



roe 



0-D''^K* + K- 
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TABLE XXXIV Dalitz plot parameters from BABAR analysis of D° 7r"7r+7r°. 



Resonance Amplitude ratio (%) Phase (deg) Fit fraction 



p+(770) 100 (fixed) (fixed) 67.8 ± 0.0 ± 0.6 

p°(770) 58.8 ±0.6 ±0.2 16.2 ± 0.6 ± 0.4 26.2 ±0.5 ±1.1 

p-(770) 71.4 ±0.8 ±0.3 -2.0 ± 0.6 ± 0.6 34.6 ± 0.8 ± 0.3 

p+(1450) 21 ±6±13 -146 ±18 ±24 0.11 ±0.07 ±0.12 

p°(1450) 33±6±4 10±8± 12 0.30 ± 0.11 ± 0.07 

p"(1450) 82 ±5±4 16±3±3 1.79 ± 0.22 ± 0.12 

p+(1700) 225 ±18 ±14 -17±2±3 4.1 ±0.7 ±0.7 

p°(1700) 251 ±15 ±13 -17±2±2 5.0 ±0.6 ±1.0 

p-(1700) 100 ±11 ±7 -50±3±3 3.2 ± 0.4 ± 0.6 

/o(980) 1.50 ±0.12 ±0.17 -59 ± 5 ± 4 0.25 ± 0.04 ± 0.04 

/o(1370) 6.3 ±0.9 ±0.9 156 ± 9 ± 6 0.37 ± 0.11 ± 0.09 

/o(1500) 5.8 ±0.6 ±0.6 12 ± 9 ± 5 0.39 ± 0.08 ± 0.07 

/o(1710) 11.2 ±1.4 ±1.7 51±8±7 0.31 ± 0.07 ± 0.08 

/2(1270) 104±3±21 -171 ±3 ±4 1.32 ± 0.08 ± 0.10 

cr(400) 6.9 ±0.6 ±1.2 8±4±8 0.82 ± 0.10 ± 0.10 

Nonresonant 57 ± 7 ± 8 -11±4± 2 0.84 ± 0.21 ± 0.12 




From lAubert et al\ (|2008bt ). 



BABAR obtains 

rn = 0.599 ± 0.013(stat) ± O.Oll(syst) 

and 

6d = -35.5° ± 1.9°(stat) ± 2.2°(syst) 
consistent with the earUer CLEO results. 



7. D" 



K+R-R 



a 



Thi s mode h as been studied by BABAR (jAubert et all 
l2005al l2008bt ) as part of an analysis for 7 determina- 
tion. BABAR uses a sample of 69,000 reconstructed 



^ KlK+K- decays. The data, shown in Fig. [Ml was 
fit to an isobar model which includes eight resonances. 
The result of this fit is summarized in Table IXXXVI In 
the fit BABAR floats the mass and width of the (/)(1020). 
The ao(980) resonance has a mass very close to KK 
threshold and decays primarily to rjir and is described 
by a coupled channel Breit-Wigner line shape. The data 
is well described by the fit, BABAR finds a reduced 
of 1.09 for 6,856 degrees of freedom. 
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3540404-009 



TABLE XXXV Dalitz pl ot parameters from B ABAR anal- 
ysis of D" -> K%K+K- (jAubert et all l2008b[ ). Errors are 
only statistical. 

Resonance Amplitude Phase (deg) Fit fraction (%) 

is:sao(980)" 1 55^8 

K%(t>{lQ2Qf 0.227 ±0.005 56.2 ± 1.0 44.9 

is:§/o(1370)° 0.04 ±0.06 2 ± 80 0.1 

iC§/2(1370)° 0.261 ±0.020 9±6 0.3 

A'gao(1450)° 0.65 ±0.09 95 ± 10 12.6 

A'-ao(980)+ 0.562 ±0.015 179 ± 3 16.0 

is:"ao(1450)+ 0.84 ±0.04 97 ± 4 21.8 

iC+ao(1450)- 0.118 ± 0.015 138 ±7 0.7 



8. L»o ^ kIt^tt" 

This decay has been studied using a 9.0 fb^^ data sam- 
ple collected using the CLEO I I. V detector i n e+e " col- 



lisions at the T(4S') resonance (|Rubin et aLl . l2004[) . The 
sample contained 155 D° K^rjir^ candidate events. 
The two large contributions to this decay come from 
i^*(892)0?7 and ao(980)°Xg. The projections of the 
Dalitz plot fit is shown in Fig. [37l Fixing the ampli- 
tude for ao(980)°i4:g to be 1 with a zero phase CLEO 
measured 



^iC* (892)0 j; 
9^if*(892)0r; 

FF(is:*(892)°r;) 
FF(ao(980)°Xg) 



0.249 ± 0.032 ± 0.013 ± 0.018, 
(259± 12±9±6)°, 
0.293 ± 0.062 ±0.029 ±0.019, 
1.19 ±0.09 ±0.20 ±0.16, 



where the errors are statistical, systematic, and model de- 
pendence respectively. For the model dependence CLEO 
considered alternative models where they added addi- 
tional resonances. They considered four different alterna- 
tive fits including: a nonresonant component, i^Q (1430)7y, 
K^{lAm)ri ± a2(1320)Xg, and kt/. The fit probability for 
these different fits were 6.4%, 19.4%, 64.7%, and 49.9% 
respectively. The fit with only two resonances had a prob- 
ability of 0.8%. From these alternative fits CLEO-c de- 
rives a fit fraction of 0.246 ± 0.092 ± 0.024 ± 0.087 for 
any additional components beyond the i4r*(892)°ry and 
ao(980)"7^g. 



9. D+ 



K- 



The decay K 7r+7r+ is one of the largest decays 

of the D+. CLEO-c has measured the branching frac- 
tion to be B{D+ K-'K+iT+) = (9.15 ± 0.10 ± 0.16 ± 
0.07)%. The Dalitz plot for this dec ay has been s tudie d 



1987 ). 



993) 



by sev eral experimen ts MARK III ("Ad ler et al 
NA14 (lAlvarez et aLl . [1991). E691 (A nios et al. 
E687 (|Frabetti et aZ.I . Il994bf ). E791 (lAitala et al\. 
I2006D . and most recently by CLEO-c (jBonvicmi et al. 



2002, 
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FIG. 37 CLEO - 
From lRubin et al\ \2Q04 \. 



KgTj-K Dalitz plot analysis. 



l2008f ). This Dalitz plot is interesting as the only clear 
resonant contribution from iC*°(892) only has a 12% fit 
fractio n and a contributio n of over 60% from Kn S-wave. 
E791 (|Aitala et all |2002| ) obtained a good fit including 
a large low-mass K~Tr'^ scalar resonance k. This fit ob- 
tained fit fractions t hat was significantl y different from 
earlier studies. E791 ( Aitala et all . \200di) re-analyzed the 
data using a model independent partial wave analysis. 
The CLEO-c analysis also uses the same model indepen- 
dent partial wave analysis. 

The CLEO-c study is based on 572 pb~^ of e+e~ col- 
lision data collected at the -0(3770) resonance. The data 
sample selected for the Dalitz plot analysis consists of 
140,793 events with a background of about 1.1%. The 
projections of the Dalitz plot is shown in Fig. [351 The 
CLEO-c analysis finds that in order to get a good de- 
scription of the data, either in the isobar model or using 
the model independent partial wave analysis for the Ktt 
S'-wave, they need to include a / = 2 n'^n'^ S'-wave. 
CLEO-c implements this / = 2 7r+7r+ S'-wave either us- 
ing an analytic form or using a model independent partial 
wave analysis. The model independent partial wave anal- 
ysis results agree with the analytic form and both give a 
good fit. CLEO-c finds a fit fraction of about 10 to 15% 
for the 1 — 2 7r+7r+ S'-wave. 



10. D+ TT + TT+TT- 

The TT+TT+TT" d ecay h as been studied b y 

E687 (|Frabetti et all . Il997l ). E691 (|Anios et all Il989f) . 
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FIG. 38 Projections of the Dalitz plot fit in the CLEO-c 
Dalitz plot analysis of D"*" — > K~-k'^tv^ in a) for m?(K-K) 
(two entries per candidate, b) for m^(7r^7r+). The data are 
shown as points with error bars. The insets on top shows the 
residuals between the data and the as points with error bars. 
The small contributions in the fit from the -R"*(1680) and 
ii'KMSO) resonances are als o shown in the in sets enhanced 
by a factor of 10. From ,Bonvicini et al\ (|2008l l. 
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TAB LE XXXVI Dalitz pl ot parameters from CLEO-c anal- 
ysis (jBonvicini et all bOOTD of D+ 



FIG. 39 CLEQ-c D+ ^ ti+-k'ti+ Dalitz plot analysis. 
From iBonvicini et al\ (120071 ) . 



Resonance Amplitude Phase (deg) Fit fraction (%) 

p(770)7r+ 1 (fixed) fixed 20.0 ± 2.3 ± 0.9 

/o(980)7r+ 1.4 ±0.2 ±0.2 12 ± 10 ± 5 4.1 ±0.9 ±0.3 

/2(1270)7r+ 2.1 ±0.2 ±0.1 -123±6±3 18.2 ± 2.6 ± 0.7 

/o(1370)7r+ 1.3 ± 0.4 ± 0.2 -21 ± 15 ± 14 2.6 ± 1.8 ± 0.6 

/o(1500)7r+ 1.1 ± 0.3 ± 0.2 -44 ± 13 ± 16 3.4 ± 1.0 ± 0.8 

a pole 3.7 ±0.3 ±0.2 -3 ± 4 ± 2 41.8 ± 1.4 ± 2.5 



E791 (lAitala et oil. I2n01al), FOCUSjLink et al\ . \2mAih . 
and CLEO-c ( Bonvicini et al I I2007D . The most recent 
analysis, with the largest data sample, is the CLEO-c 
analysis. The earlier analysis by E791 had reported the 
need to add a (t(500) Breit-Wigner to the tt+tt" S'-wave 
in order to get an acceptable fit. FOCUS analyzed this 
mode using a X-matrix description of the TT^Tr" 5'-wave. 
They obtained an acceptable fit, but did not rule out the 
need for a cr(500). CLEO-c has studied these decays with 
a sample of about 2,600 signal events, excluding the K'g 
events. The nominal fit using the isobar model supports 
the need for a cr7r+ component. The fit to the isobar 
model is shown in Figure [32] and the result from the fit 
is summarized in Table IXXXVIl 



11. Dt 



The Dalitz plot for 13+ K^K^tt^ is of interest as 
it contains the large (pn^ contribution that tra- 

ditionally has been the reference branching fraction for 
Df decays. The decay Df K'^K~tt^ has been stud- 



ied by E687 ( Frabetti et al\ . Il995al ) using a sample of 
701 events. This analysis showed evidence for a large 
Dt /o(980)7r+ contribution. FOCUS ha s also re- 
porte d a preliminary study of th is Dalitz plot ([Malvezzil . 
|2002[ ). Most recently CLEO-c (|Mitchel]| . has re- 

ported preliminary results from their study of the Dalitz 
plot in this decay. 

The CLEO-c analysis uses 586 pb^-'^ of e+e" collision 
data collected at y/s — 4.17 GeV. This sample corre- 
sponds to about 0.57 x 10*' OfD^^ pairs. In this analy- 
sis about 14,400 £)+ K^K^tt'^ candidates are recon- 
structed with a background of about 15%. The invariant 
mass distribution for the K'^K~tt'^ candidates are shown 
in Fig. [401 The Dalitz plot is shown in Fig. [HI Clearly 
visible in this plot are the (j) ^^'^ K*^ resonances. 

The data are fit to an isobar model including the 
/o(980), 0, /o(1370), /o(1710), i^*(892), and i<-o*(1430) 
resonances. CLEO-c finds that all resonances studied by 
E687 are significant, but that in order to obtain a good fit 
they need to add an additional K^K^ resonance. Sev- 
eral resonant, or nonresonant, contributions gives a sim- 
ilar improvement of the fit quality, though the /o(1370) 
gives the best fit and is used in the main result. The 
result of this fit is shown in Fig. [42l A summary of the 
amplitudes and phases extracted from this fit is shown in 
Table IXXXVIIl CLEO-c obtains a reasonably good fit, 
X^/d.o.f — 178/117, using these resonances. 
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FIG. 40 The K^K vr"'" invariant mass for tlie signal candi- 
dates in_the_CLEO^ Dalitz plot analysis of -k^ . 
From lMitcheil (|2009l ). 
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FIG. 42 The CL EO-c Da l itz pl ot fit for Dj 
candidates. From lMitcheil (|2009l '). 
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FIG. 41 The Dalitz plot for D+ ^ _ft:+A""7r+ candidates in 
the C LEO-c analysis of L>+ K+K-'k+. From ,Mitchelll 
(|2009l ). 



12. m 



The decay D+ 



E791 (lAitala aLl.l2001bl) 



'7r+ ha s been studied b y 
FOCUS (lLinkiLsI],[2Q04i), 



and BABAR (|Aubert et all l2008a[) . The BABAR anal- 
ysis selects 13,179 events with a purity of 80%. The in- 



TABLE XXXVII Dalitz plot parameters from CLEO-c anal- 
ysis of D+ ^K-K+'K+. 



Resonance 


Amplitude Phase (deg) Fit fraction (%) 




1 (fixed) 


(fixed) 


47.4 ± 1.5 


A'o(1430)A'+ 1.51 ±0.11 


146 ±8 


3.9 ±0.5 


/o(980)7r+ 


4.72 ±0.18 


157 ±3 


28.2 ± 1.9 


(/)(1020)7r+ 


1.13 ± 0.02 


-8 ±4 


42.2 ± 1.6 


/o(1370)7r+ 


1.15 ±0.09 


53 ±5 


4.3 ±0.6 


/o(1710)^+ 


1.11 ± 0.07 


89 ±5 


3.4 ±0.5 



variant mass distribution of the Z?j 



candi- 



dates is shown in Fig. [43] and the symmetrized Dalitz 
plot distribution is shown in Fig. [44l The symmetrized 
plot shows two entries in the Dalitz plot for each candi- 
date. The analysis by BABAR includes three resonances, 
/2(1270)7r+, p(770)7r+, and p(1450)7r+. In addition to 
these P- and D-wa.ve resonances the MIPWA is used for 
the 7r"'"7r~ iS-wave. This method parameterizes the am- 
plitude and phase by dividing the tt+tt" mass spectrum 
into 29 slices. The results for the amplitudes and phases 
from the fit for the parameterization of the 5- wave clearly 
show the /o(980) resonance. There is also some evidence 
for the /o(1370) and /o(1500). In Table KXXVlIIl the 
summary of the fit is given. The S'-wave parameteriza- 
tion accounts for a fit fraction of (83.0±0.9± 1.9)%. This 
decay also has an important contribution from a spin-2 
resonance, Df /2(1270)7r+. 
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FIG. 43 The n it tt invariant mass for the signal candi- 
dates in the BABAR Dalit z plot analysis of D'^ n'^n'Tr'^ . 
From lAubert et al\ (|2008af ). 



TABLE XXXVIII Dalitz plot parameters from BABAR 
analysis of Tr^Tr^Tr"*". 



Resonance Amplitude Phase (rad) Fit fraction (%) 

/2(1270)7r+ 1 (fixed) (fixed) 10.1 ± 1.5 ± 1.1 

p(770)7r+ 0.19 ±0.02 ±0.12 1.1 ± 0.1 ± 0.2 1.8 ± 0.5 ± 1.0 

p(1450)7r+ 1.2 ±0.3 ±1.0 4.1 ± 0.2 ± 0.5 2.3 ± 0.8 ± 1.7 

g-wave See Ref fAubcrt et al. 2008a') 83.0 ± 0.9 ± 1.9 



TABLE XXXIX Fit fractio ns and phases from the MARK III 
analysis (|Adler et aZ.l.llQQOal ') of the decay D" K'tt+-k'tt+. 



Resonance Fit fraction (%) Phase (rad) 
K'^p" Transverse (g-wave) 0.142 ± 0.016 ± 0.05 -1.39 ±0.09 

A'"ai(1260)+ 0.492 ± 0.024 ± 0.08 

J^i(1270)"7r+ 0.066 ±0.019 ±0.03 0.71 ± 0.25 

T^'^TT+TT- 0.140 ±0.018 ±0.04 3.07 ±0.09 

K-p°TT+ 0.084 ± 0.022 ± 0.04 -0.30 ± 0.13 

Four-body nonresonant 0.242 ± 0.025 ± 0.06 -1.07 ± 0.08 




mlTU^ n) (GeYVc) 



FIG. 44 The Dalitz plot for signal candidates i n the BABAR 



Dalitz plot analysis of 
(|2008al '). 



From lAubert et al\ 



to extract amplitudes for two-body decays. MARK III 
selected a sample of 1, 281±45 £>" K' 'TT^i: TT^ candi- 
dates. The result of the fit to this samples is summarized 
in Table IXXXIXI The largest two-body decay contribut- 
ing to this final state is if (1260)"'" with a fit 
fraction of 0.492 ± 0.024 ± 0.08. The fit gives a fit frac- 
tion of 0.242 ± 0.025 ± 0.06 for nonresonant four-body 
final states, but it is likely that this includes contribu- 
tions from other wide resonances. 

Th e decay D° K^K~ 'it^it~ has been studied b y 
E687 iFrabetti et aLl.ll995(J). E79 1 (jAitala et aLl . ll998bf) . 



and FOCUS (iLink et all l2005d) . The FOCUS study 
used 1,279 ± 48 events. They performed an unbinned 
maximum likelihood fit including 10 resonances. The 
amplitudes are summarized in Table IXLI The domi- 
nant contribution to the decay rate, about 55%, comes 
from decays to intermediate states with an axial vector 
and a pseudo scalar. The second larges contribution, 
about 30%, comes from intermediate states with two vec- 
tors mesons. The remaining contributions are from three 
body decays D VPP and D SPP. 



B. Four-body decays 



X. CONCLUSIONS 



Similar to the three-body decays discussed in the pre- 
vious section the resonant substructure can be studied in 
higher multiplicity final states. A four-body final state 
has a five-dimensional phase space which is hard to visu- 
alize. 

MARK III (|Adler et all Il990a^ studied the decay 



— !■ K TT+TT TT"*". They performed an unbinned max- 
imum likelihood fit in the five-dimensional phase space 



Charm decays remain an exciting field for both the- 
oretical and experimental investigations. In fact, most 
discoveries in heavy fiavor physics in the last five years 
involved charm quarks one way or another. These include 
mixing, new open-charm D^j states, charmonium 
states X, Y, Z states with ordinary and exotic quantum 
numbers, etc. 

In this review, we touched only a part of a vast field 
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TABLE XL Fit fractions and phases from the FOCUS (iLink et all bOOSd ) analysis of the decay D° 



Mode 

Ki(1270)+K- 



Ki{1270)+K 
Ki{1270)+K 
Ki{1270)+K 
Ki{UOO)+K- 
A:* (892)° A" (892)° 
(t>{1020)p{770f 
p{77QfK+K- 
0(lO2O)7r+7r" 
A'*(892)"A:+7r" 
/o(980)7r+7r~ 



Magnitude Phase (°) Fraction (%) 

p{770fK+ 1 (fixed) (fixed) 

isTo (1430)7r+ 0.27 ± 0.08 ± 0.06 354 ± 19 ± 19 

K * (892)°7r+ 0.94 ± 0.16 ± 0.13 12 ± 12 ± 15 



1.18 ±0.19 ±0.09 
0.39 ±0.09 ±0.11 
1.30 ±0.11 ±0.07 
0.33 ±0.12 ±0.16 
0.30 ±0.06 ±0.06 
0.83 ±0.09 ±0.10 
0.91 ±0.13 ±0.05 



259 ± 11 ± 13 
28 ± 13 ± 10 
49 ±11 ±12 
278 ± 26 ± 20 
163 ± 16 ± 15 
234 ± 10 ± 11 
240 ± 11 ± 17 



18±6±3 
2±1±0 
16±4±5 
33 ± 6 ± 4 
22 ±3 ±4 
3±2±1 
29±2±1 
2±2±2 
1±1±0 
11±2±1 
15±3±2 



of charm physics, the hadronic transitions of charmed 
mesons. We did not review many other exciting de- 
velopments in charm physics. For example, a set of 
hadronic resonant states with new and exciting prop- 
erties has been discovered in both open- and hidden- 
charm quark systems, many exciting results were ob- 
tained in theoretical (lattice) computations and exper- 
imental measurements of leptonic and semi-leptonic de- 
cays of charmed mesons, -mixing was discovered 
and used to constrain New Physics at the scales of sev- 
eral TeV (.Golowich et ai. 2007), etc. Also, experimental 
search for CP-violation in charm transitions remains one 
of the primar y ways of probing New Physics in low-energy 
interactions ( Grossman et all [2007). Finally, we did 



not dis cuss inclusive charm decays, lifetimes of charmed 
states ( Bianco et aLl . l200"3tlGabbiam et al. 1. 120041) . as well 
as charmed spectroscopy and decays of charmed baryons. 

Our knowledge of hadronic charm decays has improved 
significantly over the last few years. The _B-factory ex- 
periments, BABAR and Belle, has very large charm data 
samples that has allowed them to do very precise studies, 
including the absolute hadronic branching fractions for 
but and Df mesons. In addition, the unique CLEO- 
c data samples allow detailed studies of Z?°, D~^, and 
decays. In this review we have covered the status of the 
determination of the absolute branching fractions first for 
and mesons. These measurements are dominated 
by results from CLEO-c and BABAR and have statisti- 
cal uncertainties now below ±1% and systematic uncer- 
tainties of about ±1.8%. The determination of the 
branching fractions is dominated by CLEO-c. The pre- 
viously commonly used normalization mode Df — s- (f>TT'^ 
is not used by CLEO-c any more as it is ambiguous at 
the level of precision now obtained by CLEO-c. CLEO- 
c instead quotes partial branching fractions for a range 
of different K^K^ mass ranges around the cj) resonance. 
These partial branching fractions do not try to disen- 
tangle the contributions from the 4> or other resonance 
contributing to the rate. The CLEO-c measurement ob- 
tains a statistical precision of about 4.2% and systematic 
uncertainties of about 3% in the Dt — > K^K~tt~^ mode. 



This result should improve when CLEO-c includes their 
full data sample. In the future BES III should also be 
able to contribute to the determination of the absolute 
hadronic D branching fractions. The larger samples have 
allowed more detailed studies of Cabibbo suppressed D 
and Ds decays. Decays with smaller branching fractions 
have been explored as well as final states with tt*^ and i] 
mesons that traditionally has been harder to reconstruct, 
but thanks the excellent electro-magnetic calorimeters of 
the BABAR, Belle, and CLEO-c are now accessible. Fi- 
nally, a summary of Dalitz decays of D mesons is given. 
Many of the three-body final states have now been an- 
alyzed for their resonant substructure, and also a few 
final states with more than three particles in the final 
state have been studied. These studies show that most 
of the D decays proceed via pseudo two-body decays. A 
few inclusive measurements of D and Ds decays are also 
presented. 

We are confident that charm quarks will bring us new 
and exciting discoveries, particularly in the field of CP vi- 
olation. While the current experimental data sets are not 
large enough to probe Standard Model- level CP violation 
in the charm sector, the precision of asymmetry measure- 
ments is not limited by systematics uncertainties. This 
means that new samples of charm data from Belle, LHCb 
and Super-B factories will allow for even more precise 
measurements of CP-violating observables. In addition, 
multibody channels, which could be less prone to system- 
atic uncertainties, show great potential for searches for 
CP violation. 



Acknowledgments 

We would like to thank David Cinabro and Rob Harr 
for careful reading the manuscript and insightful com- 
ments. A.R. was supported in part by the U.S. National 
Science Foundation under Grant PHY-0757894 and CA- 
REER Award PHY-0846388. A.R. also thanks the Al- 
fred P. Sloan foundation for their support. A. A. P. was 
supported in part by the U.S. National Science Founda- 



52 



tion under CAREER Award PHY-0547794, and by the 
U.S. Department of Energy under Contract DE-FG02- 
96ER41005. 



References 

Aaltoncn, T., et al. (CDF), 2008, Phys. Rev. Lett. 100, 
121802. 

Abachi, S., et al. (HRS), 1988, Phys. Lett. B205, 411. 
Abe, K., et al. (Belle), 2007, eprint hep-ex/0701053. 

Abe, K., et al. (Belle), 2008, eprint 0803.3375. 
Ablikim, M., et al. (EES), 2005, Phys. Lett. B622, 6. 
Ablikim, M., et al. (EES), 2006a, Phys. Rev. Lett. 97, 121801. 
Ablikim, M., et al. (EES), 2006b, Phys. Lett. B641, 145. 
Ablikim, M., et al, 2007, Phys. Rev. D76, 122002. 
Ablikim, M., et al. (EES), 2008, Phys. Lett. B659, 74. 
Abrams, G. S., et al, 1979a, Phys. Rev. Lett. 43, 477. 
Abrams, G. S., et al, 1979b, Phys. Rev. Lett. 43, 481. 
Acosta, D. E., et al (CDF), 2005, Phys. Rev. Lett. 94, 
122001. 

Adam, N. E., et al (CLEO Collaboration), 2005, Phys. Rev. 

Lett. 94, 232002. 
Adams, G. S., et al (CLEO), 2007, Phys. Rev. Lett. 99, 

191805. 

Adler, J., et al (MARK-III), 1987, Phys. Lett. B196, 107. 
Adler, J., et al (MARK-III), 1988, Phys. Rev. Lett. 60, 89. 
Adler, J., et al (Mark-Ill), 1990a, Phys. Rev. Lett. 64, 2615. 
Adler, J., et al (MARK-III), 1990b, Phys. Rev. Lett. 64, 169. 
Aitala, E. M., et al (E791), 1997, Phys. Lett. B404, 187. 
Aitala, E. M., et al (E791), 1998a, Phys. Lett. B421, 405. 
Aitala, E. M., et al (E791), 1998b, Phys. Lett. B423, 185. 
Aitala, E. M., et al (E791), 2001a, Phys. Rev. Lett. 86, 770. 
Aitala, E. M., et al (E791), 2001b, Phys. Rev. Lett. 86, 765. 
Aitala, E. M., et al (E791), 2002, Phys. Rev. Lett. 89, 
121801. 

Aitala, E. M., et al (E791), 2006, Phys. Rev. D73, 032004. 
Akcrib, D. S., et al (CLEO), 1993, Phys. Rev. Lett. 71, 3070. 
Albrecht, H., et al (ARGUS), 1990, Phys. Lett. B241, 278. 
Albrecht, H., et al (ARGUS), 1991, Phys. Lett. B255, 634. 
Albrecht, H., et al (ARGUS), 1993, Phys. Lett. B308, 435. 
Albrecht, H., et al (ARGUS), 1994a, Phys. Lett. B324, 249. 
Albrecht, H., et al (ARGUS), 1994b, Phys. Lett. B340, 125. 
Alexander, ,J. P., et al (CLEO), 1990a, Phys. Rev. Lett. 65, 
1531. 

Alexander, J. P., et al (CLEO), 1990b, Phys. Rev. Lett. 65, 
1184. 

Alexander, J. P., et al (CLEO), 2008, eprint arXiv:0801.0680 

[hep-ex] . 

Alvarez, M. P., et al (NA14/2), 1990, Phys. Lett. B246, 261. 
Alvarez, M. P., et al (NA14/2), 1991, Z. Phys. C50, 11. 
Amslcr, C, et al. (Particle Data Group), 2008, Phys. Lett. 
B667, 1. 

Anisovich, V. V., and A. V. Sarantsev, 2003, Eur. Phys. J. 
A16, 229. 

Anjos, J. C, et al, 1989, Phys. Rev. Lett. 62, 125. 
Anjos, J. C, et al, 1990, Phys. Rev. D41, 2705. 
Anjos, J. C, et al, 1991, Phys. Rev. D44, 3371. 
Anjos, J. C, et al (E691), 1993, Phys. Rev. D48, 56. 
Artuso, M., B. Meadows, and A. A. Petrov, 2008, Ann. Rev. 

Nucl. Part. Sci. 58, 249. 
Artuso, M., et al (CLEO), 1996, Phys. Lett. B378, 364. 
Artuso, M., et al (CLEO), 1998, Phys. Rev. Lett. 80, 3193. 



Artuso, M., et al, 2003, Nucl. Instrum. Mcth. A502, 91. 
Artuso, M., et al (CLEO), 2008, Phys. Rev. D77, 092003. 
Asner, D. M., and W. M. Sun, 2006, Phys. Rev. D73, 034024. 
Asner, D. M., et al (CLEO), 1996, Phys. Rev. D54, 4211. 
Athar, S. B., et al (CLEO), 2008, Phys. Rev. Lett. 100, 
181802. 

Atwood, D., I. Dunietz, and A. Soni, 2001, Phys. Rev. D63, 
036005. 

Atwood, D., and A. A. Petrov, 2005, Phys. Rev. D71, 054032. 
Aubert, B., et al (BABAR), 2005a, Phys. Rev. D72, 052008. 
Aubert, B., et al (BABAR), 2005b, Phys. Rev. Lett. 95, 
121802. 

Aubert, B., et al (BaBar), 2005c, Phys. Rev. D71, 091104. 
Aubert, E., et al (BABAR), 2005d, Phys. Rev. D72, 091101. 
Aubert, B., et al (BABAR), 2006a, Phys. Rev. D74, 011107. 
Aubert, B., et al (BABAR), 2006b, Phys. Rev. Lett. 97, 
221803. 

Aubert, B., et al (BABAR), 2006c, Phys. Rev. D74, 031103. 
Aubert, B., et al (BABAR), 2007a, Phys. Rev. D76, 011102. 
Aubert, B., et al (BABAR), 2007b, Phys. Rev. Lett. 98, 
211802. 

Aubert, E., et al (BaBar), 2007c, Phys. Rev. Lett. 99, 
251801. 

Aubert, B., et al (BABAR), 2008a, eprint 0808.0971. 
Aubert, B., et al (BABAR), 2008b, Phys. Rev. D78, 034023. 
Aubert, B., et al (BABAR), 2008c, Phys. Rev. Lett. 100, 
051802. 

Aubert, J. J., et al (E598), 1974, Phys. Rev. Lett. 33, 1404. 
Augustin, J. E., et al (SLAC-SP-017), 1974, Phys. Rev. Lett. 
33, 1406. 

Augustin, J. E., et al, 1975, Phys. Rev. Lett. 34, 233. 
Bai, J. Z., et al (EES), 1994, Nucl. Instrum. Meth. A344, 
319. 

Bai, J. Z., et al (BES), 1995, Phys. Rev. D52, 3781. 
Bai, J. Z., et al (EES), 2001, Nucl. Instrum. Meth. A458, 
627. 

Bai, J. Z., et al (BES), 2005, Phys. Lett. B605, 63. 
Balcst, R., et al (CLEO), 1994, Phys. Rev. Lett. 72, 2328. 
Ealtrusaitis, R. M., et al (MARK-III), 1986, Phys. Rev. Lett. 
56, 2140. 

Earate, R., et al (ALEPH), 1997, Phys. Lett. B403, 367. 
Barberio, E., and Z. Was, 1994, Comput. Phys. Commun. 79, 
291. 

Eartch, ,1. E., et al (CLEO), 1998, Phys. Rev. Lett. 80. .3919. 
Bauer, M., B. Stech, and M. Wirbcl, 1987, Z. Phys. C34, 103. 
Beane, S. R., P. F. Bedaque, K. Orginos, and M. J. Savage, 

2007, Phys. Rev. D75, 094501. 
Beneke, M., G. Euchalla, M. Neubert, and C. T. Sachrajda, 

1999, Phys. Rev. Lett. 83, 1914. 
Eergmann, S., Y. Grossman, Z. Ligeti, Y. Nir, and A. A. 

Petrov, 2000, Phys. Lett. B486, 418. 
Bernstein, D., et al, 1984, Nucl. Instr. Mcth. A226, 301. 
Besson, . D. (The CLEO), 2009, eprint 0906.2983. 
Besson, D., et al (CLEO), 2006, Phys. Rev. Lett. 96, 092002. 
Bhattacharya, B., and J. L. Rosner, 2008a, eprint 0812.3167. 
Bhattacharya, B., and J. L. Rosner, 2008b, Phys. Rev. D77, 

114020. 

Bianco, S., F. L. Fabbri, D. Benson, and I. Bigi, 2003, Riv. 

Nuovo Cim. 26N7, 1. 
Bigi, I. I. Y., and H. Yamarnoto, 1995, Phys. Lett. B349, 363. 
Bjorken, J. D., and S. L. Glashow, 1964, Phys. Lett. 11, 255. 
Eloch, F., and A. Nordsieck, 1937, Phys. Rev. 52, 54. 
Elok, B., and M. A. Shifman, 1993, Nucl. Phys. B399, 441. 
Eonvicini, . G. (CLEO), 2008, eprint arXiv:0803.0793 [hep- 



53 



cx]. 

Bonvicini, G., et al. (CLEO), 2007, Phys. Rev. D76, 012001. 
Bonvicini, G., et al. (CLEO), 2008, Phys. Rev. D78, 052001. 
Brandenburg, G., et al. (CLEO), 2001, Phys. Rev. Lett. 87, 
071802. 

Buccella, P., M. Lusignoli, G. Miele, A. Pugliese, and P. San- 
torelh, 1995, Phys. Rev. D51, 3478. 

Buccella, P., et al, 1993, Phys. Lett. B302, 319. 
Buras, A. J., J. M. Gerard, and R. Ruckl, 1986, Nucl. Phys. 
B268, 16. 

Butler, F., et al. (CLEO), 1992, Phys. Rev. Lett. 69, 2041. 
Butler, F., et al. (CLEO), 1994, Phys. Lett. B324, 255. 
Cawlfield, C, et al. (CLEO), 2007, Phys. Rev. Lett. 98, 
092002. 

Chau, L.-L., and H.-Y. Cheng, 1989, Phys. Lett. B222, 285. 
Chau, L.-L., and H.-Y. Cheng, 1992, Phys. Lett. B280, 281. 
Chen, A., et al. (CLEO), 1983, Phys. Rev. Lett. 51, 634. 
Chen, C.-H., H.-Y. Cheng, and Y.-K. Hsiao, 2008, Phys. Lett. 
B663, 326. 

Chen, W. Y., et al. (CLEO), 1989, Phys. Lett. B226, 192. 

Cheng, H.-Y., 2003, Eur. Phys. J. C26, 551. 

Chiang, C.-W., Z. Luo, and J. L. Rosner, 2003, Phys. Rev. 

D67, 014001. 
Chung, S. U., et al., 1995, Annalen Phys. 4, 404. 
Cinabro, D., et al. (CLEO), 1994, Phys. Rev. Lett. 72, 1406. 
Coan, T. E., et al. (CLEO), 1998, Phys. Rev. Lett. 80, 1150. 
Collaboration, B. (BESHI), 2009, Nucl. Instrum. Meth. 

A598, 7. 

Cronin-Hennessy, D., et al. (CLEO), 2008, eprint 0801.3418. 
Csorna, S. E., et al. (CLEO), 2002, Phys. Rev. D65, 092001. 
Dai, Y.-S., D.-S. Du, X.-Q. Li, Z.-T. Wei, and B.-S. Zou, 1999, 

Phys. Rev. D60, 014014. 
Dalitz, R. H., 1953, Phil. Mag. 44, 1068. 
Decamp, D., et al. (ALEPH), 1991. Phys. Lett. B266, 218. 
Dobbs, S., et al. (CLEO), 2007, Phys. Rev. D76, 112001. 
Donoghue, J. F., 1986, Phys. Rev. DBS, 1516. 
Dytman, S. A., et al. (CLEO), 2001, Phys. Rev. D64, 111101. 
Dytnian, S. A., et al. (CLEO), 2006, Phys. Rev. D74, 071102. 
Eeg, J. O., S. Fajfer, and J. Zupan, 2001, Phys. Rev. D64, 

034010. 

Fajfer, S., A. Prapotnik, P. Singer, and J. Zupan, 2003, Phys. 

Rev. D68, 094012. 
Falk, A. F., Y. Mr, and A. A. Petrov, 1999, JHEP 12, 019. 
Falk, A. F., and A. A. Petrov, 2000, Phys. Rev. Lett. 85, 252. 
Frabetti, P. L., et al. (E-687), 1992, Phys. Lett. B286, 195. 
Frabetti, P. L., et al. (E687), 1993, Phys. Lett. B313, 253. 
Frabetti, P. L., et al. (E687), 1994a, Phys. Lett. B321, 295. 
Frabetti, P. L., et al. (E687), 1994b, Phys. Lett. B331, 217. 
Frabetti, P. L., et al. (E687), 1994c, Phys. Lett. B340, 254. 
Frabetti, P. L., et al. (E687), 1995a, Phys. Lett. B351, 591. 
Frabetti, P. L., et al. (E687), 1995b, Phys. Lett. B359, 403. 
Frabetti, P. L., et al. (E687), 1995c, Phys. Lett. B354, 486. 
Frabetti, P. L., et al. (E687), 1997, Phys. Lett. B407, 79. 
Gabbiani, F., A. L Onishchenko, and A. A. Petrov, 2004, 

Phys. Rev. D70, 094031. 
Gaillard, M. K., and B. W. Lee, 1974, Phys. Rev. DIO, 897. 
Gaillard, M. K., B. W. Lee, and J. L. Rosner, 1975, Rev. 

Mod. Phys. 47, 277. 
Gao, D.-N., 2007, Phys. Lett. B645, 59. 
Ge, J. Y., et al. (CLEO), 2009, Phys. Rev. D79, 052010. 
Gell-Mann, M., 1964, Phys. Lett. 8, 214. 
Glashow, S. L., J. Iliopoulos, and L. Maiani, 1970, Phys. Rev. 

D2, 1285. 

Goldhaber, G., et al., 1976, Phys. Rev. Lett. 37, 255. 



Golowich, E., ,L Hewett, S. Pakvasa, and A. A. Petrov, 2007, 

Phys. Rev. D76, 095009. 
Golowich, E., and A. A. Petrov, 1998, Phys. Lett. B427, 172. 
Gronau, M., Y. Grossman, and J. L. Rosner, 2001, Phys. Lett. 

B508, 37. 

Gronau, M., O. F. Hernandez, D. London, and J. L. Rosner, 

1994, Phys. Rev. D50, 4529. 
Grossman, Y., A. L. Kagan, and Y. Nir, 2007, Phys. Rev. 

D75, 036008. 

Harrison, e., P. P., and e. Quinn, Helen R. (BABAR), 1998, 
papers from Workshop on Physics at an Asymmetric B 
Factory (BaBar Collaboration Meeting), Rome, Italy, 11-14 
Nov 1996, Princeton, NJ, 17-20 Mar 1997, Orsay, France, 
16-19 Jun 1997 and Pasadena, CA, 22-24 Sep 1997. 

He, Q., et al. (CLEO Collaboration), 2005, Phys. Rev. Lett. 
95, 121801. 

He, Q., et al. (CLEO), 2008, Phys. Rev. Lett. 100, 091801. 
Hinchliffe, I., and T. A. Kaeding, 1996, Phys. Rev. D54, 914. 
Hoshino, K., et al., 1975, Prog. Theor. Phys. 53, 1859. 
lizuka, J., 1966, Prog. Theor. Phys. Suppl. 37, 21. 
Kamal, A. N., and R. C. Verma, 1987, Phys. Rev. D35, 3515. 
Kolbig, K. S., and H. Lipps, 1990, cERN Program Library 
E211. 

Kopp, S., et al. (CLEO), 2001, Phys. Rev. D63, 092001. 
Kubota, Y., et al. (CLEO), 1992, Nucl. Instrum. Meth. A320, 
66. 

Kuraev, E. A., and V. S. Fadin, 1985, Sov. J. Nucl. Phys. 41, 
466. 

Lange, D. J., 2001, Nucl. Instrum. Meth. A462, 152. 
Link, J. M., et al. (FOCUS), 2002, Phys. Lett. B541, 227. 
Link, J. M., et al. (FOCUS), 2003, Phys. Lett. B555, 167. 
Link, J. M., et al. (FOCUS), 2004a, Phys. Lett. B585, 200. 
Link, J. M., et al. (FOCUS), 2004b, Phys. Lett. B6G1, 10. 
Link, J. M., et al. (FOCUS), 2005a. Phys. Lett. B6G7, 59. 
Link, J. M., et al. (FOCUS), 2005b, Phys. Lett. B624, 166. 
Link, J. M., et al. (FOCUS), 2005c, Phys. Lett. B610, 225. 
Link, J. M., et al. (FOCUS), 2008, Phys. Lett. B660, 147. 
Lipkin, H. J., 1980, Phys. Rev. Lett. 44, 710. 
Lipkin, H. J., 1987, Nucl. Phys. B291, 720. 
Malvezzi, S., 2002, AIP Conf. Proc. 549, 569. 
Mitchell, . R. E. (The CLEO), 2009, eprint 0903.1301. 
Muramatsu, H., et al. (CLEO), 2002, Phys. Rev. Lett. 89, 
251802. 

Niu, K., E. Mikumo, and Y. Maeda, 1971, Prog. Theor. Phys. 
46, 1644. 

Okubo, S., 1977, Phys. Rev. D16, 2336. 
Peruzzi, I., et al., 1976, Phys. Rev. Lett. 37, 569. 
Peruzzi, I., et al., 1977, Phys. Rev. Lett. 39, 1301. 
Peterson, D., et al, 2002, Nucl. Instrum. Meth. A478, 142. 
Petrov, A. A., 2004, Phys. Rev. D69, 111901. 
Pham, X. Y., 1980a, Phys. Rev. Lett. 45, 1663. 
Pham, X.-Y., 1980b, Phys. Lett. B94, 231. 
Pham, X.-Y., 1987, Phys. Lett. B193, 331. 
Poluektov, A., et al. (Belle), 2006, Phys. Rev. D73, 112009. 
Richman, J. D., and P. R. Burchat, 1995, Rev. Mod. Phys. 
67, 893. 

Richter-Was, E., 1993, Phys. Lett. B303, 163. 

Rosner, J. L., 1999, Phys. Rev. D60, 114026. 

Rubin, P., et al. (CLEO), 2004, Phys. Rev. Lett. 93, 111801. 

Rubin, P., et al. (CLEO), 2006, Phys. Rev. Lett. 96, 081802. 

Savage, M. J., 1991, Phys. Lett. B257, 414. 

Schindler, R. H., et al., 1981, Phys. Rev. D24, 78. 

Selen, M., et al. (CLEO), 1993, Phys. Rev. Lett. 71, 1973. 

Summers, D. J., et al., 1984, Phys. Rev. Lett. 52, 410. 



54 



Sun, W. M., 2006, Nucl. Instmm. Mcth. A556, 325. 

Tian, X. C, et al. (Belle), 2005, Pliys. Rev. Lett. 95, 231801. 

Verkerke, W., and D. Kirkby, 2003, The roofit toolkit for 
data modeling, for documentation and source code, see 
http: / /roofit. sourceforge.net/, eprint 
arXiv:physics/0306116 [physics. data-an]. 



Zemach, C, 1965, Phys. Rev. 140, B109. 
Zhang, L. M., et al. (BELLE), 2006, Phys. Rev. Lett. 96, 
151801. 

Zweig, G., 1964, cERN-TH-412. 



